NS 


| Tonx #7: 


Not 


NovemMBre 1956 


LE JOURNAL DE PHYSIQUE 


ET 


LE RADIUM 


ONDES DE DIFFUSION ET PRESSION OSMOTIQUE 
Par René LUCAS, 


École Supérieure de Physique et de Chimie Industrielles et Faculté des Sciences de Paris, 


Sommaire, — L'analyse de l’agitation thermique de translation de molécules d’une solution 
peut se faire en introduisant la notion d’onde de diffusion. Ces ondes sont caractérisées par une 
vitesse de phase et une vitesse de groupe qui s’expriment à partir du coefficient D de diffusion 
moléculaire. Les actions mécaniques subies par une paroi semi-perméable faisant obstacle au 
passage d’une catégorie de molécules (pression osmotique x) peut s’exprimer par la pression de 
radiation des ondes de diffusion de cette espèce moléculaire. 

Dans le cas le plus général ou D varie avec la concentration et où des effets thermiques 
apparaissent par le mélange irréversible des corps en présence, les quantités nécessaires à 
l'évaluation de la pression osmotique sont calculées. Aux faibles concentrations on retrouve la 
loi de van’t Hoff et il est possible de prévoir des écarts sérieux à cette loi aux fortes concentrations. 

Dans le cas particulier où D est constant et où il n’y a pas d’effet thermique, la théorie donne 
une expression en bon accord avec la loi expérimentale de Raoult sur les variations de tension 
de vapeur des constituants d’un mélange binaire. La relation obtenue se traduit par une repré- 


sentation graphique intéressante pour les applications à l’osmométrie, 


Introduetion. — Les actions osmotiques jouent 
un rôle essentiel dans les propriétés physiques et 
physico-chimiques des solutions. C’est pourquoi 
la pression osmotique, véritable question-clef, a 
fait l’objet de recherches approfondies de nom- 
breux expérimentateurs et physiciens théoriciens. 
- Depuis van ’tHoff qui proposa sa relation bien 
connue, plusieurs savants tels que Boltzmann, 
Nernst, H. A. Lorentz, Planck, W. Gibbs, Ray- 
leigh, P. Langevin, Porter, Callendar, etc. ont 
consacré des mémoires théoriques à ce problème 
‘sans qu’il soit possible de l’envisager comme complè- 
tement résolu (1). 

….  Nernst, considérant la pression osmotique comme 
le phénomène de base, envisage la diffusion molé- 
culaire comme une conséquence de la pression 
osmotique, établissant une liaison directe entre les 
deux phénomènes. 

_ Après avoir été l’élément essentiel de la théorie 

- des solutions diluées (en faisant usage de la loi 
de van ’tHoff) la pression osmotique a été délaissée 
au profit de grandeurs thermodynamiques utilisant 


(2) Le lecteur trouvera dans le livre de J. Ducraux, 
Traité de chimie physique appliquée à la biologie, Hermann 
et Ge, Paris, 1938 une analyse très intéressante de ces 
mémoires. 


la notion de coefficient d'activité moléculaire. 
Sans nier l'intérêt de ces méthodes, il faut recon- 
naître qu’elles masquent notre ignorance et 
renoncent à l’analyse des mécanismes fondamen- 
taux. 

Le présent travail constitue une tentative pour 
donner une représentation de la pression osmo- 
tique en considérant la diffusion moléculaire 
comme le phénomène de base. Ce point de départ 
est donc à l’opposé de celui de Nernst. Nous intro- 
duirons la notion d’ondes de diffusion pour procé- 
der à l’analyse de l’agitation thermique de trans- 
lation des molécules. Les travaux de H. A. 
Lorentz [1] et ceux de L. Brillouin{[2]relatifs aux 
actions mécaniques exercées sur les obstacles par 
les ondes élastiques mettent en œuvre les pressions 
de radiation. On peut étendre aux ondes de diffu- 
sion les mêmes formes de raisonnements, Il est 
ainsi possible de relier étroitement les notions 
d'ondes de diffusion, de paroi semi-perméable, et 
de pression osmotique. 


.[. — Ondes de diffusion. — Dans une solution 
renfermant deux types de molécules 1 et 2 
nous pouvons analyser l’agitation thermique de 
translation de chaque espèce moléculaire soit en 
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imaginant des états de mouvements individuels 
entièrement incoordonnés (chaos moléculaire) soit 
en tenant compte de l’association des états de 
mouvements par l'effet des interactions. Dans ce 
dernier cas, particulièrement valable pour les états 
condensés, c’est le concept d’ondes qui est retenu 
pour l’analyse de l’agitation thermique. 

Si l’on considère des écarts périodiques de la 
concentration c, (de l'espèce 1 par exemple) autour 
de sa valeur moyenne, la loi fondamentale de la 
diffusion moléculaire 

de: 

DAC — s 
montre que des altérations périodiques de fré- 
quence v de c, suivant des ondes planes, se propa- 
geront à la vitesse de phase V — 2VrDv. Nous 
désignerons par ondes de diffusion de telles ondes. 
Celles-ci possèdent une vitesse de groupe 


pe D ND OV 
» (+) 


Ces ondes üe diffusion sont caractérisées par un 
amortissement considérable, mais l’agitation ther- 
mique procède d’un mécanisme d'entretien qui 
recrée de telles ondes à mesure qu’elles dispa- 
raissent. Vous admettrons, sans pouvoir le démon- 
trer, qu’il existe une relation de phase entre les 
ondes qui disparaissent et celles qui sont créées. 

Un tel processus est en tout point comparable 
à celui qui préside à l’entretien des ondes élas- 
tiques d’agitation thermique dans les solides, ou 
encore à l'entretien du rayonnement isotherme 
dans un milieu doué d’absorption pour les ondes 
lumineuses. 


IT. Actions mécaniques des ondes de diffusion sur 
une paroi semi-perméable. — Dans une enceinte 
fermée par un piston semi-perméable ne laissant 
passer que les molécules d’espèce 2 nous disposons 
sur la face externe du piston une couche liquide 
de molécules 2 et nous recherchons laction 
mécanique exercée par les molécules d’espèce 1 


_Piston 
semi-permeable 


À_ Solution (1+2) 


Fic. 1. 


sur la paroi semi-perméable. Celles-ci ne peuvent 
franchir la paroi qui les réfléchit, en sorte que nous 
pouvons dire que cette paroi semi-perméable, 
«transparente » pour les ondes de diffusion rela- 
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tives aux molécules n° 2 est parfaitement réflé- 
chissante pour les ondes de diffusion relatives aux 
molécules n° 1. 

Dans ces conditions la pression osmotique (en ce 
qui concerne son terme cinétique) peut être consi- 
dérée comme la pression de radiation des ondes de 
diffusion pour lesquelles la paroi semi-perméable 
est réfléchissante. 

La manière la plus sûre d'évaluer les actions de 


ces ondes est d’avoir recours aux propriétés de 


Pinvariant de Boltzmann-Ehrenfest ainsi que Pa 
montré L. Brillouin à diverses reprises. 


D’après Boltzmann un système périodique carac- ! 


térisé par l’énergie Ex pour un état vibratoire 
de fréquence v obéit à la condition d’invariance 
Loin 
V 
batique 
On en déduit 


— çct® au cours d’une transformation adia- 


(1) 


Pour des ondes de diffusion dont les plans d’onde « 


seraient parallèles à la surface du fond et à celle 
du piston distants de / on a 


Ve t 


(4 nombre entier, invariant lors d’une transfor- 
mation adiabatique) d’où 


5 5 à 
Le (2) 


et d’autre part 


dV 
LOF unr 
re 


lorsque nous ferons varier progressivement la dis- 
tance ! en déplaçant le piston, c, représentant la 
concentration des molécules prisonnières. Si S 
est la surface du piston 


Seite Rad où C'ee 
l C1 
D’après (1) et (2) | 
ÔEe __ v DV Ov c,0V de, Ov 
FE, V dv v Voci ci. v 
par suite 
dv ROA RE €, dV\ à 
Cr) TEEN 
Or 
es d'ÉCRRTE 
F'ovr Ù 
et 
es DST 
E, V V \5c,) ol t 
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La quantité ). étant la dérivée prise pour 


une transformation adiabatique sans échange de 
chaleur. 

Si l’on avait considéré des ondes stationnaires 
dont les plans d’onde auraient été orientés perpen- 
diculairement à la surface du piston on aurait 
trouvé, par un raisonnement semblable 


ÔE, _ _Uc/èV\ àl 
HT ae 


. On en déduit que pour des ondes parfaitement 
diffusées suivant toutes les orientations le résultat 


devient 
SE, _ Ur, afôV\ 1N 
E, ne PL I 
Tenant compte de V — 2 VD on en déduit 
M HUL: de, ÿ 
TH SE oi. + qe 


Le facteur . 


que les autres termes du second membre. Il en 
résulte que pour un ensemble d’états vibratoires 
répartis selon un spectre de fréquence on aura 


+35 (el r2Ee 


représente l'énergie d’agitation thermique 


— 2 est indépendant de v ainsi 


SE, — —2f 
ZE, 
v 
de translation des molécules d’espèce 1 


ne 3 R& 
BE = VO Di 
(Ÿ volume de la solution, constante molé- 
culaire des gaz, 9T nombre d’Avogrado, M, masse 
moléculaire, T absolue). 


nee = v: RT 
M; 
d’où 
SE, @Te 3 c (BD) 1 
PLUS 17, SIA (3) 


La pression & agissant sur le piston semi- 
perméable sera, en remarquant que 


dE, = — 5 SÔ 
ARTE: 3 c, [D 
sf M; Re Se A. | A: (a) 


On reconnaît de suite l’expression de la loi de 
van ’t Hoff dans le cas limite des petites concen- 
trations et l’expression précédente que j’ai donnée 
il y a quelques années [3] permet de prévoir des 
écarts à cette loi lorsque l’on aura affaire à des 
concentrations moyennes ou fortes. 

L'étude de la pression & sur la base de cette 
relation (4) présente de sérieuses difficultés. En 
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effet la solution est soumise à une pression hydro- 
staiique égale à s. Lors d’une variation adia- 
batique du volume, la quantité a. devra 
s’exprimer à partir de 
dD 


De is EE 


dD dD 
= d& + sr 17 


. dD ; HET Ps 
où représente la dérivée de la diffusion par 


rapport à la pression p. 
Par suite 


dD\ |_ fèD|* D ds) , dD ar 
ESA Lee dp Gel TT E 
La condition exprimant le caractère adiaba- 
tique du changement s’obtiendra en formant la 
chaleur élémentaire dQ accompagnant une trans- 


formation par variation des paramètres indépen- 
dants et en faisant dQ = 0. 


Cette méthode introduit le paramètre ee sur 


lequel nous n’avons à l'heure actuelle aucune 
donnée expérimentale ; d’autre part on est ainsi 
? 


Te : à Dee : do 
conduit à une équation différentielle en de "2 
C1 


fermant ce paramètre ce qui rend difficilement 
utilisable la relation (4). 

Du point de vue expérimental il est important 
de remarquer que les propriétés des solutions sont 
(de manière quasi constante) étudiées à pression 
constante ou s’écartant peu de la pression atmo- 
sphérique. Il est hors de doute que la pression 
osmotique dépend de la pression hydrostatique 
exercée sur la solution — on peut par exemple 
citer des mélanges binaires dont la zone de démix- 
tion peut être modifiée par action {de la pression. 

Il est donc préférable d’éviter le rôle d’une 
pression hydrostatique variable appliquée à la 


solution pour l’évaluation de la quantité Se 
1/Q 
On peut tourner la difficulté par la méthode sui- 


vante. 


III. Pression osmotique d’une solution main- 
tenue à pression hydrostatique constante. — 
Constituons un récipient cylindrique muni d’une 
surface semi-perméable S solidaire du corps de 
pompe. Au-dessus de S nous disposons la solution 
étudiée et au-dessous de S se trouve un piston 
étanche P mobile. Par déplacement de P vers le 
bas nous pouvons faire passer les molécules n° 2 
de la solution vers la partie inférieure. Les molé- 
cules n° 2 supposées à l’état liquide seront sou- 
mises à une pression hydrostatique négative — x, 
dans un état métastable (absence de formation de 
vapeur) que nous supposerons maintenu pendant 
la transformation élémentaire considérée. C’est 
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cette quantité x qui mesurera la pression osmotique 
de la solution étudiée. 

Nous conviendrons de compter positivement la 
hauteur /, de la solution au-dessus de S dans le 
sens vertical ascendant et nous compterons positi- 


Paroi S 
semi - permeable Î + 


Fic. 2. — La paroi $ semi-perméable 
est solidaire du corps de pompe. 


vement les*déplacements Ô/, du piston dans le 
sens opposé. 

Lors d’un déplacement adiabatique Ô/, du pis- 
ton P (5, > 0) l'opérateur fournit un travail 
mécanique xS'Ô/, qui se retrouve sous deux formes : 
d’une part création (pour les molécules 2) d’une 
énergie potentielle récupérable de dilatation, et 
d'autre part modification de l’énergie thermique 
des ondes de diffusion des molécules 1 dans le 
volume occupé par la solution. 

Pendant le déplacement 37, la quantité Z, subit 
un accroissement Ô/.. 

Si l’on retranche du travail xS3l, la quantité 
NS M 
2 2 1 + TKk2Q 
de dilatation (x2Q compressibilité adiabatique des 
molécules 2) l’utilisation de la relation (3) donne 


12Q 
rs |, RCE Me pr 
Tv NERO) 


_RTe 3 c; PD\ 15, 
ne +55 Goal à 
% désignant le volume de la solution. Or Ÿ = SI, 
d'où 


_ RTe,; 3 c; (D 
‘Aa [: T3D en 


IV. — Calcul des chaleurs de transformation 
accompagnant une variation élémentaire de con- 
centration. — Soit une solution constituée de 
la masse », de molécules 1 et de la masse m, de 
molécules 2. Le volume de cette solution disposée 
dans le récipient muni du fond hémi-perméable 


qui représente l’énergie potentielle 


1+%2QT ( n) 
1 ÿl 
1+ 32%: 7 j 


M + m 5e 
est p, — PT (pe masse spécifique). La concen- 
M P 
tration du corps n° À est ç, = —2©—, 
P 137 mi Ms 
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1) Ajoutons, à température constante, la masse 
Sm, du corps 2 à la solution, c, varie de dc. Nous 
avons 


ANR ne LATTES 
DE nm SPAS ide C1 
d’où 
M1 © Se] = MP ÿ 
ne Ê Mi + M6 dl (me + m2)? ce 


= mare 
s (M + ma)? 
Soit _ la quantité de chaleur à fournir par. 
L d > D j 
unité de masse du mélange et par unité de chan- 
gement de concentration afin de maintenir cons- 
tante la température. Il faudra céder la quantité 
de chaleur 


dc; [1 2%] = 


F dc, OM. 


Ôq — 45.07 me 


(m1 + mi) e, de Mage 


2) Retirons de la solution la même masse Ôm2 
du corps n° 2 en déplaçant vers le bas le piston P, 
l'opération étant faite à température constante. 


te doit être effectué 


É KR no 
Le travail mécanique 
2—7) 

Par-" désignant la masse spécifique de 2 sous pres- 
sions — x. k 

3) Pendant cette opération réversible il a fallu. 
fournir au système une certaine quantité de chaleur 
afin de maintenir la température constante. Soit: 
og 
DC 
masse du mélange et par variation unitaire de. 
concentration du corps 1, nous avons ici (M, “+ mm) 
ÿq' n A 5 
_ (— dc,) pour cette quantité de chaleur fournie. 

1 

4) La mise sous pression — x de la masse Ôm 
à température constante a nécessité la cession au 
système d’une certaine quantité de chaleur. Si A, 
est le coefficient thermodynamique usuel cette 


la quantité de chaleur à fournir ; par unité de. 


quantité sera ha dm,(— x) = dm T (2%) %, 


Ta 


SOit ÔMo (x, coefficient de dilatation ther- 


p2 
mique du corps 2). 


5) On ramène de manière isotherme la masse 
Ôm, de la pression — x à la pression extérieure (de 
valeur supposée négligeable vis-à-vis de x). 

L'opérateur reçoit une certaine énergie due à la 
libération de l’énergie potentielle accumulée, ce 
qui peut s'exprimer en écrivant qu’il faut fournir 
au système l’équivalent changé de signe soit 


— X2 T° ÔM 


2 P2 
(42 compressibilité isotherme). 


ol 


6) La transformation précédente a exigé l’apport 
de la quantité de chaleur Ôm,h,r égale et de signe 
opposé à celle du 4).. 

Le bilan des transformations restituant l’état 
initial est donc 


0 — a _ A : T 
NAT 1 AT) 
(1 — 2 S) (m3 + m9 
Mi P dq” 1 X2 2 
a re 
C1 d Ôc 2 
He) qurnne 
m 
Tenant compte de c, — Dee 
M + Mo 
OL or. 
{ 0 = — — NE il TS © 
: dy j dC1 Ne ep \Dp/TI Lox—m)C1 263201 


_ d’où l’on déduit le coefficient calorifique massique 
de changement unitaire de concentration de la 
solution par voie réversible. 


= —h—21(%) | | 
A dC1 e \oc,/T Pot—n)C1 2201 


V.— Transformation adiabatique élémentaire. 
— Soit C la chaleur spécifique à pression constante 
dela solution Ale coefficient calorifique défini précé- 
demment. Une transformation élémentaire réver- 
sible retirant de la solution à travers la paroi 
semi-perméable la masse dm, et s’accompagnant 
d’un changement ÔT de température exigerait la 
quantité de chaleur 


80 Um, + m2) (CÔT + N5c;) —h, x ms. 


Si la transformation se fait sans échange de 
chaleur 5Q = 0 d’où 


ÔM 
CÔT + AdC: — he mu Em == 0) 
De €, = UP on déduit pour une transfor- 
APE Ne 

mation adiabatique 

= M M de ] 

Be (mm +m) * Mi + M Fs de 

d’où 


ù Ps 
En remplaçant ri par cette valeur il vient 
À 2 


itTm 


à met he Lee CA 
Q 
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VI. — Évaluation de la dérivée adiabatique du 
cœfficient de diffusion. — La solution demeurant à 
pression constante, le coefficient de diffusion molé- 
culaire D dépend de deux paramètres qui sont la 
température T et la concentration c, du corps 1. 

Lors d’une variation élémentaire de ces deux 
paramètres 
DD 
di 


(ee = Gale + (or. (are 


D’après l'expression évaluée en V de (5) 
1/Q 


2D 
aD des +57 AT 


d’où 


on en déduit 


+ rfi 


D dD dD pe \oc, 
(se Fe Fa T Ée Ce; 
a e Ôl 
VII. — Évaluation du rapport SJ. accompa- 
2 


gnant une transformation élémentaire adiaba- 
tique. — Ce rapport peut s’obtenir en exprimant 
que la masse totale des fluides reste constante lors 
de la transformation retirant la solution de la 
masse Ôm, du corps 2 à travers la paroi semi- 
perméable par déplacement du piston P. 

La masse spécifique initiale po, de la solution 
devient par suite de la transformation 


de do 
Po + ST ÔT Lie as 
La masse spécifique du corps n° 2 soumis aux 
valeurs — x de la pression et T + ÔT de la tempé- 
rature devient 


û ù 
rm Pros T dE ÔT.. 
La conservation de la masse totale s'exprime par 
de 


Poli = Pot Oo + (60 Se DIRE de. 3e) (l, + dl). 
1 


dT 
- Par suite, en négligeant les termes du second 
ordre 
1 ONE à Li ds di de dp J 
2 Po [e Se Gror | je 2)] 
: : noce 
En ‘introduisant la quantité SA évaluée au para- 
graphe V 


Re 
dl Po Fa dp 


fe AC + hT Ê . En | à 


dc, Ce or. 


918 
or 
DST Pl | ne 
h5, = à, Fe 5, 1 
car 
ici rcre 
d’où l’on déduit 
"(a 
ôl; qe Po Fe dP 
: 5 Gal 
2 1 M T À 
140 , Fè den ha r | pe \oc/a 
PodCr  PodT C 
; . *#. / D 
L'évaluation des quantités figurant dans ( Se. 
1/Q 
2: ë dl 
ainsi que la connaissance de À et de S permettent 
2 


en définitive de relier la pression x à un ensemble 
de grandeurs accessibles par l’expérience. Parmi 
ces grandeurs la dérivée isotherme du coefficient 


HULISE ONDES 
de diffusion avec la concentration d’une part, 
1 


a dérivée de la chaleur de mélange _ d’autre 
1 
part, jouent un rôle fondamental. 


VIII: — Cas des solutions idéales. — Une cir- 
constance particulièrement favorable se présentera 
pour simplifier les relations précédentes si les 
espèces moléculaires 1 et 2 sont très voisines l’une 
de l’autre — ce serait par exemple le cas de molé- 
cules ne différant l’une de l’autre que par des 
atomes isotopes. 

Dans ce cas p et p°,, sont pratiquement égaux ; 
d’autre part il n’y a pas d’effet thermique sensible 
par mélange irréversible des deux espèces molé- 
ÿq 
dC: 
culaire ne dépend pas de la concentration c, et 


culaires — 0. Le coefficient D de diffusion molé- 


JOURNAL DE PHYSIQUE 


N° 11 


les pressions de vapeur partielles des constituants 
1 et 2 du mélange varient suivant une relation 
linéaire des concentrations moléculaires relatives. 
Cette loi de Raoult n’a été, jusqu'ici, justifiée que 
par l’expérience. Il est facile de constater que la 
loi de Raoult implique une loi de croissance de la 
pression osmotique avec la concentration beaucoup 
plus rapide que ne l’indique la loi de van ’t Hoff. 
Il est donc particulièrement intéressant de voir ce 
que donnent les relations précédentes, dans le 
cas des solutions idéales, et d’en faire compa- 
raison avec la loi de Raoult. 

La quantité À $e simplifie. À des termes du 
second ordre près 


Mr vis Lo Te 
C1 Pom) 26201 
d’où 
T {2 T? 
AC, + ho T = — RER BI 
A : Pam) 26 ÿ 


Introduisant les paramètres thermodynamiques 
usuels 


et remarquant que les propriétés thermodyna- 
miques de la solution sont les mêmes que celles de 
l’espèce moléculaire n° 2, on pourra supprimer 


l'indice 2. La relation donnant = se simplifie 


dl 

Ô 
car 2 0. 
dC; 
PPDRE, 1— xx 

Ôle C2) CRE Te pr 

FETE 

C\pm  P 26 


C1 


/dD 

À c n Cr En remplaçant Fe par sa valeur dans 
tend à s'identifier avec le coefficient d’autodiffusion, dC1/0 
par suite dD = V, laquelle on fait (Se) = 0 et remplaçant — par 

dC; C1 Le P(—r) 
On sait d'autre part que les solutions idéales 1 + 7x HODtIONS 
vérifient bien la relation de Raoult selon laquelle p 
34 0D x x 
__arat Tibet + Ta +25) (HET) fox “ 
M, QT XT XT 


en ayant remplacé %e par À (y rapport des cha- 


leurs spécifiques). Si l’on veut dégager les facteurs 
principaux agissant sur la variation de x avec la 
concentration on peut, en première approximation, 
négliger les termes en x ainsi que le terme 


AT T 


devant l’unité, ce qui conduit à la relation approchée 
(et dans le cas des solutions idéales). 


PT ER 


Lust al bn dé. in dû mé ns À dé es à 


‘4 


Lévis 4H, es ét ot. cn Le © 


so DL 2, di 


No 11 
RTC; 
M; 
T (6) 
121 ATe, 


La pression osmotique se présenterait comme 
une fonction hyperbolique de la concentration et 
doit croître avec celle-ci beaucoup plus rapidement 
que ne l’indique la loi de van ’t Hoff (en effet 


est toujours positif). 


D oT 
En vue d’une vérification graphique commode on 
peut mettre la relation (6) sous la forme suivante : 


M; 3 12D 1 
RE 2 TDi 


DIT PC 1 (7) 


C1 
T 


d’après laquelle (toujours pour les solutions idéales) 


ne 3 ee 
la quantité 7 est fonction linéaire de la concen- 
TT 


tration. 


IX. Pression osmotique déduite des tensions 


de vapeur. — Représentons selon le schéma ci- 


dessous la solution et le solvant (molécules 2) séparés 
par la paroi semi-perméable et en équilibre osmo- 
tique. Dans le voisinage de la paroi le solvant est 
sous pression négative dont la valeur absolue 
maximum représente la pression x définie antérieu- 
rement. Il est facile de constater le caractère 
stable d’un tel équilibre osmotique. 


E semi-perméable 


Fic. 3. — La solution surmonte la paroi semi-perméable S. 


Si l’ensemble est à température uniforme; on 
peut relier la pression de vapeur p, des molécules 2 
au niveau de la solution à la pression de vapeur 
saturante p,° et à la pression osmotique x. Réci- 
proquement, si l’on connaît p, on pourra aisément 
en déduire rx. 

De la définition de la compressibilité x du sol- 
vant [2] on déduit la relation entre sa masse spéci- 
fique b et la pression 

» ôe e 
500: sets FE US ni 
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e, désignant la masse spécifique à pression normale 
du liquide 2. D’autre part 


dp = — pgdz = pog ex? dz 
(g accélération de la pesanteur). 
De 
dp ex? — p5 g dz 


on déduit 
s Li 1 fP ; 
= — D NE XD. — 
Po 8 Æ pe = e - xp) 


P représentant la pression négative du solvant 
au niveau de la paroi semi-perméable (r — — P) 


d’où 
P 
[er] 
0 


En développement limité on aura par suite 


eXT — 1 
P08X ; 


h = 
P0oX8 


1 Lee 
d'ARE RER À r2). 
208 (= < 2 rt) 

Si la vapeur du solvant (2) suit la loi des gaz 
parfaits on pourra relier sa pression p, au niveau 
de la solution et à la pression p,, au niveau de 
cote O par le calcul élémentaire classique 


dPy = — Pr28d2 (pr2Masse spécifique de la vapeur). 
Or 
Ms Pas RT d’où Po 8 dz 
po P22 RT 
par suite 
Ps _ — Mg, _—M,; À 
ue R&T RTE) (= Fe ] 
M 


Pa = Por e— Rés (+ 27) 

La loi de variation des pressions de vapeur, 
selon Raoult est p; — pos (1 — x) relation ‘dans 
laquelle x est la proportion moléculaire du corps 
dissous n° 1 


C1 
M, 


er GX 


M; M, 


ane 


Si le mélange est constitué de masses M, et M, 
très voisines (MW, M) on aura | 


A 


C1 


C1 T C2 


Posant 


Por Te 


Fo 


Caen 0 CT ler 

En identifiant les relations donnant p, d’après 
la valeur de x et par application de la loi de Raoult 
on déduit 


920 
ee M Xr2 
Poell — 2) = Poz Fa = Poe Re (r+ 57) 
d’où 
£ X 2 — ATPo RUE (8) 
RS RC Le 
2 M; Épia ù 


Afin de prendre un exemple d’application numé- 
rique j'ai supposé que les espèces 1 et 2 en pré- 


c. (g par cm) 0 0,2 0,4 

rr (CGS) 0 0,708.108 1,64 .108 
c1 (GGS) 0,32.10—8 0,282.10—8  0,243.10—8 
TR 


Le graphique (fig 4) ci-contre indique les varia- 


2 C À 

tions de 2 en fonction de c;. On remarque de 
TR € 

suite que dans un domaine considérable de varia- 


Fic. 4.— Courbes représentatives de 1 et de 2 d’après la 
TR T 


relation approchée de la théorie, courbe de ZE, En 
C 


. « 2 
abscisses concentrations (gr par cm% de solution). 


tion de la concentration ec; (jusqu’à plus de 

k AC 
80% du maximum de celle-ci) la quantité 1 
. . . . . TR 
varie de manière sensiblement linéaire comme le 


prévoit la relation (7) 
C1 _ M: 


3 


19D 1 

DT PC‘! 

et conduit à attribuer au coefficient de température 
EE 

de la diffusion D3T une valeur voisine de 18.10—3 

par degré C. 

Cette quantité est tout à fait de l’ordre de celles 
des quelques données expérimentales connues et 
s'accorde également avec l’ordre de grandeur déduit 
de la relation de Stokes-Einstein 


PR 
7 Grand 
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sence ont les constantes physiques du benzène 

pris à 200 C et pour lequel 

M, & M, = 78 Po = 0,879 y = 1,44 

x — 0,934 10-10 CGS & -== TMS 
C=.0,43 NA AS AOTE=A IIS 


Dans le tableau suivant, kr désigne la pression 
calculée suivant (8) 


0,6 0,7 0,8 0,85 
3,09.108 ,278.108 6,4.108 8,95.108 
01910 0,164-10—$ 0,125:103 00 0950 

ce qui conduit à 
1 
d — 
DT = Fe (n viscosité). 


Il est remarquable de constater que ce coefli- 
cient varie relativement peu d’un liquide à un 
autre. 


À titre de comparaison, le quotient de la 


1 
pression calculée d’après (8) à la concentration €; 
a été porté sur la même figure. On voit combien ce 
dernier quotient s’écarte d’une représentation 
linéaire. On remarque également que la relation (7) 
pe fait pas prévoir une valeur infinie de la pression 
osmotique au voisinage du maximum de concen- 
tration comme l’exige la relation de Raoult. A 
vrai dire, une valeur infiniment grande pour + ne 
paraît pas physiquement acceptable. L'utilisation 
de la relation (5) plus rigoureuse réduirait les écarts 
aux grandes concentrations. 


TR (C.GS)x10° | 


Ta (Raouit) 


= TT (théorie) 


T,, (Van t'Hoff) 


0 02 0,4 06 08 1 


Fic. 5. — Courbes représentatives de la pression osmo- 
tique suivant la loi de van’t Hoff xy, d’après la loi de 
Raoult rg, et d’après la théorie x : abscisses, concentra- 
tions (g par cm de solution). 


MENO' 1 


La figure 5 représente les variations de pression 
osmotique d’après la loi de van ’t Hoff d’après la 
loi de Raoult, et d’après la relation (6). Il y a 
accord satisfaisant entre les deux dernières valeurs 
calculées dans un large domaine de concentra- 
tion. Ce n’est que pour les valeurs de c, dépas- 
sant 0,75 ‘g/cm$ que la valeur calculée par la 
théorie. des ondes de diffusion s’écarte (par 
valeur plus faible) de celle prévue par la loi de 
Raoult. 


X. Intérêt de la représentation graphique de 
. en fonetion de ce. — Les constatations précé- 


A 


dentes conduisent à examiner, pour diverses 
valeurs expérimentales, l'intérêt présenté par la 
représentation graphique proposée. 


F1c. 6. — Courbes représentatives de 2 (trait plein) et 


de — (pointillé) en fonction de c;. En abscisses con- 


1 
centrations (en g par cm° de solution) de saccharose 
dans l’eau. 


Les figures 6 et 7 montrent les courbes repré- 


; c T Û 
sentatives de- et de = en fonction de c dans le cas 
TT 


de solutions aqueuses de saccharose et de glucose 
pour un domaine considérable de variations de 
concentrations. Ces valeurs sont empruntées à 
Berkeley et Hartley (4). 


On constate que la loi de forme = A + Bec 


est beaucoup mieux adaptée à la représentation 
des résultats que la disposition habituelle qui 


. x T AS . 
consiste à représenter FALn fonction de c. Du point 


de vue de la détermination des masses molécu- 
laires par osmométrie nous pouvons donc insister 
sur les avantages de ce nouveau mode de repré- 
sentation, Les valeurs expérimentales connues 
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avec le plus de précision (solutions plus con- 
centrées) jouent un rôle plus efficace et plus certain 
pour obtenir, à concentration nulle, le point 
d’intersection fournissant la masse moléculaire 
cherchée. 


C, (gr/cm*) 
4 1, (gr/cm 
0 01 02 03 04 05 06 


Fic. 7. — Courbes analogues à celles de la figure 6. 
Corps dissous, glucose en solution aqueuse. 


XI. Conclusions. — La méthode d’étude pro- 
posée pour la pression osmotique conduit à donner 
de celle-ci une représentation nouvelle beaucoup 
plus précise et satisfaisante que l’analogie entre 
les gaz et le corps dissous utilisée par van 
+ Hof. 

Tout en retrouvant la loi de van ’t Hoff aux 
grandes dilutions nous avons montré les raisons 
d’être des écarts importants entre cette loi et celle 
valable pour les solutions concentrées idéales. 

Une comparaison a été faite entre la pression 
osmotique d’après notre théorie (cas des solutions 
idéales) et celle prévue par application de la loi 
de Raoult. Cette comparaison conduit à une concor- 
dance très satisfaisante dans un large domaine de 
concentration (jusqu’à 80 % du maximum de 
celle-ci.) 

La comparaison entre les modes de représen- 


. 7 et 5 -& 
tation des résultats expérimentaux, soit 7 en 


à ASC : 
fonction de c, soit - en fonction de c, montre un 
T 


avantage sérieux en faveur de ce dernier. Dans 
le cas des solutions idéales on peut le recommander 
pour les applications de l’osmométrie à la déter- 
mination des masses moléculaires. Le problème 
des solutions non idéales, plus complexe, peut être 
abordé à partir des résultats contenus dans le 
présent travail. Nous pensons fournir un article 
sur ce sujet dans le cadre de la présente théorie. 


Manuscrit reçu le 20 juillet 1956. 
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Towxes (C. H.) et SaawLow (A. L.), Spectroscopie des 
ondes ultracourtes. (Microwave Spectroscopy.) 1 vol. 
15 x 23 cm, 698 p., McGraw-Hill, New-York, 1955, 
12,50. 

Cet ouvrage admirable à tous égards, est dû à l’un des 
grands maîtres de la spectroscopie microhertzienne, Townes, 
dont les travaux fondamentaux dans ce domaine sont 
universellement connus, et à l’un de ses élèves. Les auteurs 
ont réussi à présenter une extraordinaire richesse de 
matériel‘avec un don pédagogique allié à un souci de 
rigueur digne de tous éloges. 

Il suffira d’indiquer les principales têtes de chapitres pour 
se rendre compte de l’intérêt exceptionnel de ce livre. 
Les spectres de rotation de molécules diatomiques font 
l’objet du chapitre 1. Le chapitre 2 est consacré aux molé- 
cules polyatomiques linéaires. Les chapitres 3 et 4 traitent 
des molécules assimilables à des gyroscopes symétriques et 
asymétriques. Le chapitre 5 est consacré aux spectres 
atomiques. Chapitre 6 : Structure hyperfine quadrupolaire 


dans les molécules. Chapitre 7 : Molécules douées de moment 


cinétique électronique. Chapitre 8: Structure hyperfine 
magnétique dans les spectres moléculaires. Chapitre 9: 
Interprétation des constantes de couplage hyperfin à l’aide 
de la structure moléculaires et des moments nucléaires. 
Chapitre 10 : Effets Stark dans les spectres moléculaires. 
Chapitre 11 : Effets Zeeman dans les spectres moléculaires. 
Chapitre 12 : Spectre de l’ammoniac et mouvements entra- 
vés. Chapitre 13 : Profils et largeurs des raies spectrales. 
Chapitre 14 : Éléments et technique des circuits micro- 
hertziens. Chapitre 15 : Spectrographes microhertziens. 
Chapitre 16 : Ondes millimétriques. Chapitre 17 : Mesure 
et contrôle de fréquence. Chapitre 18 : Utilisation de la 
spectroscopie microhertzienne dans l’analyse chimique. 
Les 7 appendices sont consacrés respectivement aux inten- 
sités des composantes de structure hyperfine et aux énergies 
dues aux interactions quadrupolaires nucléaires, aux inter- 
actions quadrupolaires nucléaires dans les molécules linéaires 
et dans les gyroscopes symétriques, aux coefficients relatifs 
aux niveaux énergétiques d’un gyroscope légèrement asymé- 
trique, aux niveaux énergétiques d’un rotateur rigide, aux 
transitions rotationnelles, aux constantes moléculaires mises 
en jeu dans les spectres microhertziens et, enfin, aux pro- 
priétés des noyaux stables (abondance, masses et moments). 
Une bibliographie fort complète et un tableau de constantes 
fondamentales et des facteurs de conversion terminent cet 
ouvrage singulièrement riche et d’une clarté exceptionnelle. 
Nous recommandons ce livre aux lecteurs français en 
soulignant une fois de plus l’intérêt qu’il présente au double 
point de vue scientifique et didactique. $ 
d T. KHAN. 


Kiaces (Gerhardt), Introduction à la physiques des micro- 
ondes. (Einführung in die Microwellenphysik.) 4 vol. 
x1 + 280 pages, 15 X 23,5 cm. Éditions D. Steinkopff, 
Darmstadt, 1956, 29 DM.) 

Les hyperfréquences trouvent actuellement des appli- 
cations de plus en plus étendues, non seulement dans les 
diverses branches de la physique, mais également en chimie, 
voire en biologie et en médecine. Aussi l'ouvrage présenté 
ici se propose-t-il de familiariser le chercheur avec cette 
technique, hybride de l’optique et de la radioélectricité. 

Les six chapitres du livre comportent les sujets suivants : 
Ch. 1 : Lignes de transmission. Ch. 2 : Guides d’ondes. 
Ch. 3 : Circuits résonnants et cavités. Ch. 4 : Oscillateurs et 
détecteurs. Ch. 5: Appareils et méthodes de mesures. 
Ch. 6: Radiateurs d’ondes (cornets, miroirs, lentilles, 
antennes diélectriques et fentes). Le tout paraît former un 
ensemble assez complet et équilibré : au chapitre 5 par 
exemple, on trouve les lignes de mesure (ligne à sonde et 
ligne à écrasement), les ponts (T, Y et T magique), les 
coupleurs directionnels, les transformateurs d’impédances, 
les déphaseurs, les atténuateurs et les ondemètres, avec 
leurs descriptions et leurs modes d’emploi, puis quelques 
paragraphes sur la mesure des fréquences (étalons de fré- 
quence, spectrographes hyperfréquences, stabilisation de la 
fréquence, horloges moléculaires) et enfin les mesures des 
constantes diélectriques. à 

Une bonne bibliographie termine l’ouvrage, dont on peut 
penser qu'il atteint son but, en évitant à la fois d’être trop 
élémentaire et superficiel, ou au contraire trop théorique et 
trop mathématique. 

M. JesseL. 


Essic (Jean), Douze, notre dix futur. (1 vol., 14 X 22 em, 

172 pages, Dunod, Paris, 1955, 450 F.) 

L'auteur passe en revue les avantages présentés par un’ 
SE de numération à base douze par rapport au système 
actuel. 

Il insiste sur le fait que non seulement serait nécessaire 
une refonte des tables de calculs mais également une trans- 
formation des/systèmes d’unités, qu’il étudie d’ailleurs. 

Examen très fouillé de la question, ce livre nous présente 
à la fois les avantages incontestables et les inconvénients 
de ce système sous une forme qui ne rebute jamais. 

R. Buver. 
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THÉORIE CINÉTIQUE DES PLASMAS HOMOGÈNES FAIBLEMENT IONISÉS 
III. L'opérateur de collision dans le cas du gaz de Lorentz imparfait. 


Par Micuxez BAYET, 


Faculté des Sciences, Toulouse. 


Jean-Loup DELCROIX, 


École Normale Supérieure, Laboratoire de Physique, Paris. 


et JEAN-François DENISSE, 
Observatoire de Paris, Section d’Astrophysique. 


Sommaire. — On poursuit l’étude entreprise précédemment [1], [2], en tenant compte, au 
premier ordre en m/M, des échanges d’énergie entre les électrons et les molécules (modèle du 


« gaz de Lorentz imparfait »). 


Les nouveaux termes ainsi introduits ont une action négligeable sur la relaxation des anisotropies, 
mais déterminent la forme de la partie isotrope de la fonction de distribution des vitesses élec- 
troniques f(v). On définit alors les fonctions propres du nouvel opérateur de collision J(f) relatif 
a une fonction isotrope, et on donne la méthode de calcul de ces fonctions et des valeurs propres 


associées. 


L’application de ces résultats à l’étude de la partie isotrope de la fonction f fait l’objet 


de l’article IV. 


I. Généralités. — Dans déux articles pré- 
cédents, [1], [2], nous avons développé de façon 
approfondie une théorie des plasmas homogènes 
faiblement ionisés, basée sur le modèle dit du gaz 
de Lorentz parfait. Dans ce modèle, on suppose 
que la masse m des électrons est infiniment petite 
devant la masse M des molécules, et, de’ce fait, 
m 
È M° 
On néglige notamment les échanges d'énergie entre 
les électrons et les molécules lors des chocs élas- 
tiques. Dans ces conditions, si un champ électrique 
fournit chaque seconde aux électrons une énergie P, 
le gaz d'électrons va s’échauffer continuellement 
puisqu'on n’a prévu aucun mécanisme de perte 
d'énergie, et il est donc impossible d’atteindre un 
état d'équilibre. Il est certain que ceci est contraire 
à l'expérience. De plus, il résulte de cette théorie 
qu’au bout d’un certain temps, la fonction de 
distribution des vitesses électroniques est pro- 
fondément modifiée de sorte qu’on ne peut plus 
considérer l’action du champ électrique comme 
une perturbation faible et que la méthode de 
résolution de l’équation de Boltzmann que nous 
avons utilisée tombe en défaut. 

_ En fait, les électrons perdent de l’énergie : par 
choc élastique sur les molécules de gaz; par choc 
inélastique avec excitation ou ionisation de ces 
molécules ; par rayonnement électromagnétique. 

Il est certain que c’est le premier de ces trois 
phénomènes qui est le plus important si l’on se 
borne à l’étude des gaz pas trop chauds (10 000 0K 
par exemple). C’est d’ailleurs le moins difficile à 
étudier de façon quantitative ; en effet, l'influence 
des chocs élastiques peut être décrite sans faire 


on néglige tous les termes contenant le rapport 


intervenir d’autre propriété particulière du gaz 
considéré que la masse de ses molécules et la loi de 
force électron-molécule. 

C'est pourquoi, dans ce troisième article, nous 
étudierons de façon approfondie les phénomènes 
dus aux chocs élastiques en négligeant les autres 
causes de perte d'énergie. Cette étude sera faite de 
de façon classique en supposant que l'interaction 
électron-molécule est décrite par une loi de force 
centrale. Le but de cette théorie est essentiellement 
de répondre aux deux questions suivantes : 


1) Les chocs élastiques sont-ils toujours suffisants 
pour permettre un état d’équilibre quelles que 
soient la loi de force considérée et l’intensité du 
champ électrique appliqué ? 

2) Quand il y a équilibre, quelles sont les lois que 
suivent cet équilibre : « Température électronique » 
d'équilibre ? Fréquences de relaxation ? 


IT. Principes généraux de l’étude des chocs élas- 
tiques dans le modèle de Lorentz imparfait. — Le 
m 
M 
non nul (gaz de Lorentz imparfait), l'intégrale des 
chocs J n’a plus les propriétés remarquablement 
simples que nous avions trouvées dans le cas du 
gaz de Lorentz parfait et qui se traduisaient par les 
formules : 


rapport — étant maintenant supposé petit, mais 


J(a) ==) (1) 
J(aCim) re aViCim (2) 
J(aSim) = — 4YySim (3) 


a étant une fonction isotrope quelconque de la 
vitesse. 
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Nous devons maintenant exprimer cette inté- 
grale au moyen de la formule générale : 


an = [FT WF —1mepap deav. (1) 


Le calcul nécessite donc la connaissance de la 


= - 
fonction de distribution r(Ÿ) des vitesses des 
molécules. Nous supposerons que Fest isotrope et 
maxwellienne à la température T. 


r(Ÿ) == »( # }” exp. ( 


_ er) à (5) 


TASER 

Ceci veut dire que le gaz des molécules est pra- 
tiquement en équilibre thermodynamique à la 
température 7, c’est-à-dire que la puissance trans- 
mise des électrons aux molécules, puis des molé- 
cules à une source extérieure de température Test 
relativement faible. Ceci n’impose pas que les 
électrons soient aussi à la température T parce que, 
d’une part, leur nombre est faible; d'autre part, 
l'échange d’énergie entre les électrons et les molé- 
cules se fait mal. 

Il est évident que le calcul approché de J, fait 


m | ‘ 
en supposant que -7 est petit mais non nul, va 


conduire à des termes correctifs en — dans les 


M 
formules (1), (2) et (3) et donnera : 
Ja = 0(ÿ) (6) 
H{aCim) = — an Cim[t + 0%), 1>0 0) 
Ta = Ne [1 1 (5): 1>0. (8) 
0,, 0, 0, étant trois termes correctifs infiniment 
petits avec > voit que pour les fonctions 


sphériques anisotropes la correction est faible en 
valeur relative tandis que pour les fonctions 
isotropes, elle est fondamentale du fait que pour 
ces fonctions, on a 

m 


J{a) = 0 pour U— 


0. 

Il paraît donc raisonnable, et c’est ce que nous 
ferons, de négliger 0, et 0, et de ne tenir compte 
que de 0,. Physiquement, cela revient à admettre 
que la relaxation des anisotropies s’opère de la 
m 
M 
petit, et que la seule modification que nous intro- 
duisons concerne l’évolution de la partie isotrope 
de la fonction de distribution. 

L'expression de J(a) = 0, peut se déduire de la 


[TJ 7) à donnée par 
Chapman et Cowling [3] : 


même façon, que — soit nul ou simplement très 


valeur de l'intégrale 
te] 
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On obtient par dérivation : 
J({a) = ass [av(a + (9) 


En fait, nous donnerons dans une publication 
ultérieure un calcul plus direct de cette formule. 
Dans la parenthèse qui y figure, le premier terme 
existerait seul si les molécules, tout en étant de 
masse finie étaient immobiles avant le choc, le 
deuxième terme tient compte de leur agitation 
thermique. 

On peut maintenant former le système d’équa- 
tions aux dérivées partielles équivalent à l’équation 
de Boltzmann de la même façon que précédemment 
mais tout en tenant compte de la formule (9). Si on 
néglige comme nous l’avons dit les termes correctifs 
tels que 0, et 0,, on voit que les seules équations 
modifiées sont les équations (28), (40), (68) de nos 
articles précédents [1], [2], relatives aux termes 
isotropes qui s’écrivent maintenant : 


dé = J(a) 


boo = J (Bo) 


(10) 
(11) 
Co = Do — 2 © Evo Sin 201 + 26. Fy COS 2 wt + J(c00). (12) 


Pour étudier ces équations aux dérivées partielles 
il semble indiqué d’utiliser encore la méthode qui 
consiste à rechercher les fonctions et les valeurs 
propres de l’opérateur J, lequel n’est cette fois 
défini au moyen de la formule (9) que sur l’espace 
des fonctions isotropes. Une méthode analogue a 
d’ailleurs été également indiquée par Holstein [4] 
à propos de la résolution de l'équation de 
Boltzmann. 

Nous montrerons que dans les cas les plus 
courants ces fonctions propres forment une suite 
discrète constituant une base complète et que les 
valeurs propres correspondantes (— à) sont toutes 
négatives sauf la première d’entre elles À, qui est . 
nulle ; la fonction propre associée à À, est la 
fonction maxwellienne. Par suite, étant donné une 
distribution isotrope quelconque, on pourra 
l’analyser en la décomposant suivant la base des 
fonctions propres, c’est-à-dire en une série com- 
prenant un premier terme maxwellien et une série de 
fonctions non maxwelliennes. À chacun de ces écarts 
à la distribution maxwellienne correspond une loi de 
relaxation en eXt; cette méthode permet donc 
d'étudier l’évolution dans le temps d’une distri- 
bution quelconque et de prévoir en l’absence de 
champ comment elle tend vers la distribution 
maxwellienne. De façon générale, on montrera que 
cette évolution dépend fortement de la forme de 
l’interaction entre les électrons et molécules. 

En présence d’un champ électrique, nous déve- 
lopperons chacune des fonctions &5o, 090; Co 200; 
E5o en série de fonctions propres et nous rempla- 
cerons ainsi chacune des équations (10), (11), (12) 
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par une suite infinie d'équations différentielles 


ordinaires par rapport au temps. La solution de ce 
système décrit complètement le comportement de 
la fonction de distribution des électrons en présence 
de champ électrique. 


III. Étude du nouvel opérateur J; Valeurs 
propres, fonctions propres. — L’équation aux 
valeurs propres du nouvel opérateur de collision 
s’écrit en appelant y les fonctions propres et — À 
les valeurs propres : 


J{y) = — Ày. (13) 
Soit d’après (9) : 
1 d RT dy\l 
Do de av oral + 0. (ti) 


Pour mettre l’équation sous une forme canonique 
avec une variable sans dimension, faisons les 
changements de variable et de paramètre : 


D Du D — (15) 
m 
_— Le 0 16 
À by (16) 
Il vient : 
1 dy 
ÈxA 3 Em scs 2 — 
a [ru (u +) + uu?y (17) 
ou 
d?y 2 d y 
de (hu + or): 


+ (6 rs Qu À Log CAQe +) ÿ= 0. (18) 


On remarque tout d’abord que la forme de cette 
équation ne dépend plus que de la fonction v.(u) ; 


elle ne dépend pas de ; il en résulte que les 


valeurs propres initiales À sont proportionnelles 
UE 

ré 

Ceci dit, pour définir des fonctions propres et des 
valeurs propres, il nous faut imposer à y dans (18) 
des conditions aux limites surabondantes. L’inter- 
valle d'intégration s'étendant de u — 0 à u = co 
qui sont deux points singuliers, la discussion 
dépendra essentiellement de la forme que prendra 
Ja fonction v,(u) pour u — 0 et u — co ; on sait [5] 
que si on admet entre les électrons et les molécules 


4 4 à ; ï 
une loi de force centrale en 7 depuis r = 0 jusqu’à 
r — oo, On trouve une expression de v : 

(19) 


V1 = vou? 


vo étant une constante et l’exposant nr dépendant 

de s selon la relation 

AA Ses 9 
SL 


(20) 
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les valeurs correspondantes remarquables de 7 et s 
sont représentées dans le tableau I : 


TABLEAU I 
s 1 ? 3 


D | — co 


5 + oc 
— 1 0 4 


— 3 


En fait, on ne peut pas admettre que la loi de 
force ait une expression aussi simple dans tout le 
domaine de r = 0 à r = © et il en résulte une 
expression de v, en fonction de w certainement plus 
complexe que celle représentée par la formule (19). 

Pour analyser d’un peu plus près les phénomènes 
de façon générale, on peut énoncer les deux règles 
suivantes : 


I) L’EXPRESSION DE y, AU VOISINAGE DE u = Ü 
est déterminée par la forme de la force d'interaction 
à très grande distance. 

Cette règle tombe en défaut dans le seul cas où la 
force d’interaction est. constamment attractive 
(interaction électrons-ions positifs par exemple) ; 
elle est satisfaite dans tous les autres cas car les 
électrons de faible vitesse ne pénètrent jamais 
dans les régions de potentiel même faiblement 
répulsif. 


IT) L’EXPRESSION DE y, AU VOISINAGE DE U — co 
se détermine par la forme de la force d’interaction 
à très courte distance sauf quand la force d’inter- 
action à grande distance dérive d’un potentiel 
coulombien ou d’un potentiel à variation encore 
plus lente (s < 2). Cette proposition est démontrée 
en appendice. 


Dans le cas de l’interaction à grande distance, la 
force électron-molécule est en général due au dipôle 
induit et correspond à s — 5. On pourrait aussi 
envisager le cas s — 3 (molécules ayant un moment 
permanent). Mais on doit exclure le cas s < 2 (en 
particulier molécules ionisées) ; nous référant au 
tableau de correspondance entre nr et p nous ferons 
l'hypothèse générale au voisinage de u = 0 


V,T vu quand u— 0 avec — 3 < n° < 0, (21) 


Le cas de l’interaction à courte distance est plus 
difficile à étudier de façon générale ; c’est d’ailleurs 
la région où une théorie purement classique est en 
défaut ; on sait toutefois que dans cette région la 
force est à variation rapide ; nous reportant à 
nouveau au tableau de cosrespondance entre n et p, 
nous ferons alors une deuxième hypothèse : 


V, uw quand u — © avec 0 < n° < 1. (22) 
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Nous pouvons maintenant étudier les deux 
points singuliers de l’équation (18). Tout d’abord 
au voisinage de w — 0, nous pouvons écrire (18) 
sous la forme : 


d 
dit PRE OQÿ = 0 (23) 
avec 
jp Shen ve 
u 
D 


Nous pouvons dans Q négliger le terme = 
0 


u” 
devant 6 + 2n' sauf quand n' = 0; ceci est 
d’ailleurs mathématiquement nécessaire car 
quand n' n’est pas entier, ce terme introduit de 
graves complications analytiques. En effet, Q n’est 
alors pas développable en série au voisinage 
de u — 0; remplaçons donc Q par l’expression 
approchée : 


Q = 6 +2n D CL 0: 


pour 


Quand n'—0, nous garderons l'expression 
complète de Q. 

Dans les deux cas, nous sommes amenés à dis- 
cuter une singularité régulière par la méthode de 
Fuchs ; l’équation aux indices est : 


r? + (n° + Ar = 


dont les racines sont : 


Il convient de faire une discussion suivant les 
valeurs de n’. 

Cette discussion étant fastidieuse, nous la 
résumerons dans le tableau IT avec les notations 
suivantes : 

rx plus grande racine de l’équation aux indices, 

ru plus petite racine de l’équation aux indices, 

Yr = Tao + au + ....) intégrale de (23) cor- 
respondant à 71, 

Yu deuxième intégrale de (23) indépendante 
de Yi; on sait que quand r; — ri n ’est pas entier, 
Yu s'obtient à partir de rx de même la façon que yr; 
tandis que Si 7; — rn est entier, on obtient en 
général yrx par une méthode de variation des cons- 
tantes à partir de y, sous la forme 


cf =: — [Pau 
I = YI - 
* : yi 


= Nous désignerons enfin par O(yr) et O(yx) des 
quantités équivalentes à yr et Yr au voisinage 
de u = 0. 


du. (24) 
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TABLEAU II 
n' TI TITI O(yr)  O(yrr) 
= 0 — 1 1 ie 
u 
. 1 
0 >n > —1 0 — (n' +1) 1 ES 
n = —1 0 0 1 Logu 
—(n+1) 1 


—3 <n <—1 —(n +1) 0 


Nous pouvons résumer ces résultats de la façon 
suivante : | 
0 

y est finie pour w = 0 
Yu est infinie pour u = 0; 
— pour —3 <n'<—1 
yr et y sont toutes deux finies pour 4 = 0. 


— pour —1 En 


Pour définir des fonctions propres, nous impo- 
serons donc la condition suivante : 


—1<n <O0 y reste finie pour u = 0 


—3<n <—1 = ky pour u — 0 


k constante choisie « priori pour chaque valeur de n° 


(25) 


Dans les deux cas, la condition homogène choisie M 


détermine y à un facteur près. Il suffira d'imposer 
une deuxième condition homogène pour uw = co 
pour que les conditions soient surabondantes et 
définissent des valeurs propres et des fonctions 
propres. 

Étudions donc maintenant ce qui se passe au 


voisinage de # — co et pour ceci, faisons d’abord 
le changement de variable : u — 7 ce qui ramène 


la singularité à l’origine, l'équation (18) devient : 


dy [d dy 6 d 
dé (Los K L ee Fa Lo DE k 
(26) 
qu'on peut encore écrire : 
5 + PY + Qy = 0 (27) 
et compte tenu de ui ren (22), on a : 
# | 
ro 
6 + 2n” 2u qe 
Q Rd à VG 4" 


Pour les mêmes raisons que précédemment, nous 
remplacerons (27) par une équation approchée en 
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mettant © sous la forme simplifiée 
= 6 + 2n” 


ja error 0 
Quand »#” = 0, nous gardons l’expression com- 
plète de Q. 
Dans les deux cas, P est un infiment grand 


d'ordre = et Q d'ordre - ; d’après les règles géné- 


rales relatives aux singularités des équations diffé- 
rentielles [6], on en conclut que l’une des deux 
intégrales de (26) est développable en série au 
voisinage de t — 0 et que l’autre présente en ce 
‘point une singularité essentielle ; on cherche donc 


la première solution sous la forme 


yr = (ao + at + aot? + ....). 


En portant cette expression dans l’équation (27) 
simplifiée, on obtient une équation déterminante 
pour r ; r étant ainsi choisie on peut calculer les 
coefficients de yr de proche en proche; ensuite Yr 
s'obtient par la formule (24) qui s’écrit ici : 

nt 
e un 


‘È 
ME SAUVE 
4 ES Y$r 


Les résultats obtenus de cette façon sont résumés 
dans le tableau LIT : 


(29) 


TABLEAU III 


1/4 
n r O(yx) O(yrr) 
1 
BR 1 e 15 
n' = 0 - + 38 'ATCARS ms : eu uv” 
é uv + 1Vs” 
O<n<1 n+3 mi BA 6 


Nous pouvons donc les énoncer de la façon 
suivante : : 

Quand w — co, yrest équivalent à une certaine 
puissance de % ; yrr contient un facteur et». 

Or, il est facile de voir que des fonctions y ayant 
l’allure asymptotique de yr ne peuvent être des 
fonctions de distribution réelles ; elles donneraient 
en effet une valeur infinie à la valeur moyenne 
d’une certaine puissance de la vitesse. Vous impo- 
serons donc à la fonction y au voisinage de u — © 
la condition suivante : 


n'— 0 
DRE A 


C2 
y et uw” quand u— co 


(30) 
y et quand u — co, 


* 
* * 


En conclusion, les fonctions propres et les 
valeurs propres sont définies par les conditions (25) 
et (30). Nous allons maintenant énoncer quelques 
règles générales les concernant : 
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4) La fonction ys — e7% est toujours fonction 
propre avec la valeur propre zéro. — En effet, on peut 
voir que 7, est intégrale de l’équation (17) si on 
fait u — 0 et par ailleurs y, satisfait évidemment 
aux conditions (30) (en y faisant u — 0) et (25) 
(en choisissant 4 = 0 quand — 3 < n° < —1). 


2) Toutes les autres valeurs propres sont positives. 
— En effet, on peut ramener l'équation aux valeurs 
propres (17) à une équation de Sturm-Liouville 


en y faisant le changement de fonction 
Jet (31) 


il vient : 


(32) 


2 

_ PS eu a + pu? e—w ÿ = 0. 

La fonction d, — 1 qui correspond à 7, est 
fonction propre avec la valeur propre uo = 0; 
elle ne possède aucun zéro dans l'intervalle de 
variation O ...:..... + oo ; il en résulte d’après 
les propriétés générales des équations de Sturm- 
Liouville [7] que u, est la plus petite des valeurs 
propres ; autrement dit, toutes les autres valeurs 
propres a, sont positives. 


3) Les valeurs propres forment un spectre dis- 
continu. — Pour le voir, nous transformons l’équa- 
tion (32) en une équation de Schrédinger à une 
dimension en faisant un nouveau changement de 
fonction et de variable en posant : 


(33) 


(34) 


1 
je. L'AT it 
0 2 


On obtient après un calcul assez long mais direct 
que nous ne reproduisons pas : 


de de. VIO 0 (35) 
avec 
(ON AE À CLP NE en te 
% eu [le FE e?]| 
PR : : (36) 
(2) 
2 4\e 2 


(35) est une équation de Schrôdinger à une dimen- 
sion (1); V est l’énergie potentielle de la particule 
correspondante. Pour discuter la nature du spectre, 
il nous faut connaître le comportement asympto- 
tique de V. Pour l’obtenir, nous pouvons d’abord 
développer (36), il vient : 


(°) Le changement de variable transformant l’équation 
de Sturm-Liouville en équation de Schüdinger à une dimen- 
sion est indiqué par Morse et Feshbach (loc. cit, p. 730) 
mais l’expression du potentiel que donnent ces auteurs est 
inexacte, ce qu’il nous prient de signaler. 


CE 
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y = | 141% 1 1 Revenons aux fonctions y; 

ee 2 L du u e 

1 d'Logv, 1 f u? eu" y; eu y; ew* du = 0 
3dLog", 1 DETUE Ha 
& du à dE F LE BORN TRE (9) co … UV 
PE dde be si u?ew yiy; du = 0 iÆj. (40) 
Quand v, est de la forme y, = vu”, cette c."q: Lde 


formule devient : 


Dre EAU M Du ] 
Pre 2 un | ra + u 
n 1 
3n + 4 MA ere F 
Lu Cu s Enee, 
n u 


Quand w— co, cette dernière expression est 
équivalente à 


Vo Vo 
—— y ul— u) = —'untr2, 39 
2 EME R (39) 
Or, nous avons supposé qu’au voisinage 


de 4 = co, on avait v, — vu" avec 0 < nr < 1. 

Donc, V tend vers l’infini quand w tend vers 
l'infini ; dans ces conditions, il ne peut y’avoir 
d’états libres de la particule représentée par l’équa- 
tion (35) ; tous les états sont des états liés et le 
spectre des niveaux d'énergie, c’est-à-dire des 
valeurs propres est discontinu, (ce. q. f. d.). On 
remarquera au passage que V tendrait vers zéro et 
que l’on aurait un spectre continu si l’on avait 
n" < — 2 ; ce serait en particulier le cas si la seule 
force d’interaction était une force de Coulomb 
correspondant à n° = —3 ; peut-être en résulte-t-il 
que les gaz de Lorentz électrons-ions ont une rela- 
xation de l’échauffement des électrons différente 
de celle des gaz de Lorentz électrons-molécules 
neutres ; cependant nous laisserons maintenant de 
côté ces cas-là pour rester dans le cadre des hypo- 
thèses (21) et (22). Les valeurs propres u; (ou À) 
forment alors une suite discrète, ainsi que les 
fonctions 7, ; On sait aussi d’après une propriété 
générale des équations de Sturm-Liouville que la 
fonctions % possède 2 zéros dans l'inter- 
valle 0 + co. 


4) Les fonctions propres y sont orthogonales entre 
elles sur l'intervalle O .... co avec la fonction 
poids u? e“, — En effet, on sait que dans l’équa- 
tion de Sturm-Liouville sous sa forme classique (32), 
la fonction poids est le coefficient r de puy et que 
l’on a la relation d’orthogonalité : 


b 
L'ruwau-o is; 
a 
soit ici 


CO 
s à utew Mhdu=0 if 


5) La contribution des fonctions propres non 
maxwelliennes y,(j 0) à la densité électronique est 


nulle. — Il suffit pour le voir de remplacer dans (40) u 


Y par Yo = e* etil vient 


CA 
*É u? y; du = 0 


ou 


pe y; & mu? do = 0 CANCO EL, (41) 


6. Normalisation. — À cause du résultat pré- 


cédent, il est inutile de normaliser les fonctions 


propres ; on peut même remarquer que leur signe 
n’est pas déterminé par les conditions que nous 
leur avons imposées. Toutefois, il peut y avoir 
intérêt à choisir ce signe de façon que l’intégrale 


. co . 
d'énergie correspondante ÆE; — +} ut y; du soit 


positive. 


7) La recherche des fonctions propres y, et des 
valeurs propres }, correspondant à une loi de choc 
déterminée v,(u) est un problème que l’on ne peut 
en général résoudre que de façon numérique. Pour 
terminer, nous indiquerons toutefois un cas parti- 
culier où on peut mener analytiquement le calcul 
jusqu’au bout ; c’est celui où la fréquence v, est 
indépendante de la vitesse : 


V1 = Cste. (42) 


L’équation aux valeurs propres (32) se ramène à 
une équation d’Hermite par le changement de 
fonction © — uw et on obtient finalement pour les 
fonctions propres et les valeurs propres associées 
les expressions : 


eu? 


(23) 


Yi = Hi +1(u) 


1/2 


& — i( “) 
où Hx,.1(u) est le polynome d’'Hermite d’ordre 
2i + 1. 

Les premières fonctions propres ont alors les 
expressions suivantes : 


(44) 


Yo—= 2e“ (45) 
= 8e" (° mn D (46) 
Ya = 32 e—u* (u — 5u? + =) (47) 


L 
Il 
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Appendice : Pour les électrons de grandes vitesses, 
la fréquence de collision v, dépend de la forme 
du potentiel d'interaction à courte distance. 
— Supposons que le potentiel d'interaction o(r) 
puisse être considéré comme étant à courte distance 
pour toutes les valeurs de r inférieures à une 
certaine valeur p,. 


8 IR ë 
Fré. 1. 


Pour un électron de vitesse #, on peut définir 


une distance p, telle que : mo? = œ(p,) et si le 
potentiel est répulsif pour r < p, (fig. A), on peut 
toujours choisir v suffisamment grand pour que p, 
soit très inférieur à p.. 

v, se calcule par la formule : 


CA 
= 2rNo [ (1 — cos x) p dp 


dans laquelle la déviation y est elle-même fonction 
de p et v selon la loi F. 


——— 
(l AIT 2p P° 
| re gro Cire 


Tm) Tayon vecteur “à sommet de la trajectoire 
considéré (plus courte distance d’approche) est 
solution de l’équation : 


ARALIENS Der (B) 


On peut écrire : 
VE 2 No [” (4 — cos x) p dp 
+ 2rNo [” (4 — cos x) p dp + ... 
Po 


+ 2rNo [2 No (1 — cos y) 
Di 


ou 
Va VW + v + v, 

Pour tous les chocs proches qui forment v’(p <p), 
la déviation y est finie et on peut écrire en ordre de 
grandeur 

v 7x Nov pi. 
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Pour les chocs lointains, on a de façon générale : 


(7 a) << Te 
Soit d’après (B) 
mT—P 
ou mieux 
ft œ(p) 
P [1 ie mv? Pl] 


Calculons l’ordre de grandeur de v” 
En combinant (A) et (B), on obtient pour la 
déviation dans un choc quelconque : 


co Th dr 
x = r—2 Fin. P[7m) - 
LB" T a tu 2 
F7 ER (r?—7r2 13 mv 


Or pour un choc lointain, dans tout l’inter- 
valle r,, + co la quantité 


r° rm _Plrm) — Ur) 


est du même ordre de grandeur que 
E(Tm) 
: mv? 
2 


donc petite devant 1 ; on a donc pour y l RAA 
approchée : 


CO Tm dr 09 rih ®(7m) p(r) 
==4ÿ} ROM NS np | DA ER en dr 
x T rm rVrry ;. rm P? (r—rh)#2mvi 
© FE(rm) — se rin r dr 
x m +2 RER | mu? (rè— rn)#2 p? 


ou encore : 


rm o(7m) en | 
1 
x D? MW? Jrm 


o(r) ] 2r dr 
o(rin) ie (- SE 15e 
rm 


De manière générale, on peut voir que l'intégrale 


r RERGESS 
ci-dessus est finie et puisque F est très voisin de 1, 


on en conclut que pour les chocs lointains, y est un 


. 7 . Q A r 
infiniment petit du même ordre que _ ; on peut 


d’ailleurs obtenir une expression plus simple pour y 
quand au delà de r,, le potentiel est de la forme : 


k 


To 


œ(r) = 


car on a alors 


63 


930 


“et en faisant le changement de variable : 


on obtient : 


L'intégrale A(s) ne contient plus r,, ; c’est une 
fonction de s, nulle pour s = 1, positive et croissant 
indéfiniment avec s pour s > 1 ; de plus, elle se 
calcule exactement pour s entier, et l’on obtient les 
valeurs suivantes : 


s AND ES A 
à Gi ees LE LS) 
A(s) (DE DEA cu HS Frs T 


On peut écrire finalement : 


HS CL A(s)\? fæ ; 2 
ÿ 2 No f pdp=rNo(S) [° œ(p)?pdp (C) 


avec 


EE 
o(p) 7 ps 
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on obtient 

pes A(s)A? 1 1 

V 21| “) DL pi 
pour s > 2 
et 


pour s & 2 car l'intégrale (C) est alors divergente. #4 
Dans le cas du potentiel coulombien, cette diver- 
gence est logarithmique [9]. 
Dans le cas des potentiels à variation plus rapide 4 
que le potentiel coulombien (s > 2), on voit que v” 
est fini et peut s’écrire : 


Ne Late) 
V x Nops. FAR 
Da RE à ("> A(s)k ) 
1 Nas A 


Quand v — co, Po —+ 0 et v” < v'; la contri- 


L 
bution des chocs lointains, c’est-à-dire de ceux qui 
dépendent de la portion de la courbe de potentiel 
pour laquelle r > p,, devient donc négligeable par 
rapport aux chocs dont la déviation dépend de la 
forme de la courbe d'interaction à courte distance. 
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ÉTUDE DES PLASMAS EN RÉGIME TRANSITOIRE 


Par J. SALMON. 


Sommaire. — L'étude de l’évolution d’un gaz électronique de Lorentz, lorsqu'on supprime les 
forces extérieures qui maintiennent la température des électrons à une valeur supérieure à celle 
des molécules, conduit, lorsque la fréquence de collision est constante, aux résultats suivants : si 
la densité électronique est uniforme la fonction de distribution des vitesses initialement Maxwel- 
lienne garde cette forme avec une température «instantanée» qui décroît exponentiellement avec 
le temps. Il semble donc possible d'étendre la notion de distribution Maxwellienne à des régimes 
non stationnaires. D’autre part, si la densité électronique n’est pas uniforme, la fonction de distri- 
bution des vitesses prend une forme toujours simple tandis que l’équation de diffusion des élec- 
trons, tout en gardant l’aspect classique, contient un coefficient qui varie avec le temps. 


I. Introduction. — Nous avons cherché à cal- 
culer la fonction de distribution des vitesses des 
particules d’un gaz de Lorentz lors du retour de 
celui-ci à l’équilibre thermique avec les molécules 
du milieu. Dans un plasma ionisé par exemple, la 
coupure du champ électrique, qui maintient les 
électrons à une température plus élevée que celle 
des molécules provoque le refroidissement de ceux- 
ci. Nous allons voir qu’on peut, à l’aide de formules 
assez simples, traduire l’évolution de la fonction 
de distribution des vitesses f et de la densité élec- 
tronique 7 lors de ce phénomène à condition que 
certaines hypothèses soient satisfaites. 

Nous imposons au milieu d’être faiblement ionisé 
c’est-à-dire que les densités nr et n°’ des électrons 
et des molécules sont faibles devant la densité des 
molécules V,. Aussi la fonction de distribution f, 
de celles-ci n’est pas affectée par les chocs dûs aux 
électrons ; son:expression est maxwellienne 


3/2  _ Mat” 
| e 2XT, (4, 


hi = Mr 


V1, My, T; vitesse, masse et température des 
molécules. . 

Le gaz d'électrons est traité comme un gaz de 
Lorentz imparfait ; on ne fait intervenir que les 
collisions électrons-molécules mais lors de celles-ci 
on tient compte des variations en direction et en 
module des vecteurs vitesse. 


Les notations adoptées sont les suivantes: trepré- 
Le > 

sente le temps, r la position de l’électron et v sa 

; RARE A ae ; 

vitesse. 1, 7, X sont les vecteurs unitaires portés 


par les axes x y Z. #, «, @ sont les coordonnées 
> 
sphériques du vecteur v. On pose enfin. 
dr = dx dy dz dC = v? sin x da do. (2) 


Dans une collision y représente le module 
commun des vecteurs vitesse relative et 0 leur 
angle. o(y, 0) désigne la section efficace différentielle 
de choc à propos de laquelle nous faisons une 


hypothèse fondamentale. Nous supposerons en 
effet dans la suite des calculs que o(7, 8) est propor- 
tionnelle à y, ce qui signifie que la fréquence 
de collision est constante. Il est alors important 
de remarquer que la quantité À définie par 
l'intégrale 

NT [4 — cos 0) y sin 0 40 (3) 

0 


ne dépend que de la nature des molécules [1]. 


IT. Résolution de l’équation de Boltzmann. — 
La fonction de distribution f des électrons obéit à 
Péquation de Boltzmann qui s'écrit en l’absence 
de forces extérieures 


df +7 def 
e 


si l’on désigne par 2 ta contribution des chocs 


élastiques et si l’on néglige les chocs inélastiques 
ainsi que la recombinaison. On choisit une solution 
du type 


ne A t) aie all ii U; Jos 


— 
F, et p ne dépendent de v qu’à travers son 


module ». Il est aisé de se rendre compte que 


JÉJPET +7 +R 1pac = 0 (6) 


tandis que n représente la densité électronique si 


[LT FAC Re (7) 


Désignons par s la quantité 


0 
Al EC) 


et reportons (5) dans (4), il vient : 


ne + For 
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def 


Le calcul de est très complexe si l’on veut 


tenir compte de l’agitation des molécules. L’utili- 
sation de la formule de Chapman et Cowling [2] 
jointe à celle de l’équation (3) mène au résultat 


def, nNIX DAT: 06 me ]- " 
53 3 

DR nn) [ra] | | 

TA PAS GET Pen Fo + pe S 


Dans 10 groupons d’une part les termes pos- 
sédant x en facteur et d’autre part les termes 
restants. Il vient l'important système 


LÉ F |] — 0 (11) 


m 


n ES 


> 
chere 
LE Er “dr” dt 


Une telle séparation se justifie car si n est indé- 


N,A à FEAT 
dv. dv 


2 
v M 


Fe. 
ralen U2) . 


— 
pendant de r et de t la fonction de distribution 
doit être isotrope. p est alors nul et seule l’équa- 
tion 11 subsiste. Nous choisirons pour celle-ci une 
solution du type 


me 
F, = dlt)e 2XTW (13) 
laquelle conduit à 
__ _m#t nd m?bdT Nm d 5e 
0O—=e 2xT + SAT De ta de. (1-7) 


dt 
équation qu’on sépare en deux pour obtenir le 
système 


aT m 
dy m ANSE 
œ M0 (: —7) 0: (15) 


Résolvant (14) et reportant dans (15), il vient en 
désignant par T, la température électronique 
initiale 
de 8p2- mo? = 
Fe ren si] + RIT RS 


Cette 


formule représente une distribution 


m 
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maxwellienne où la température au lieu de rester 


fixe varie selon la loi 
JE == Ts + (To — T;) e— ef (17) 


w, étant donné par la formule 
&9 = 2 Mix (18) 


Il peut sembler peu correct d’utiliser le mot 
«température » alors que ce terme implique une 


notion d'équilibre mais du moins doit-on signaler 


l’analogie. 

Remarquons enfin que quelle que soit la forme 
de la fonction p l’énergie cinétique moyenne prend 
une forme simple. : 


En effet 
[frac == (1) 
d’où 


sm? £ 2 [ff Fac =; ÆIT, + (To — Ti) eve]. 


Si l’on consent à adopter l’expression « tempé- 


rature instantanée » déduite de la formule 


moi (t) = : KT(t) 


bw'r= 


on aboutit à l’énoncé suivant. 

« Lors du refroidissement d’un gaz électronique 
de Lorentz par le jeu des collisions élastiques entre 
les électrons et les molécules, la température de 
ceux-ci décroît de manière exponentielle si la fré- 
quence de collision est constante. » 


La loi de décroissance est donnée par la 
formule 17. 
He T; + (To Fe T';) e— ot (17) 
avec 
2m 
O9 — M NN, (18) 


N,se confond avec la fréquence de collision v 
lorsque la section efficace est isotrope ©, est alors 
égal au produit de v par deux fois le rapport des 
masses de l’électron et des molécules en présence. 

Enfin, si la densité électronique n reste cons- 
tante, la fonction de distribution des électrons 
prend une forme maxwellienne. 


mo? 


(19) 


f=nFi=n| 


2R KIT, + (To — Tije—®t 


BR — 
] e 2KÎT;+(To—Tie ot]. 


III. Évolution de la densité électronique. — 
Lorsque n varie il faut chercher son expression à 
l’aide de (12) et de la solution de (11). Dans ce but 
nous séparons (12) en deux parties. 


Tirant p de (20) et reportant dans (21), il vient : 


HAE ESA du. ra Ne 
n n U n L 
r(%) E r [en 0) ee 


7 Nat 


No11 


ts Base RS © Eee ie he SAP 2e PT 


> hi . = 
Re AE rec 0 ve ie > je 


Dr Nr 


No 11. 
Multiplions par dC et mtégrons 


on v? 3e (TS PACS 15) eTCof] An. (22) 


K 
TA SS 

L’équation 22 a la forme d’une équation de 
diffusion mais le coefficient de” diffusion D est 
fonction du temps 


K 
D\t) = mN,x À (F1 + (To Ex T;) ému (23) 
On se ramène à la forme classique en posant 
Oot 
sit) (24) 
1/00 


| car ce changement de variable conduit pour 22 à 


l'équation 


ùr = ET 7 An. 


(25) 


IV. Application de l’équation (22).— L’équation 
(22) entend représenter en particulier l’évolution 
de la densité électronique lors de la coupure d’une 
décharge mais il y a lieu de noter qu’il existe au 
sein du plasma une densité ionique n’ initialement 
égale à n. Les électrons diffusant vers les parois du 
tube à décharge plus vite que les ions n’ va tendre 
à devenir supérieur à n d’où l'apparition d’une 
densité de charge (n° — n)e dont l'effet est de 
retarder la diffusion des électrons. L’équation (22) 
n’est donc valable qu’au début du phénomène. 

On peut, par contre, dans d’autres domaines 
envisager des applications possibles de (22) avec des 
conditions moins restrictives. Par exemple lors de 
l'émission d'électrons par un filament à l’intérieur 
d’une enceinte où se trouve un gaz à la pression P, 
et à la température T,, la densité électronique est 
régie par une équation du type (22) avec en plus un 
terme traduisant la production d’électrons par le 
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filament 
dn KAn 
op à mN;À [F1 


Au cas où le filament est lui-même maintenu à 
une température constante, le terme en eT°%# a 
peu d'intérêt. Il n’en est pas de même en régime 
transitoire ou en régime pulsé. Prenons même le 
cas où une brusque surintensité dans l’alimentation 
du filament provoque la rupture de celui-ci. On 
peut alors se servir simplement de l’équation (22) 
en lui imposant des conditions convenables. 


+ (To — Tije—So] + g. (22 bus) 


Appendice 


Dans leurs travaux sur les plasmas Bayet, 
Delcroix et Denisse [3] ont utilisé, pour résoudre 
l'équation de Boltzmann, une méthode basée sur 
le développement en harmoniques sphériques de la 
fonction de distribution f. Cette méthode s’est 
révélée d’un très grand intérêt car elle permet 
d'obtenir un haut degré d’approximation dans des 
cas très généraux. 

Nous allons montrer qu’il y a identité entre nos 
résultats et ceux obtenus par l’application de cette 
méthode à notre problème particulier. 

On part de la formule mathématique 


(Cr FSC z? 
(1 de PU pe dx > 22?p ! Hp + 1 (x). à (26) 
H2p + {(x) désignant le polynôme de Hermite 


de degré 2p + 1. Posant 


eo ee | ae RONA 
2= 1) e 0 e= (5) 


la fonction de distribution F, donnée par la for- 
mule 16 s’écrit 


æ2 
RON TRE 
Fr ee. (1 + 792 


d’où à l’aide de 26 


5 mr? 
FANS [ m2 l e 2KITi+(Ty—Ti)e Dot] 
27 Ki + (To —T;) e 2%] 
RU EEE TR 
= m Fee e 2KT Ni TEA : RES ( Te (27) 
27 KT; | mo° 2 É | 22np | Hogtis KT) 
/2KT, 
On retrouve précisément le résultat de Bayet, 
Delcroix et Denisse. Manuscrit reçu le 15 février 1956. 
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LA DISSYMÉTRIE DU TENSEUR DES EFFORTS ET SES CONSÉQUENCES 


Par Yves LE CORRE, 


Laboratoire de Physique Théorique du Collège de France, Paris et École Supérieure 
de Physique et Chimie, Paris. 


Sommaire. — Nous dénonçons les deux postulats suivants sur lesquels repose la théorie classique 
de l’élasticité : f, 

19 Le système des forces appliquées dans un solide en équilibre, de part-et d'autre d’une surface 
élémentaire, est réductible à une composante unique, le couple résultant étant nul. 

20 Dans un solide en équilibre chaque parallélépipède élémentaire peut être considéré comme 
séparé du reste et soumis comme tel aux principes de la mécanique. 

Nous montrons que le tenseur des efforts peut être dissymétrique. Le nombre des constantes 
élastiques d’un cristal triclinique s’élève donc à 45. Nous analysons certaines conséquences nouvelles 
qui en résultent, elles sont toutes conformes au principe de symétrie de Curie. En particulier, 
il a été montré que, contrairement à la théorie de Voigt, aucune des oscillations transversales qui 
se propagent suivant l’axe binaire d’un cristal de phosphate monoammonique n’a la même vitesse 


que la vibration transversale dégénérée qui se propage suivant l’axe quaternaire. 


La première théorie de l’élasticité cristalline est 
due à Cauchy [1] qui a défini les 21 coefficients 
d’élasticité estimés jusqu'ici suffisants pour rendre 
compte de toutes les déformations conformes à la 
loi de Hooke. Voigt [2] a renouvelé et approfondi la 
théorie de Cauchy par le calcul tensoriel dont il est 
l'inventeur ; Max Born [3] a réédifié l’élasticité 
cristalline sur la structure atomique. Il a rejoint 
toutes les conclusions de Voigt mais il suppose que 
les forces de rappel entre les atomes sont centrales. 
Or on sait que cette hypothèse n’est pas justifiée. 
Au contraire, J. Laval [4] définit simplement les 
forces de rappel entre deux atomes au moyen d’un 
tenseur dérivant des diverses énergies mutuelles, 
biatomiques, triatomiques, tétraatomiques, etc. 
règnant dans le milieu cristallin. Il trouve 
de la sorte 45 coefficients d’élasticité distincts, 
IMyes SYmMétriques seulement en «y et 65 tandis 
que les 21 coefficients de Cauchy-Voigt C,,8s le 
sont de plus en «et y et en Bet à. 

Nous avons montré, dans notre thèse [5], qu’il 
est possible de retrouver les résultats de Laval par 
une théorie macroscopique. Il suffit de se rendre 
compte que le tenseur des efforts T,, n’est pas 
nécessairement symétrique. Dans une série de notes 
récentes, Sir G. V. Raman et Viswanathan [6], [7], 
après avoir repris le point de vue de Laval [4] 
admettent à leur tour que le tenseur des efforts 
n’est pas symétrique. Mais les raisons qu'ils invo- 
quent à ce sujet sont insuffisantes. Nous nous pro- 
posons donc de reprendre, en les complétant, les 
arguments que nous avons avancés auparavant. 


I. La démonstration classique de la relation 
Toy = Tya eSt faisse. — La démonstration classique 
de la symétrie du tenseur des efforts est bien 
connue (fig. 1). On considère un parallélépipède 
élémentaire dont les arêtes sont parallèles au trièdre 


trirectangle OË, £, £,. On représente par des flèches 
la résultante des efforts appliqués par unité de 
surface à chacune des six faces du parallélépipède. 
Les composantes 7,, de la densité tensorielle des 
efforts sont définies en décomposant les forces 
appliquées à la surface unitaire normale à Oë, 
parallèlement aux axes OË,,. 


Fic. 1. — Représentation classique du tenseur des efforts, 


Le calcul de la condition d'équilibre paraît, 
simple : la résultante du système de forces est 
évidemment nulle, le couple ne s’annule que si 7. 
est symétrique en « et y. Mais ce calcul est inutile ; 
il est même sans fondement car la figure sur laquelle 
on raisonne est incomplète. En effet, les tensions 
auxquelles sont soumises chacune des six faces sont 
toutes représentées par des vecteurs polaires (du 
type gradient). Or, il est bien connu qu’un tenseur 
antisymétrique du second ordre doit être représenté 
par un vecteur axtal (du type rotationnel). C’est-à- 
dire qu’un tenseur du second ordre doit en général 
être représenté par trois vecteurs polaires à six 
composantes indépendantes et un vecteur axial à 
trois composantes indépendantes. La repré- 
sentation des tensions par trois vecteurs polaires 


À 
# 
: 
4 

ë 
4 
& 
ë 
we 
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Not 


conduit nécessairement à la symétrie de 7,,. Dans 
le cas de la densité tensorielle des efforts, le vecteur 
axial représente une densité de couple en volume. 
Cette grandeur n’est pas nouvelle : on crée une 
densité de couple en volume dans un corps magné- 
tique (ou diélectrique) anisotrope à l’aide d’un 
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champ magnétique (ou électrique) et dans une : 


L Jame absorbante à l’aide d’une vibration lumineuse 


polarisée circulairement. 

Ainsi la démonstration classique de la symétrie 
du tenseur des efforts est fausse parce qu'elle 
s'appuie sur un postulat que nous récusons mais 
que l’on aurait pu encore énoncer ainsi: « Le 
système des forces appliquées dans un solide en 
équilibre de part et d’autre d’une surface élé- 
mentaire est réductible à une pure unique : 


le couple résultant est nul. » 


Ayant exposé l’erreur théorique de la démons- 
tration classique nous allons donner quelques 
exemples où le tenseur des efforts peut être dis- 
symétrique. 


IT. Possibilité d’un tenseur d’efforts dis- 
symétrique en élasticité statique. — Considérons 


- d’abord un cylindre dont la base est collée sur un 


support rigide et dont la face latérale est soumise 
à un effort tangent uniforme normal à l’axe de 


F1G. 2. — La base supérieure du cylindre est collée sur un 
support rigide, la face latérale est soumise à un effort 
tangent. La réaction de la base doit compenser le couple 
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révolution. Le système des forces appliquées ainsi 
est réductible à un couple. La réaction de la base 
rigide doit annuler ce couple. Cette compensation 
peut se faire d’une infinité de manières. Deux cas 
extrêmes sont représentés sur la figure 2 : la surface 
élémentaire de la base est soumise soit à un 
couple ®, soit à un effort tangent ®. I] semble 
évident que les forces atomiques seront distribuées 
dans la réalité selon un mode intermédiaire, chaque 
surface élémentaire étant soumise à la fois à un 
couple et à un effort tangent. 

Pourtant la théorie classique nie la possibilité 
d’une densité de couple superficielle. [1 est exact 
que si le champ de force interatomique est central 
il en est ainsi : la structure atomique se déforme de 
façon à annuler tout couple en. surfacé ou en 
volume. Mais si nous voulons édifier une théorie de 
l’élasticité sans hypothèse aucune sur le champ de 
forces, nous n’avons pas le droit d’écarter a priort 
la possibilité d’une densité de couple, à moins que 
les conséquences qui en résultent soient absurdes. 
Il se trouve justement que certaines conséquences 
prêtent à discussion. 

Appliquons à un cristal appartenant àl’hémiédrie 
centrée cubique (1), taillé parallèlement à ses axes 
quaternaires, une succession d’efforts homogènes, 
d’abord T;,, puis 7, (de même intensité), 
ensuite 7, — 7, etenfin 7%} + 7,,. La figure 3 
et le tableau I rendent compte des efforts appliqués 


Ta, di re 
nas {fr 8 mil | 2e Ti] [a Tll 
Ta ie Die 
a b e d 


F1c. 3. — Parallélépipède soumis à des efforts T',, ou T;; ; 


résultant ; deux cas extrêmes sont envisagés : chaque les renseignements complémentaires sont donnés dans 
surface élémentaire de la base est soumise : a) à un le tableau I. Les cercles représentent les densités de 
couple normal, b) à un effort tangent. couple normales au plan de la figure. 
TABLEAU I 
a b G d 
Te = T ie 0 Te ET 1 — T 
Toy = 0 Ta =T Tu = —T Tai = 7 
lo = Syosio T 2 = Syooi TL 190 = (S1212 7 Sup21) D tie = (Samio + Sio1) T 
toy = Simor T' loy — Soyso1 ly2 — (S12,21 — So1s21 la = Gage S21521) T 
“ : 12 
W — g‘iva Ie W — 3 °2121 TE W — 3 Sins + Sos21 — 281001) T? 14 = 3 (1912 + Sons21 + 251921) 1? 


ainsi que des déformations et des énergies élas- 
tiques résultantes en fonction des compliances $,+.88 
(de même symétrie que les coefficients 9,5). 
L'expression de la densité d'énergie élastique 


1 NES 
>> Say. RS Toy Tes 


JET ER 
2 «v,p5 


impose des énergies différentes dans les cas a et b 


et une énergie non nulle dans le cas c. Ceci provient 
de ce que les tensions 7, et T,, travaillent sur des 
glissements 4 et { différents. La seule critique 
que l’on puisse opposer à ces figures est de nier la 


(*) Cet exemple et les suivants portent sur des cristaux 
aussi simples que possible, choisis d’après le tableau des 
constantes élastiques affectées aux 32 classes cristallines 
établis dans la référence 8. 
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possibilité de tels efforts. Effectivement, ces efforts 
n’ont aucune existence si les forces interatomiques 
sont centrales, mais rien d’autre ne nous autorise 
à les récuser. 

Il est plus délicat de discuter le cas où le même 
cristal est soumis, non plus à des tensions mais à des 
déformations de, dons 1e — ay €t enfin #o + fn 
(fig. 4, tableau IT). On constate avec surprise que 
les énergies des déformations a et b sont différentes 
alors qu’elles correspondent à la même déformation 
apparente du cristal. De même l’énergie de la défor- 
mation c n’est pas nulle alors qu’elle ne correspond 


PHYSIQUE Ne 11, 
à aucune déformation apparente du cristal. L’expli- 
cation est la suivante : dès que les forces ne sont 
pas centrales, la forme extérieure du cristal déformé 


[] (1 OS 


FiG. 4. — Parallélépipède soumis à des déformations t;2 Ou 
t»1 ; les renseignements complémentaires sont donnés, 
dans le tableau II. 


TABLEAU II 


a b (o d | 
RENE ti = 0 lis, 14 bin EE 
ta — 0 oi Et EEE loi 
Tio = Mina t Ts = Iipoit Die = (Moro — yo o) t Tia = (Oise + Misos) t 
Ta = AMisat Tai = Aopait Tai = (Miosor — Tossos) t Ta — Pas + Io1sa1) t 

1 dL 1 | 
W- = PMRTIe PO W => Mont HR 2 (Mos1o + Monsor — 225 21) À es 5 Pass + Mosso1 +209 21) À 
d’une façon homogène ne suffit pas pour définir dire au tenseur T,, = T, T,;, = — T ; cette défor- 


l’énergie de déformation élastique ; il faut tenir 
compte des déformations de densités électroniques 
et des déplacements relatifs des sous-réseaux 
atomiques, lesquels sont différents dans les cas a 
et b car les efforts exercés sont différents. De même, 
dans le cas c on peut réaliser un déplacement 
relatif des sous réseaux et une déformation des 
densités électroniques sans modifier la forme appa- 


7 à 


ie 
PE 
Da 
d'a 


F1G. 5. — Déformation atomique du quartz soumis à une 
densité de couple en volume parallèle à l’axe optique. 
Les deux sortes d’atomes d'oxygène sont représentées 
par des carrés et des losanges. Les niveaux relatifs sont 
indiqués (blanc, croix, noir). Cette déformation est 
compatible avec une polarisation de l’axe optique qui 
est normal au plan de la figure. 
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rente du solide. Dans notre thèse [5], nous avons 
proposé la figure 5 pour représenter la déformation 
d’un cristal de quartz soumis à une densité 
uniforme de couple suivant l’axe ternaire, c’est-à- 


mation est compatible avec une polarisation piézo . 
électrique de l’axe ternaire interdite par la théorie h 


de Voigt. 


Dans la théorie atomique on considère l’atome « 
comme ponctuel et la théorie de l’élasticité revient , 


à étudier les modifications d’un réseau de points 
en équilibre stable sous l'effet des déformations. 
Dire que les forces interatomiques sont centrales, 
c’est dire que les forces de cohésion ne dépendent 
que de la distance entre les atomes. Il suffit alors 
de connaître les modifications de distances entre 
ces points pour évaluer l’énergie des déformations. 
Si l’on admet que les forces entre les atomes ne sont 
pas centrales, cette connaissance est insuffisante, 
il faut se donner dans un repère quelconque les 
neuf composantes t,, du tenseur des déformations 
et non plus les six composantes t1,, to, loss las + tgo, 
ten + is €t fo + tx Le tenseur des défor- 
mations f, étant un tenseur mixte, il est connu que 
l’on peut, par une rotation des axes, le rendre 
symétrique. Mais les composantes du tenseur des 
efforts et les constantes élastiques sont modifiées à 
leur tour dans ce changement d’axe. Étant donné 
que la matrice de changement de repère dépend 
elle-même des composantes t,,, une telle rotation 
ne fait généralement que compliquer inutilement 
les calculs sans modifier évidemment la valeur de 
la densité d’énergie 


1 
= > Moy,e5 tey 188. 


2 oy.8 

III. Expression correcte du couple auquel est 
soumis un parallélipipède élémentaire. — Il faut 
noter qu’il n’est pas nécessaire d'admettre la notion 
de densité de couple en volume pour montrer 
l’erreur de la théorie classique. Celle-ci s’appuie en 


D 


aps 


ie ho 


Héfrege.re de 


2 Verde 


Not 


effet sur un autre postulat que l’on doit appliquer 
avec circonspection : 

« Toute partie d’un solide, quelle que soit la 
façon dont elle est limitée, peut être envisagée à son 
tour comme un corps auquel les théorèmes généraux 
de la mécanique sont applicables.» Ce postulat 
paraît évident dès que l’on imagine le solide divisé 
en petits parallélipipèdes élémentaires séparés les 
uns des autres. Mais cette division est délicate ; 
d’une part elle nous conduit à représenter les 
actions entre les différentes parties macroscopiques 
d’un corps par un système de forces centrales car les 
forces centrales sont les seules que l’on puisse se 
représenter simplement ; d'autre part, elle incite à 
certaines confusions entre forces et tensions. 

Considérons un cristal comme le chlorure de 
sodium dont tous les atomes sont situés en des 


| centres de symétrie. Des tensions homogènes, 
_ définies comme étant les forces agissantes de part 


et d’autre d’une surface unitaire, induisent une 
déformation homogène du réseau cristallin. Les 
atomes restent situés en des centres de symétrie, 
ils ne sont donc soumis à aucune force et restent en 
équilibre, mais cela n’exige nullement qu’il y ait 
une relation entre les composantes du tenseur des 


. efforts. Notre raisonnement ne porte que sur les 


atomes situés assez loin de la surface, dans une 
région où le réseau est bien triplement pério- 
dique, c’est-à-dire au delà de quelques dizaines 
d’angstrôms de la surface. Dans cette couche super- 
ficielle les atomes ne sont en équilibre que grâce à 
l’application des forces extérieures qui créent la 
déformation. La seule condition d'équilibre qu'il y 
a lieu d’écrire est une relation entre les forces 
extérieures appliquées qui doivent former un sys- 
tème réductible à zéro ; ce n’est pas une relation 
entre composantes du tenseur des efforts. Il importe 
de ne pas confondre la force à laquelle est soumise 
un atome avec la tension locale. 

Le fait se généralise aisément si chaque atome 
n’est pas situé en un centre de symétrie car c’est la 
stabilité même du milieu cristallin qui exige que 
l'énergie interne soit stationnaire par rapport à 
toute déformation homogène. Ainsi la condition 
d'équilibre statique est simplement la suivante : 

Le système des forces extérieures appliquées 
au solide doit être identiquement nul. Chaque 
atome doit être soumis à une force nulle dans le 
solide déformé. 

Il en est ainsi dans les déformations représentées 
figure 6 où un cristal orthorhombique taillé paral- 
lèlement à ses axes binaires est soumis tantôt à une 


À 2 
déformation hétérogène de la forme f,, = tsin _ 
Le ; À ; 
tantôt à {, — t sin es Les énergies des défor- 


mations sont respectivement 


[Atrane fe do et IE 1 dv. 
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La nouvelle théorie impose, en accord avec l’ani- 
sotropie du cristal mais en désaccord avec la théorie 
de Voigt, une énergie de déformation différente 
dans les deux cas. 

Lorsqu'on isole un parallélépipède élémentaire 
de volume dv dans un solide soumis à des' défor- 


Fe 


o S, 


F1G. 6. — Parallélépipède soumis à deux déformations qui, 
d’après la théorie de Voigt, nécessitent la même quantité 
d'énergie élastique. Nous estimons, au contraire, que ces 
énergies sont différentes. 


mations hétérogènes, il faut calculer la résul- 
tante F, dv et le couple F, dv auxquels il est soumis 
en sommant les forces exercées sur chaque atome 
et les moments de ces forces. Le calcul de la résul- 
tante F, dv ne présente pas de difficulté : 


_ Z Tes 
Fa dv = 2 ed 
Pour calculer le couple, plaçons-nous dans le cas 
particulier où seul T,, n’est pas nul et décomposons 
le volume dv en volumes infiniment plus petits dv, 
la composante suivant OË, est : 


dE 
72 DT 
D de dé, a f à, Fe 4 Da 
2 
soit 
2 ak 2 d°Ti 
T', dv FD de? Pésae 


Ainsi le couple s'exprime en fonction des dérivées 
secondes des tensions et est négligeable contrai- 
rement à la formule classique 


Tdi 7, dp, 


V. Possibilité d’un tenseur d'efforts dissy- 
métriques en élasticité dynamique. — On observe 
des déformations comme celles de la figure 6 lors- 
qu’une onde plane se propage dans un milieu 
cristallin indéfini : la condition à respecter est que 
chaque atome satisfasse au principe de l’inertie. 
Du point de vue macroscopique chaque volume 
élémentaire doit y satisfaire à son tour mais les 
couples qu'il y aurait lieu de considérer s'ils 
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n'étaient négligeables sont ceux dont nous avons 
donné l'expression en fonction de la dérivée seconde 


Ad Do 
des déformations. Si l’élongation est u , la densité p 
et le temps £, il vient : 


Cherchons les solutions correspondant à la pro- 
pagation d’une onde plane : 


; 2 
Ux = Ux exp (ot T4 Es) 


w est la pulsation, À la longueur d’onde acoustique, 
SE 


= . . 

U le vecteur élongation de module maximum et q 
un vecteur unitaire normal. aux plans d'ondes. Il 
vient : 


É LT? 
po? Us = 5e 2 Toy.es 47 48 Us 


et si V est la vitesse de l’onde élastique : 


PVUs = ZE Moy.ps 4v 98 Un. 
YBR 


C’est l’équation fondamentale de la propagation 
d’une onde plane. Remarquons que les coef- 
ficients Moyes et user sont multipliés par le 
même facteur qg, gs et n’interviennent donc que par 
leur somme 


2Moy,es = Moyes + os,py- 


Ce sont les coefficients 9,85, symétriques en «y 
et 83, en y et à, donc en « et B, au nombre de 36 
pour un cristal triclinique, qui déterminent la 
vitesse de propagation d’une onde plane. Si l’onde 
n’est pas plane, ‘il faut faire appel aux 45 coe- 
fficients élastique, comme dans les problèmes 
statiques. 


VI. Conséquences de la dissymétrie du tenseur 
des efforts. — Nous avons publié par ailleurs [8] le 
tableau des nouvelles constantes élastiques et piézo 
électriques pour les 32 classes de symétrie du 
milieu cristallin. Notre tableau des constantes élas- 
tiques a été repris sous d’autres notations par 
Sir C. V. Raman et Viswanathan [7]. L’examen de 
ces tableaux conduit aux conclusions suivantes qui 
sont toutes en accord avec l’anisotropie cristalline 
et le principe de symétrie : 


19 La partie centrale d’un cristal taillé en forme 
de cylindre peut tourner légèrement lorsqu'on 
appuie sur ses bases ; c’est le cas de la wulfénite 
comprimée suivant son axe quaternaire ou du 


Note 


quartz, de la tourmaline et de la calcite comprimés 
normalement à leur axe ternaire. 


20 Il est possible de créer une polarisation piézo-, 


4 


électrique parallèle à certaines directions cristal-m 


lographiques auxquelles la théorie de Voigt ne 
reconnaît pas cette possibilité : c’est le cas de l’axe” 
quaternaire du sulfate de nickel et de l’axe ternaire, 
du quartz « ou 6. De plus, les cristaux de la 


classe 29 (0,432) appartement à l’hémiédrie énan-" 


tiomorphe cubique qui, quoique dénués de centre MA 
de symétrie, n’étaient pas considérés comme piézo= 

électriques, peuvent être polarisés par des efforts M 
dissymétriques. 


30 Chaque fois que suivant deux directions L 
cristallographiques perpendiculaires de natures 
différentes, OË, et OË, peuvent se propager deux 


6 E | 
o tre F 


F1G. 7. — Propagation de deux vibrations transversales M 
pures suivant deux directions différentes. S'il s’agit 
d’ondes se propageant dans un cristal de phosphate 
monoammonique orienté de telle sorte que OË, soit un 
axe binaire et OË, l’axe quaternaire inverse, nous avons 
montré que, contrairement à la théorie classique, les 
vitesses de propagation ne sont pas égales. 


vibrations purement transversales dirigées, la pre- 

mière suivant O6, la seconde suivant OË, (fig. 7),! 
la théorie classique leur affecte la même vitesse de 
propagation. ‘4 

Nous estimons, au contraire, que ces vitesses sont. 
différentes ; l’expérience a été réalisée sur le phos- « 
phate monoammonique [9]. Les vitesses de pro- 
pagation mesurées par la méthode de Bergmann- 
Schaefer sont les suivantes : 2,06 km/s suivant 
l’axe binaire et 2,40 km /s suivant l’axe quaternaire! 
Dans cette expérience simple il est évident que les 
perturbations apportées à l’édifice cristallin sont 
différentes pour les deux types d’ondes et il est 
logique que’leurs vitesses de propagation soient 
différentes. 

Il est probable que l'incertitude qui règne actuel- 
lement sur les valeurs des constantes élastiques … 
cristallines provient, en grande part, d’une théorie. 
insuffisante. Il importe que les résultats expéri- 
mentaux, soient donnés dorénavant dans le cadre 
de la nouvelle théorie, laquelle est plus simple que 
l’ancienne car elle est plus proche de la symétrie 
des phénomènes. 


Manuscrit reçu le 15 juin 1956, 
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REVUE DES LIVRES 


Gamow (G.), Monsieur Tompkins explore l’atome. (1 vol. 
16 X 22 cm, 120 p., 24 fig., Dunod, Paris, 1954, 440 F.) 
M. Tompkins, dont nous avions fait la connaissance au 

pays des merveilles relativistes et quantiques, rencontre ici 

le démon de Maxwell, d’illustre mémoire, puis le Révérend 

Père Paulini, chargé de surveiller la circulation des électrons 

et enfin le vieil artisan sculpteur, qui tourne, avec du bois 


||} quantique, des protons et des neutrons et les assemble sans 


colle pour en fabriquer toutes sortes de nuclides. 
Comme dans les autres ouvrages de cette série, les diva- 


| gations de M. Tompkins ont pour origine les conférences 


de vulgarisation de son beau-père, le professeur à barbe 


Ë blanche. Ici le cycle de conférences porte sur l’histoire dela 


L des merveilles a conquis l’immortalité.… 


. physique atomique et de la chimie nucléaire. Il est traité en 


quatre points : « Réalité de l’atome » (expériences de Jean 
Perrin et de Bragg) ; « L'intérieur de l’atome » (théorie de 
Niels-Bohr) ; « Des trous dans le vide » (théorie de Dirac) ; 


. « Le monde intime du noyau » (réactions nucléaires). 


On peut penser que les travaux scientifiques de G.Gamow 
pèseront plus lourd, aux yeux de la Renommée, que les 
aventures de M. Tompkins. Mais en cette matière, la pro- 
phétie est incertaine, comme on peut le voir par l’exemple 
de G. L. Dodgson, dont les traités de mathématiques parais- 
sent aujourd’hui bien oubliés, alors que son Alice au pays 


M. JEssEL. 


La chimie des hautes températures. (Colloque National du 
CNRS: Paris, mai 1954.) 1--vol., 16 X 24 cm, 
264 pages, C. N. R. $., Paris, 1955, 2 000 F. 

Ce livre réunit les communications effectués au Colloque 
National de mai 195% sur la Chimie des Hautes Tempé- 
ratures. : 

Les techniques modernes de production et de mesure des 


. hautes températures ainsi que les matériaux réfractaires 
- nouveaux constituent le sujet d’une partie des commur- 


nications ; l’autre étant (consacrée à l’étude de phénomènes 
physiques et électrolytiques à haute température. 
R. Buver. 


HozLAzNper (Alexandre), La biologie des radiations. 
(Radiation Biology), volume III. Lumière visible et 
presque visible. (1 vol. relié de 765 p., 15 X 283 cm, 
McGraw-Hill, New-York, 1956, prix s 10 (en anglais). 
L'action biologique de la lumière constitue un sujet 

passionnant non seulement pour le biologiste, mais aussi 

pour le physicien. La table des chapitres de ce volume, 
traités chacun par des auteurs différents, est rendue 
attrayante par l’énoncé des problèmes traités et par les 
personnalités qui ont été appelées à les rédiger. Malheu- 
reusement, le lecteur physicien (et non biologiste) éprouve 
une certaine déception ! Il se trouve en présence d’un long 
catalogue de résumés de mémoires expérimentaux, et 
éprouve une grande peine à se faire une idée exacte des 
faits acquis avec certitude. Bien entendu, cela n’est pas 
général ; la définition de la photosynthèse, phénomène 
fondamental de la biologie végétale, est claire ; et le phy- 
sicien est surpris par la valeur élevée des rendements d’un 
phénomène mettant en jeu plusieurs quanta de lumière par 
molécules. Table des matières : 1. Échange d’énergie dans 
les photoréactions. 2. Structure électronique et excitation 
des polyènes et porphyrènes. 3. Génération, contrôle et 
mesure de l’énergie de radiation visible et presque visible. 

4, Rendement énergétique en photosynthèse. 5. Mécanisme 

de la photosynthèse. 6. Absorption, action et spectre de 

fluorescence des pigments de photosynthèse dans les cellules 
vivantes et en solution. 7. Formation et accumulation de la 
chlorophylle. dans les plantes. 8. Réduction des nitrates. 

9. Phototroprisme. 10. Photopériodisme. 11. Germination 

des semances. 12. Viscosité, perméabilité et écoulement 

photoplasmique. 13. Phénomènes électriques dans la vision. 

14. Photorécepteurs invertébrés. 15. Action photodynamique 

et ses effets pathologiques. 

J. Winter. 
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VALEUR LIMITE, POUR LES TRÈS FAIBLES ANGLES, 
DE L'INTENSITÉ OBSERVABLE DU RAYONNEMENT X DIFFUSÉ PAR LES FLUIDES 


Par G. FOURNET, 
École Supérieure de Physique et Chimie. 


Sommaire. — La valeur limite de l’intensité du rayonnement X diffusé par un fluide sous de 
très faibles angles a été récemment contreversée. Après avoir donné un critère permettant d’extraire 
« l'intensité observable » de l’intensité totale, nous montrons que les formules classiques sont en 
accord avec nos résultats (en ce qui concerne l’intensité observable) obtenus au moyen d’une 


analyse plus fine des phénomènes. 


1. Introduction. — L'expression de l'intensité 
du rayonnement X diffusé par un ensemble de 
particules identiques, à symétrie sphérique, est 
généralement mise sous la forme 
Lots) = Le(s) N F°(s) 

sin 2Trs 


1 + m [PP m1 


kr? dr} (1) 


où s désigne le rapport 2sin0 /À (20 etX sont respec- 
tivement l’angle de diffusion et la longueur d’onde 
du rayonnement), 

I. (s), l'intensité diffusée par un électron dans les 
conditions considérées, 


N,le nombre moyen de particules en examen 
situées dans le volume irradié, 

F(s), le facteur de structure des particules. 

. Les notations n, et P(r, n.,) sont classiques ; les 
définitions correspondantes sont rappelées dans le 
corps du présent article. 

Nous avons déjà indiqué [3], en perfectionnant 
les idées de quelques auteurs, qu’il faut ajouter 
un terme à l’expression (1) pour obtenir l'intensité 
totale diffusée Z:(s). Ce terme complémentaire ne 
présente de valeurs notables que pour les très 
faibles angles, largement inférieurs au degré ; il 
est pratiquement négligeable pour les angles obser- 
vables et c’est pour cela que nous avons désigné [3] 
l'expression (1) par «intensité observable » Z,(s). 

La notion de valeur limite Z,(0) de l’intensité 
observable provient des difficultés expérimentales 
et théoriques de mesurer correctement les phéno- 
mènes de diffusion aux très petits angles. Cette 
valeur limite est assez bien définie expérimen- 
talement par le prolongement « au sentiment » 
d’une courbe expérimentale. 

La relation (1) permet d'obtenir pour cette 
valeur limite 


150) = (0) N F2(0) [1 + n; st [P(r, n;) —1] 4rr2dr$- 
(2) 


Cette éxpression n’est pas uniformément admise. 
Pour les auteurs qui en sont partisans la valeur 
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numérique de (2) est encore sujet de discussion ; 
certains indiquent pour cette valeur | 


I5(0) = 1e(0) N F2(0) n, KT y (3) : 


où y est le coefficient de compressibilité isotherme, 
tandis que d’autres fournissent zéro [1, ?, 5]. 

En travaillant sur ce sujet nous nous sommes 
aperçus que cette question est très délicate et que 
la plupart des calculs qui ont, été effectués dans ce 
domaine ne sont pas assez rigoureux pour'entraîner 
une conviction universelle. Nous nous proposors 
dans cet article de montrer que la valeur limite est 
définie par (2) et (3) pour le modèle de fluide ré | 
nous définissons dans le paragraphe 3. 


2. Définition de l’intensité observable. — Avart 
de définir ce modèle nous désirons revenir sur la 
notion d'intensité observable pour en donner une“ 
définition précise. Les expériences de diffusion des 
rayons X et leur dépouillement sont toujours 
effectués en supposant que, les corrections d’absorp- 
tion étant faites, les intensités observées sont 
proportionnelles aux nombres de particules en 
examen. Nous définirons donc l’intensité obser- 
vable comme étant la portion de l'intensité totale 


qui varie linéairement avec le nombre moyen W de 
particules irradiées. Avec cette définition 
l'influence de la forme du volume Y (la notion de 
forme excluant lanotion d’étendue) ne peut se faire 
sentir que dans le terme complémentaire 2, (s). 


3. Définition du modèle de fluide. — Nous consi- 
dérons un ensemble de AV, particules contenues 
dans une enceinte de volume V,. Cet ensemble est 
défini de la façon suivante : les particules sont 
toutes identiques ; la statistique classique de 


- Boltzmann est utilisée ; l’énergie totale du fluide 


peut s’exprimer sous la forme 


Be Sims Sut! ÿ ÿ du: (} 


k=1 K+j=1 


la première somme concerne l’énergie cinétique 
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| (m est la masse d’une particule) ; elle comprend 


| VPaction des parois limitant le volume V,;. 


N, termes. La seconde somme comprend également 
N, termes, Uz étant le potentiel d’où dérive le 
champ de force extérieur agissant sur la kiè®e par- 
ticule ; nous négligerons l’action de la pesanteur 


! de sorte que Ü ne sert qu’à symboliser sous forme 


d’un «champ répulsif très localisé mais intense» [7] 
Le 
troisième terme de (4) comprend effectivement 


| D | valeurs du potentiel ®,; chaque valeur 
correspondant à chacun des  — couples 


formés par les N, particules. Les particules pos- 


|: sédant la symétrie sphérique, nous supposerons 


que le potentiel ®,; ne dépend que de la distance r,, 
entre les particules Æ et 7. Les valeurs de la fonc- 
tion ®(r) sont extrêmement élevées pour les petites 
valeurs de r (traduction de l’impénétrabilité des 
particules) et très faibles et même sensiblement 
nulles pour les valeurs de r grandes devant les 
dimensions des particules en examen mais extrême- 


ment petites devant les dimensions de l’enceinte. 


Pour décrire les propriétés du fluide nous n’utili- 


—_— 
serons que deux fonctions : la fonction n. (R,) 


telle que 
n:(R) dÆy 


est la probabilité de trouver un centre de particule 
Les 
dans l’élément de volume d A,entourant l'extrémité 


— î > — 
du vecteur À, ; la fonction ni À) telle que 


À, 


est la probabilité de trouver à la fois le centre d’une 


À) ak ak 


— à 
particule dans l’élément d À, et le centre d’une autre 


_— 
particule dans l’élément dA,. 
Nous allons maintenant donner les résultats que 


_ Yvon [7] a obtenu sans aucune approximation, 


: — 

pour ce modèle. La fonction n(R.) possède une 
valeur constante 7, dans à peu près toute l’étendue 
de l'enceinte V, et ne présente des valeurs diffé- 
rentes de cette constante que sur les bords de 
l’enceinte V, ; dans notre modèle ces variations 


sont dues au potentiel U. 
. TE > . 
La fonction nA\R;, R.) peut s’exprimer sous la 


forme 
> > 
né, R; 
> > > > 
= (A) n, (À) [PG, À) +R o(, RÀ;) + | (5) 
0 0 É 
Dans un domaine V, un peu plus restreint que 


—> 
celui où n(R,) est constant et égal à », la fonc- 
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tion PR, R R) ne dépend. que de »., et de la dis- 

tance r;; entre les centres des particules en examen ; 

nous noterons cette fonction par P(r,;, n;) ; sur les 
+ 


—> 
bords du volume V,1la fonction p(r,, À) présente 
quelques variations par rapport à celles qui seraient 
déterminées à partir de P(r,;, n;). 


La fonction QRR, R;)n’a pas été complètement 
déterminée par Yvon ; sil a été simplement montré 
que pour les distances : r; grandes devant les dimen- 
sions des particules en examen on obtenait 


QU, À) = —_ (8) o(R). (6) 


> 
La fonction (R,), dont la valeur est constante 
dans le domaine V, (cette valeur ne dépend que 
de ñ.), présente des variations sur les bords de V,. 
Pour résumer les résultats d'Yvon que nous 
venons d'exposer nous pouvons indiquer que dans 


ES 
une très grande partie V; de V, la fonction n (#) 
présente une valeur constante n, très voisine 
de W,/V,. Dans ce domaine V,; la fonction 


— — 
—. À) ne dépend que de la distance r;; entre les 
centres des deux particules # et 7 considérées ; 
pour les valeurs de r;; comprises entre 0 et plusieurs 
diamètres de particules (20 par exemple) nous 
avons 


= nf Pre m3) + 5 Qu m3) + ..]5 O0) 


pour les grandes valeurs de r;; nous avons plus 
simplement 


naG, À) = nana) = ni — 5 6%) +]. (6) 


Le terme en 1 /N, des développements (7) et (8) 
n’a jamais été considéré dans des calculs de dif- 
fusion ; nous montrerons qu’il peut être essentiel 
aussi donnons nous la valeur de b?(n,) déterminée 
par Yvon 


bü(ns) = 1 + m [© Pin, nj-4jérrédr +.) (ds) 
0 , 
le premier terme négligé étant en 1 /N,. 
4. Calcul de l’intensité du rayonnement diffusé. 
— L'expression générale et complète de l'intensité 
moyenne diffusée par un ensemble de NW, particules 


contenues dans un volume V, dans le cas où le 
volume irradié est V < Vi est “fournie par [4] 


AN(S) = Le) F4(s) fo) ar 


> > > 
+ ff, cos [re Re R)ln(Re, À) deal (10) 
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ES 
où le vecteur s , de module s, est dirigé suivant la 
bissectrice intérieure des directions des rayon- 
nements incident et diffusé. 
Nous allons considérer le cas où le volume irradié 
est compris dans le volume V, (cf. paragraphe 3) ; 
nous obtenons à partir de (10) 


is) = Lis) Fè() far 
E RE ee Fes (HR) n | Pr n) 
1 mare 
+7 A na) + “| aa | (11) 


Nous pouvons maintenant considérer le pouvoir 
diffusant 


tt) 
Ta{s) F?(s) 


comme constitué d’une somme de cinq termes que 
nous désignerons par À, B,C,DetEÆE 


4 J dr, 
La Pre % 


SR US TS 

B= nt ff Pins n1) —11coslax $ (RER) dx dk 
SES A > 

Con JL 1-c0slex s (R—K)]lak à; 


Æ (QU; 1) + b#(n)] 


AE D) li x 
cos [ar S (28) d y dR; 

i > > > > > 
+ LE bn)] cosl2x s (Hi —#)]aH; dr; 


Le calcul du terme À est très simple et fournit 
la valeur r,V ; cette valeur représente le nombre 
moyen { de particules situées dans le volume . 

Les termes C et E sont proportionnels à l’inté- 
grale double 


PE ARE 
[fees Lors (fK)la ar, 


qui peut s’interpréter comme représentant l’inten- 
sité diffusée (cf. [3] par exemple) par une particule 
homogène, de densité électronique égale à l’unité 


(12) 


A 
et de volume V. Nous désignons par V? 1 la 
valeur de cette intégrale double ; d’après cette 
définition la valeur de &(0) est l’unité, nous savons 
d'autre part que pour une « particule » d’un 
volume V dont l’ordre de grandeur est de plu- 


ee 
sieurs mmÿ la fonction i( ) décroît extrêmement 
rapidement avec s. Avec ces notations les termes C 
. et E peuvent maintenant s’écrire sous la forme 


Cns Pe SNA) 
2 72 
= LE bs{nlh 2iCe) a im) ils). 
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Le calcul du terme B s'effectue par l’inter- 


médiaire des stades suivants : 


— nous prenons comme nouvelles variables 


d'intégration 


3 Fe _ > 
Tkÿ —= Ry— h; et K; 


> > 1 —> 

Fi nà [ar J n;) — 1] cosl2r s .rgl'aré. 
a LS 
volume permis à r#; 


— la fonction P(r;,;, n;,) — 1 étant pratiquement, 


nulle pour des valeurs de r;; très faibles devant les 
dimensions de V nous pouvons remplacr le 
—# 


« volume permis au vecteur 74; » par un volume 
infiniment étendu. 
L'intégrale par rapport à r,; devient ainsi 


; > 7 —> 
de [P(rr5, n1) — 1] cos[2x s 6] dry; 
CO 


où nous pouvons maintenant supprimer l'indice kj 
et considérer simplement l'élément  d’inté- 


— 
gration dr ; une intégration effectuée sur des 
_ 


vecteurs r de même module mais de toutes les 
orientations possibles fournit pour (13) 


co sin 2Tsr ; 
[ [P(r, A rm Ne dr 


ce qui permet tous calculs faits, d'écrire B sous la 
forme 


sin 2Trrs 


B N {1 1 
= ni [PE n) IS 


4knr? dr. 


Par des raisonnements analogues la valeur du 
terme D est 


D = M [ULQU m) + br] =. 


krr?dr 


puisque la fonction Q(r, n;) + b?(n;) est prati- É 


quement nulle pour des valeurs de r petites devant 
les dimensions du volume Y. 


5. L’intensité observable. — Grâce au critère 
que nous avons précédemment défini nous voyons 
que l’intensité observable, groupant les termes À, 
B et D, est 


Ils) = N Le(s) F?(s) 
ES) sin 27rs : 
É + ra Î [Ptr, LE en ps dr 
ni sin 27Trs 


+R Qi m3) + bUm)) amrtdr+ | (14) 


27 rs 


— 
tandis que le terme complémentaire A 5) s'écrit 
sous la forme 
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(13) 1 


smittintié à mm 


| 
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1,05) = Lt) F5) N2i(5) [1 — 


b?(n;) Le 
0° 


ne (15) 


L'intégrale en Q(r, n,) étant du même ordre de 


* grandeur que l'intégrale en P(r, n;) nous voyons 


que l’expression de l’intensité observable est bien 
celle que nous avons donnée (1) au début de cet 
article, à des corrections tout à fait négligeables 


| en 1/N, près. 


Les termes en Q(r, n,) (et sa limite b?(n.,)) jouent 
un rôle quand on considère un angle de diffusion 
extrêmement faible et tendant vers zéro. Dans ce 
cas le deuxième terme de Z, (0) 


__ 24 EN .NI4(0) F2(0) b2(n;) 


— 1410) F*0) N?i(0) T, 


peut être du même ordre de grandeur que la valeur 


_ limite Z,(0) puisque b?(n.) est de l’ordre de l’unité 


l’unité (cf. (9)). Si nous nous limitons dans l’expres- 


. sion de l'intensité totale 1,(0) aux termes dont les 


ordres de grandeur sont ÆV° et N nous pouvons 
donc négliger le terme en Q(r, n,) et b?(n,) pro- 
venant de l'intensité observable mais nous devons 


| prendre en compte le terme en b?(n,) du terme 
complémentaire Z,(0) pour écrire 


I4(0) 


FA) F0 | +m f [P(r, n;) —1]4nr? dr + ... 


SN 
HE PEN. bn) + . (16) 


Cette expression va nous permettre de montrer 
comment certains auteurs [1], [2], [5] ont pu croire 
que la limite Z,(0) (qui correspond à la première 
ligne du second membre de (16)) était nulle. Consi- 


_dérons le cas où le volume V, contenant les parti- 


cules est totalement irradié ; l'amplitude diffusée 


à l'angle nul est V,V/1.(0) F (0) et nous avons donc 
alors : 
I:(0) — 


N$ Le(0) F0). (17) 


Les calculs que nous avons développés dans le 
paragraphe 4 ne sont valables en toute rigueur que 
lorsque le volume irradié V est compris dans le 
volume V;, lui-même plus petit que V,. La façon 
dont nous avons défini V, montre que le 


. rapport V,/V, est d'autant plus près de l’unité 


que V, est plus grand ; nous pouvons donc dans 
le cas où V = V, employer nos résultats, et en 
particulier l’expression (16), en ne commettant sur 
les termes écrits de cette expression (16) que des 
erreurs aussi faibles que nous le désirons. Nous 


obtenons ainsi avec V — N, (cf. (17)). 
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I4(0) Ca Es Le PS 2 
EL FQ = Vi + [ LP(r, n1)— 1] &rrt dr +. 
NE Ne bn) ce: 


Pour satisfaire cette égalité quand le terme 
en b?(n.,) n’a pas été pris en compte, il ne reste plus 
qu’à écrire que 

1 + cs (PinaleeAllanrs dr 


est égal à zéro ; c’est ainsi que cette relation que 
nous condamnons a été obtenue. Si par contre on 
considère le terme en b?(n.,) et sa valeur (9) on 
obtient l'identité N, = N,. 

On trouvera dans de nombreux ariicles le passage 
de (2) à (3) (cf. par exemple [6], [7]et [31). 


6. Conclusions. — Nous pouvons maintenant 
voir pourquoi le calcul de la limite de l’intensité 
observable pour les très faibles angles est particu- 
lèrement délicat. Si l'intensité observable Z,(0) 
est de l’ordre de W, le terme complémentaire Z,(0) 
et l’intensité totale Z:(0) sont de l’ordre de A? ; 
cf. les relations (14), (15) et l’égalité 


In(0) = 100) + 1)(0). 


Une erreur relative de 1 /N sur le terme complé- 
mentaire Z,(0) soit une erreur absolue de l’ordre 


de NV? x += N est donc du même ordre de 


grandeur que l'intensité observable ; cette 
remarque permet de nous rendre compte de la 
précision et de la rigueur nécessaires dans le calcul 
du terme complémentaire et nous explique le rôle 
essentiel, contrairement à ce que l’on pourrait 
penser, du terme en Q(r, n.)/N, du développement 
de r,(r) dans ce calcul ; Yvon lui-même reconnait 
que la « différence entre 


na(r) = n£ P(r, m3) + ni 


Q(r, n:) | 
Se ARE t 


et nr? P(r, n.) est un peu subtile ». 
L'expression (16) nous montre que le rôle 
de b?(n;,) peut être négligé quand NW, est grand 


devant N ; c’est pour cela qu’il est relativement 
facile d'obtenir [3] des conclusions exactes dans le 
cas où le volume irradié V est petit devant le 
volume V,;. 

Nous pensons montrer dans un prochain article 
comment il est possible d'appliquer un certain 
nombre des idées que nous venons de développer 
à la diffusion des rayons X par les substances 


amorphes. 
Manuscrit reçu le 17 avril 1956. 
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LE RENDEMENT DFE FLUORESCENCE DE LA COUCHE x. 
MESURES SPECTROMÉTRIQUES SUR %Tc* (6,04 h) ET ‘5ln* (4,5 h) 


Par JEeannE LABERRIGUE-FROLOW et Pierre RADVANYI, 
Laboratoire de Physique et Chimie Nucléaires, Collège de France. 


Sommaire. — Nous avons étudié avec un spectromètre 8 les électrons Auger Æ et les électrons 


de conversion! Æ dans les transitions isomériques de °3Tc* et 1} 


$In* afin de déterminer le rendement 


de fluorescence de la couche K, Rx, de ces deux éléments. LLé valeurs obtenues sont : pour Te, 
Rx = 0,70 + 0,03 et pour In, Br — 0,87 + 0,03 ; les intensités relatives des raies Auger Æ-LL, 
K-LX et K-XY sont respectivement : 1: 0,48 : 0, 05 pour Te, et 1 : 0,42 : 0,05 pour In. En nous 
basant sur les déterminations de RK publiées depuis 1946 nous avons Calculé une nouvelle courbe 
semi-empirique, par la méthode des moindres carrés, en utilisant la formule théorique de Burhop 


( 
1 RK 


I. Introduction. — Il est nécessaire dans divers 
problèmes de physique nucléaire de connaître avec 
exactitude plusieurs constantes atomiques, tel le 
rendement de fluorescence de la couche X, Rz ; 
sa valeur intervient notamment d’une manière très 
sensible dans certaines déterminations du rapport 
capture L 


capture X 
tivistes et relativistes de RZ (et des intensités rela- 
tives des raies Auger XÆ) ont été faits, en 1935 et 
1936, avec des fonctions d'onde hydrogénoiïdes en 
appliquant les règles de Slater pour tenir compte 
de l’efflet d'écran [2-3-4-5]; mais la précision 
obtenue avec ces fonctions d’onde paraît insuf- 
fisante pour les besoins actuels des expéri- 
mentateurs. Ceux-ci tirent donc en général À; de 
différentes courbes semi-empiriques ou empiriques. 
Cependant les résultats expérimentaux sur Rx 
présentent relativement une grande dispersion. 
Aussi nous a-t-il paru intéressant, de déter- 
miner quelques rendements de fluorescence de 
la couche X. 
La méthode que nous avons employée consiste 
à étudier à l’aide d’un spectromètre £, les électrons 
Auger À émis à la suite de la conversion interne 
dans la couche X d’un rayonnement y nucléaire. 
Le rapport de l'intensité des raies Auger XÆ à celle 
de la raie de conversion X donne la valeur du 
rendement Auger K:R,x, d’où on déduit le 
rendement de fluorescence X : Rx = 1 —R,x. 
Cette méthode a pour avantage de procéder par 
comparaison directe ; elle permet ainsi des déter- 
minations précises. Elle a déjà été utilisée par 
différents auteurs [6-7-8] surtout pour des éléments 
lourds, mais dans ce domaine la valeur de Rx 
varie peu avec Z ; nous l’avons appliquée à des 
éléments de Z inférieur à 50, domaine de grande 
variation de À4 en fonction de Z. 


[1]. Des calculs théoriques non rela- 


IT. Choix et préparation des sources. — Nous 
avons choisi pour notre étude deux isomères : “Tec 
de période 6,04 heures et *5ln de période 


1/4 
| — À + BZ + CZ; les valeurs trouvées pour les coefficients sont : 
B — 0,0332 


C—=—1,14 X 10. 


4,5 heures. Le rendement de fluorescence X du 
technétium n’avait pas encore été mesuré ; celui de 
l’indium avait été déterminé sur le radio: 
isotope 113 In par Harrison, Crawford et Hopkins [9 
au compteur proportionnel. 


STc* (&ouh) Rain* G,5h) 
À + (2 834 keV 
(2) Sie 0,142 MeV Me 
CR+) 
5% 
92+) 2 
0 de DIN 
0,63 
MeV 
VA 
38 Ru (stable) 5 Sn (stable) 


Fic. 1. — Schémas de désintégration 
de 2$Tc* et 1151n* d’après [10]. 


Il nous a semblé préférable de ne pas choisir 
pour notre expérience un nuclide se transformant 
par capture : il aurait fallu alors tenir compte de la 
capture L et des embranchements éventuels vers 
les différents niveaux et ceci aurait introduit une 
incertitude supplémentaire. 

%Te et 15[n dans leurs états fondamentaux sont 
radioactifs 5, mais avec des périodes très longues, 
et cette activité ne peut donc jouer aucun rôle dans 
notre travail. L’isomère *!5]n présente un faible 
embranchement $- d'environ 5 % donnant un 
fond continu peu intense. 

Le technétium et l’indium ont été préparés et 
séparés au Commissariat à l'Énergie Atomique à 
partir de molybdène et de cadmium irradiés en 
neutrons, Te a été séparé parfaitement sans 
entraineur ; dans l’extraction de In un peu d’en- 
traîneur restait présent. 

La difficulté principale de cette méthode par 
spectrométrie $B réside dans la détermination 
correcte de l'intensité des électrons Auger. Cette 
difficulté est d'autant plus sensible que l’élément 
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étudié a un Z plus petit. Il faut veiller particu- 
lièrement à réduire l’épaisseur de la source et de 
son support et à bien connaître l’effet de la fenêtre 
du compteur servant de détecteur. 

Nous avons préparé nos sources par sublimation 
à l’air sur des feuilles de LC 600 d’épaisseur envi- 
ron 20 ug/cm?. Les sources, d’un diamètre bien 
défini de 2 mm, n’avaient pas d'épaisseur mesurable 
Le dépôt de Tc ainsi obtenu était absolument 
invisible à l’œil nu ; le dépôt de In était visible et 
devait donc avoir une épaisseur plus grande ; cela 
est confirmé par la forme des raies Auger K 
du côté des faibles énergies. 


III. Dispositif expérimental. Facteur de correc- 
tion dû à la fenêtre du compteur. — Nous avons 
utilisé le spectromètre 8 type Slätis-Sieghahn de 
notre laboratoire avec une ouverture de dia- 
phragme correspondant à une résolution de 1,8 % 
pour nos sources. 

Il est nécessaire de connaître exactement le 
comportement du détecteur pour les électrons 
Auger examinés, c’est-à-dire dans un domaine 
d'énergie allant de 10 à 30 keV. Le détecteur utilisé 
était un compteur Geiger comportant sur sa paroi 


latérale une fenêtre de 8 mm de diamètre. Nous: 


avons préféré ne pas utiliser de grille de support 
pour la fenêtre ; celle-ci était alors constituée par 
une feuille de formvar d'épaisseur — 60 ug /em?. 
L'influence de la fenêtre a été déterminée en 


- relevant au début et à la fin des expériences le 
» spectre 8 de %S dont la forme (permise) est bien 


connue jusqu'à moins de 5 keV [11]. L'écart de la 
forme du spéctre 8 que nous obtenions par rapport 
à la forme réelle, qui est évaluée par le tracé de la 
droite de Fermi correspondante, donne le facteur 
de correction à apporter à nos mesures. Le fait 


œ 


Facteur de correction 


+ 


5 ——— - - RSS à D 


0 7 20 keV 30 


Fic. 2. — Correction due à la fenêtre du compteur 
(— 60 ug /cm?). 


d’encadrer ainsi notre expérience par deux étalon- 
nages nous assurait que la déformation de la 
fenêtre de formvar n’était pas importante au cours 
du temps. La source de %$, séparée parfaitement 
sans entraîneur par une méthode pratiquée et mise 


LE RENDEMENT DE FLUORESCENCE DE LA COUCHE X 


945 


au point par M. Marius Chemla [12], était déposée 
par sublimation sur une feuille de LG 600 
de + 20 ug /em? comme les sources de Te et de In. 

La figure 2 donne le facteur de correction obtenu. 
Les points relatifs aux deux séries de mesures se 
recouvrent bien, compte tenu des marges d’erreur ; 
l'erreur sur le facteur de correction dans le domaine 
d'énergie qui nous intéresse est de l’ordre de 3 %,. 
Le facteur de correction est pratiquement égal 
à Î au-dessus de 45 keV ; la fenêtre arrête tout 
rayonnement au-dessous de 5 keV environ. Nos 
résultats sont légèrement différents des valeurs 
indiquées par Saxon [13], ce qui indique une fois 
de plus la nécessité, pour des déterminations 
quantitatives précises, d’étalonner chaque appa- 
reillage utilisé. 


IV. Résultats expérimentaux [14]. — La figure 3 
donne les spectres obtenus pour les électrons 


W Raies Auger a) 89°.  Aaies de conversion 
KL [WA f 
[rs fusée keV Esoo te 


Fic. 3a. 
Raies Auger b) 15% 1n Raïes de conversion 
N/I W/ K 
306 keV 
L150 
H-LL 
Fig 20,1 keV 
+12 
‘e f-100 
H-LX 

Lg 23,2 keV 
Le 

+50 
L4 

K-XY 
26,5 keV 
2 : 
( l'(Ampères) 4 ; 1 (Ampères) 
CDR PRE LL EN OP OMR CE TTC ASE et EE 
Fic. 3b. 


F1G. 3. — Raïes Auger Æ et raies de conversion de **Tc et 
de 115*[n, mesurées au spectromètre B (corrigées de la 
décroissance des sources). 


6% 
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Auger K et les électrons de conversion de **Tc et 
de *l15]n compte tenu de la correction due à la 
décroissance des sources pendant la durée des 
mesures. On voit que les 3 groupes ÆÀ — LL, 
K—LX, K— XY de raies Auger X sont bien 
résolus, ainsi que les raies de conversion Æ et 
L + M du y de 140 keV dans le cas de Te et du y 
de 334 keV dans le cas de In (les électrons de 
conversion du y de 142 keV (1 %) de *°?Tce ne se 
distinguent pas de ceux du y de 140 keV). 


Fe 


N/1 x C 
L100 


K-LL 
15,4 keV 


a) Te 


K-LX 
17,7 keV 


H-XY 
20,1 keV 


1 (ampères)_, 
3 10 


Fic. 4a. 


K-LL 
b) tn 20,1 keV 


K-LX 
23,2 keV 


K-XY 
26,5 keV 


nn, n ñ VERS 1 à 1(ampeses) 
RICE DER TEE MENT PRE NU UE RD A AR La L 


Fic. 4b. 


FiG. 4. — Raïes Auger Æ de ‘’*Tc et de 115*[n corrigées 
de l’absorption dans la fenêtre du compteur. 
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La figure 4 reproduit les groupes Auger X de Te 
et de In, après avoir effectué la correction due à 
l'influence de la fenêtre du compteur. 

On constate sur le spectre des électrons Auger 
une queue vers les énergies basses due en grande 
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partie à l’absorption dans la source. Cette queue, 
peu importante dans le cas de °’Tc (elle repré- 
sente — 10 % de l’ensemble des Auger), lest 
davantage dans le cas de In par suite de l’épaisseur 
plus grande de la source. La surface de chaque 
spectre a été mesurée au planimètre. Dans le cas 
de In on a défalqué le fond continu dû au 
spectre 8. Le rapport de l’aire des raies Auger X 
à l’aire de la raie de conversion X donne pour le 
rendement Auger : 
pour 

22Tc : (R4x) = 0,30 + 0,03 
et pour s 

1151n : (Ra) = 0,13 + 0,03. 


donc pour les rendements de fluorescence : 


(Rx)z = 43 — 0,70 M 0,03 
et 


(Rr)z = 49 = 0,87 + 0,03. 


Pour ‘ln, Harrison et al. [9] avaient obtenu 
Rx = 0,82 + 0,02 

Si l’on compare nos valeurs à celles que l’on peut 
tirer par interpolation des mesures faites par 
d’autres auteurs sur des éléments voisins, on 
constate que notre valeur de RÀZ est pour Tec au- 
dessous de la moyenne des autres valeurs et pour In 
au-dessus. Ces deux écarts — en sens inverse — 
sont supérieurs à nos erreurs expérimentales. 

Dans le cas de Te, les électrons de conversion 
étudiés suivent l’émission d’un rayonnement y 
de 2,0 keV fortement converti. Il ne semble pas que 
cela puisse avoir un effet notable sur la valeur 
de RZ : le rayonnement de 2,0 keV ne peut, énergé- 
tiquement, être converti que dans la couche M et 
d'autre part la durée de vie du réarrangement 
atomique doit être beaucoup plus courte que celle 
du niveau nucléaire de 140 keV. 


INTENSITÉS PARTIELLES DES RAIES AUGER X. — 
Nos résultats expérimentaux permettent, faci. 
lement dans le cas de Te, un peu moins facilement 
dans le cas de In (voir fig. 4) de décomposer les 
groupes À — LL, K — LX et K — XY de raies 
Auger. Le tableau suivant donne les intensités 
relatives que nous avons obtenues ainsi que 
quelques valeurs théoriques et expérimentales dues 
à d’autres auteurs pour l’argent et l’indium : 


TABLEAU I 


INTENSITÉS 


RELATIVES DES Nos valeurs 


VALEURS EXPÉRIMENTALES 


VALEURS THÉORIQUES 
le G 


Autres auteurs 


GROUPES AUGER as LD ain [7] Ag [14], [15] Ag [8] In [2], [4] Ag [16] Ag 
Æ-LL 1 1 1 1 4 n 1 
K-LX 0,48 + 0,06 0,42 + 0,10 0,493 0,73 0,417 0,715 0,312 
K-XY 0,05 + 0,01 0,05 + 0,02 0,054 0,07 0,076 0,102 0,022 


à 
ù 


. 


Non 


On sait que d’après la théorie ces intensités 


- relatives ne devraient pas varier beaucoup avec Z 
” et on peut les comparer avec les résultats obtenus 
—…. pour des éléments voisins. Nos valeurs sont en 
- accord avec les mesures spectrométriques de [7] et 


[8] pour Ag-et In, mais par contre différent. des 
mesures de [15] et des résultats théoriques de 
Burhop et Pincherle [2], [4] et de Rubenstein et 
Snyder [16]. Nous pensons que le calcul théorique 
de ces intensités n’est pas encore satisfaisant et que 
lPaccumulation de données he Aentalee est 


+ intéressante pour ce problème. 


RAIES DE cONvERSION À ET L + M.— Les raies 


de conversion À et ZL + M que nous avons 


observées, entièrement séparées no falc ot 


…— aisément décomposables pour *9in; nous ont 


donné, comme valeurs des rapports Ana les 
résultats suivants : 


K 
(x) Oops d 0e 


de *®Tc en accord avec la valeur de 6,08 + 0,50 
donnée par [17]; ces mesures indiquent que le 


- rayonnement de 140 keV est du type Y.. Et : 


K 
ny) = 3178 L 015 pour.le y de 334 keV 


de *t51n en accord avec les valeurs de 3,85 données 
par [18] et 3,76 par [19], indiquant que le rayon- 
 rement de 334 keV de ln est du type 44. 


V. Discussion générale sur Rz7. — Différentes 


- courbes, donnant les valeurs de RL en fonction deZ, 


soit théoriques, soit semi empiriques ou empiriques 
ont été publiées. Il nous a semblé que même les 
courbes les plus récentes et les plus utilisées étaient 
discutables : en effet, Broyles et al. [8] ont tracé 
leur courbe empirique d’après l’ensemble des 
valeurs expérimentales depuis 1925 sans faire de 
tri. Burhop [20] pour sa dernière courbe semi- 


empirique a fait un certain choix et a utilisé un 
système de pondération qui ne nous ont pas semblé 
tout à fait justifiés. De plus, depuis la parution de 
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ces courbes plusieurs nouveaux travaux expéri- 
mentaux sur À4 ont été publiés dont un certain 
nombre pour des éléments lourds. 

Nous avons donc calculé une nouvelle courbe 
semi-empirique par la méthode des moindres 
carrés d’après la forme théorique proposée par 
Burhop [20] 

Re 
Ë ere LUE 


1/4 
] A LU R7Z C7 


où À correspond à une correction pour l'effet 
d'écran et C à l'effet relativiste. 

Nous n’avons retenu que les valeurs expéri- 
mentales de À4 obtenues depuis 1946 et effecti- 
vement publiées. Ce choix comporte sans doute une 
part d’arbitraire ; plusieurs mesures antérieures 
à 1946 sont fort probablement excellentes. Mais 
un choix et une pondération convenable sont 
difficiles, et, à condition de recouvrir un nombre 
suffisant de mesures, notre procédé doit être correct. 
Nous avons écarté les mesures indiquées comme 
incertaines par leurs auteurs. Sur la figure 5 sont 
portées les 42 valeurs expérimentales que nous 
avons utilisées pour le calcul des 3 coefficients de la 
formule théorique de Burhop. Nous avons donné le 
même poids à toutes les valeurs retenues, à l’excep- 
tion de celles de Laskar [21] auxquelles nous avons 
jugé nécessaire d’attribuer un poids moitié en 
raison de leur caractère semi-théoirique. 

La méthode des moindres carrés nous a donné 
pour les coefficients 


À = — 0,0217. 
.B = + 0,03318. 
C = —1,14.10—56. 


Ces valeurs de À, B et C semblent raisonnables 
du point de vue théorique [20]. 

Nous avons pu ainsi calculer une nouvelle courbe 
semi-empirique donnant À} en fonction de Z qui 
est tracée sur la figure 5. 

On peut comparer sur le tableau II les valeurs 
de R4 pour Z = 20, 30, 40, 50, 60, 80 tirées de 
notre courbe avec celles tirées de différentes courbes 
données antérieurement par d’autres auteurs. 


TABLEAU II 
VALEURS DE RK TIRÉES DE : pour Z — 20 30 40 50 60 80 
4 [2] Burhop (1935) et Pincherle (1935). Courbe théorique non relativiste 0,13 0,41 0,685 0,84 0,91 0,97 
(3] Massey et Burhop (1935). Courbe théorique relativiste HS Un, 41 0,68 0,82 0,87 0,93 
[22] Tellez-Plasencia (1949). Courbe empirique 0,14 0,43 0,69 0,83 0,88 
[4] Burhop (1952). Courbe semi-empirique 0,125 0,42 0,695 0,85 0,92 0,975 
[8] Broyles, Thomas et Haynes (1953). Courbe empirique 0,17 0,475 0,73 0,845 0,89 0,945 
[20] Burhop (1954). Courbe semi-empirique 0,141- 0,456 0,708 0,840 0,903 0,945 
(5] Bergstrôm (1955). Courbe empirique 0,17 0,475 0,695 0,84 0,895 0,955 
[23] Gray (1956). Courbe empirique 0,16 0,47 0,715 0,83 : 0,905. : 0,96 
Nos résultats (1956). Courbe semi-empirique 0,138 0,442 0,698 0,833 0,898 0,946 
pour Z — 18 36 47 
16] Rubenstein et Snyder (1955). Calcul non relativiste 01960-67085 


Malgré des valeurs relativement différentes pour les constantes À, B et C, notre courbe est assez proche de celle de 


Burhop (1954) [20]. 
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1,00 he — à 
L | BA x 8 
RUE M Zi 
| 4 | 
RARE 
de F2 
0,75 
Bi EE PNA de 
À 27 cet 
L:) 
Bien v 3 x 12 
| é # 218 
REA ie RE o 5 x 14 
0,501! | sl > 6 + 15 
PA RIRES eo 7 x 16 
LÉ ÉRPEERES PRES a 8 MT 
ÉTAT ET 4-9 . 18 
(ME D LE Ur VW. 
( kk ÿ=A+BZ+CZ° 
L- DE LR, 
DS EE A=-0,0217 
FAN LES LA BU 
| B =+0,03318 
RE) Lt APN C2 1e. 10S 
LATE 
0 25 5 75 2 100 
Fic. 5. — Rendement de fluorescence de la couche X. 
1. (Spectromètre 8). Srerren, Huger et HumBeL, Helv. 10. (Spectromètre GB). AzumA (T.), J. Phys. Soc. Japan 
Phys. Acta, 1949, 22, 167. 1954, 9, 443. 
2. (Spectromètre $). BercsrrôM et THULIN, Phys. Rev, 11. Roos (C. E.), Phys. Rev., 1954, 98, 401 (mesures 
1950, 79, 538. directes sur un faisceau de rayons X avec un 
3. (Émulsions photographiques). GERMAIN (L. S.), Phys. compteur à scintillations). 
Rev., 1950, 80, 937. | 12. (Spectromètre 8). Prusrr (C. H.) et WiLriINsoN 
&. CR proportionnel). Inscx, Phil. Mag., 1950, R. G.), Phys. Rev., 1954, 96, 1340. 
5. (Compteur proportionnel). Wesr et RoraweLr, Phil. 14 BEN Ê}:; WAPSTRA PASSE LES 
Mag., 1950, 41, 873. A CR DRE One n 
6. (Spectromètre B et coïncidences). HuBer, HumBEL 14. (Compteur proportionnel). Harrison, CrawronD ef 
É ; ; * Hopkins, Phys. Rev., 1955, 100, 841. 
ScHNEIDER et de SaLir, Helo. Phys. Acta, 1952, ï ; 
25. 3 : 1951, 24 627. 15. (Compteur proportionnel). Heinrze (J.), Z. Physik, 
7. (Compteur proportionnel). BerToLINI, Bis et Zappa, 1955, 148, 153. 
Nuovo Cimento, 1953, 10, 1424. 16. (Spectromètre $). Mzapsenovic (M.) et SLaris (H.), 
8. (Spectromètre 8). BroyLes, THomas et HAyNES, Phys. Arkie Fysik, 1955, 9, 41. 
Rev., 1953, 89, 745. 17. (Spectromètre 8). Gray (P.), Phys. Rev., 1956, 102, 
9. Lasxar, C. R. Acad. Sc., 1953, 236, 2149 (valeurs 1306. 
semi-empiriques tirées des mesures de largeurs de 18. (Spectromètre 8). LABERRIGUE-FRoLOw (J.) et Rap- 


raies X de Gokhale, Ann. Physique, 1952, 7, 852). 


VANYI (P.), C. R. Acad. Sc., 1956, 242, 901. 
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Il est intéressant de constater que près des deux 
tiers des valeurs expérimentales sur lesquelles nous 
nous sommes basés n’avaient pas pu être incor- 
porées par Burhop [20] dans ses calculs, car elles 
ont été publiées ultérieurement. 

_ La courbe que nous obtenons présente un 
“maximum pour Z — 98,5 ; R% atteimdrait ainsi 
- une valeur maximum de 0,956. 


* 

* *% 
On peut enfin examiner les modifications 
apportées par ces nouvelles valeurs du rendement 


de fluorescence X dans les déterminations [1], [24] 
du rapport d’embranchement capture ZL/cap- 


Hture: K — re les expériences faites sur 7°Kr et 


| 71Ge avaient donné respectivement avec des 
valeurs de RÀ4 de 0,57 et 0,45, des rapports AL /AK 
de 0,26 + 0,03 et de 0,30 + 0,02, alors que la 
théorie ne donne que 0,10 dans les deux cas. Les 
nouvelles valeurs sont pour 7°Kr : Rz — 0,586, ce 
qui conduit à AL/AK — 0,22 + 0,03 et pour 
71Ge : Rx — 0,473, ce qui conduit à AL/1K — 
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Il nous semble, en conclusion, que la précision 
obtenue pour les valeurs de he n’est pas encore 
tout à fait suffisante devant les nécessités de 
certaines expériences de physique nucléaire. Mais 
nous pensons que cette précision peut être obtenue 
avec les appareils et les méthodes expérimentales 
dont on dispose actuellement dans les laboratoires, 

Nous sommes heureux de remercier M. le Pro- 
fesseur Frédéric Joliot pour l’intérêt qu’il a pris à 
notre travail. Nous exprimons nos remerciements à 
M. Chemla qui a bien voulu procéder pour nous à 
l'extraction de %5S sans entraîneur, et à M. David- 
Boyer qui nous a aidé pour la bonne marche du 
spectromètre 6. Nous avons pu effectuer ce travail 
en tant que Chercheurs du Centre National de la 
Recherche Scientifique et grâce à une subvention 
accordée à notre laboratoire par le Commissariat à 
l'Énergie Atomique. 


Manuscrit reçu le 31 mai 1956. 


Note sur épreuves. — Les résultats tout récents de 
Parronis Jr., BRADEN et WvLy (Bull. Amer. Phys. Soc., 
1956, II, 1, p. 258) et ceux de Roos (qui nous ont été aima- 
blement communiqués par l’auteur avant leur publication 


0,25 + 0,02. nl subsiste cependant un Cart ans Phys. Rev.) sont en très bon accord avec notre 
appréciable entre expériences et théorie. courbe. 
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CARACTÉRISATION DE SPECTRES AU MOYEN DE FONCTIONS INDICATRICES 


Par FLorence AESCHLIMANN. 


Sommaire. — La notion de fonction indicatrice d’un spectre introduite par J. L. Destouches 
est utilisée pour caractériser les spectres usuels, notamment ceux qui se rattachent aux fonctions 
sinus et I. La relation entre fonctions indicatrices de spectres de grandeurs fonctions l’une de 
lPautre est précisée et utilisée. Un exemple de spectre rattaché à la fonction de Bessel J,/3 est 
traité en mécanique ondulatoire et en théorie de la double solution. Enfin la notion de fonction 
indicatrice est étendue au cas de systèmes radioactifs à niveaux quasi-stationnaires. 


1. Fonction indicatrice d’un spectre. — Soit U,o 
le spectre d’une grandeur À, J. L. Destouches [1] 
a appelé fonction indicatrice du spectre toute fonc- 
tion S(€) remplissant les deux conditions sui- 
vantes : 


19 S(€) est une fonction uniformé 


= 0 pour Ce A1Q 
ni s(&) ZÆ 0 pour CeX —X40. 


où X désigne le corps des nombres complexes. 
Quand dans l’ensemble de ces fonctions S(€), on 
fait choix d’une certaine fonction S,(£), cette fonc- 
tion est appelée fonction indicatrice typique du 
spectre. On fait un tel choix quand parmi les 
fonctions S(€) se trouve une fonction usuelle 


St). Une telle fonction permet de caractériser 
un spectre d’une manière simple, et les fonctions 
indicatrices fournissent une classification des 
spectres. 


2. Extension à une grandeur de rang quel- 
conque. — Au lieu d’une grandeur de rang 1 
(c’est-à-dire telle qu’un résultat de mesure 
s'exprime par un seul nombre a), on peut consi- 
dérer une grandeur À de rang quelconque » (un 
résultat de mesure de cette grandeur À s’expri- 
mera par un ensemble de z nombres, coordonnées 
du point représentatif de la valeur de À dans 
l’espace des observations (R,)). Cette fois un 
résultat a de mesure de À est un n-uple ordonné 
formé par nr nombres constituant les coordonnées 
de a dans (R,). Si le spectre de À est discontinu, 
on peut fixer une ordination aux points a de (R,). 
On peut aussi considérer chaque coordonnée de a 
comme la valeur d’une grandeur de rang 1 ; soit 
À; la grandeur correspondant à la 1°" coordonnée 
de a, alors si le spectre de À est discontinu, celui 
de À; l’est aussi et pour À; on est ramené au cas 
des grandeurs de rang 1. Pour À, on définira une 
fonction indicatrice comme une fonction de point, 
après ordination des points du spectre. En parti- 
culier dans le cas d’une grandeur de rang au moins 
égal à 2 qui se décompose en deux grandeurs, le 


spectre de À ; alors 
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spectre de À est le produit direct des deux spectres ; 
dans ce cas on considérera les. fonctions indica- 
trices des spectres des grandeurs composantes. 


3. Cas d’une grandeur fonction d’une autre 
grandeur. — Considérons deux grandeurs quanti- 
fiées À et B et supposons B fonction de À, soit 

B — f(4) | 


et supposons le spectre de À discontinu. Alors à 


toute valeur propre a; du spectre de À correspond 
la valeur propre b, du spectre de B avec 


bi = f(æ). | 
Soit alors S,(€) une fonction indicatrice du 


& = f(G) 


est un zéro de la fonction indicatrice S3 de B si … 
&, est le ième zéro de S,. Dans le cas d’un spectre 
sans point d’accumulation pour À et dans le cas où 
(f(a)l est croissante pour [al suffisamment grand on a 


Sal = te À (1 —) De 


îi=1 


et par suite 


Sn(e) = galo* À (1 — 
i=1 


a) .ehi(®) 


si tous les f(a;) sont non nuls. Pour les f(a;) nuls, 
le facteur sera supprimé dans le produit et rem- 
placé par &“ 

Enfin 
— Ex 


= (1 —") a feet 


si f(0) = 0 
go (£) | 


On agirait d’une façon semblable dans le cas de 
spectres ayant des points d’accumulation. 

On voit alors que S,(f(&)) est une fonction qui 
s’annule pour chaque zéro du spectre de À et 
aussi en général pour des valeurs étrangères si 


f(e) = à] 


4% passe de S, à S3 par la formule 


 No11 
admet plusieurs racines C1, Ga, Gin Valeurs de € 
satisfaisant à l'équation précédente, d’où des fac- 


_ teurs 1) qui seront les valeurs spec- 


Gi 


… trales de À. Donc 


B = f(4) — Sg(f()) = S4(0) Tar(Ÿ) 


“…_ Si fest une fonction algébrique, alors B sera 
&… dite grandeur fonction algébrique de À. 

| « ExemPLre : Soit B — A4?, alors “+a, et — a; 
M. donnent la même valeur b, — aÿ. Si S, comprendles 


(4 C a (0 +97 (© 
2 (: +2). alors S, 


contient le facteur 


/ 
(1— À) one. 
(23 


\ Â/ 


_ | facteurs 


On voit, si + ae &, et — a; cÆX,, que l’on 
SB(C?) — S4(C) S4 (— €) ent) 


Soit G(É) une certaine fonction entière. Nous 
dirons qu’un spectre d’une grandeur À est rat- 
taché à la fonction entière G si l’on peut trouver une 
fonction algébrique f telle que les valeurs spec- 
trales de À puissent s’écrire sous la forme : 


A; = f(G) 


où les €; sont les zéros de G. La fonction S, sera 
exprimable à partir de Gsilon a 


SE) = GIG(E)) G(G(E)... G(C(E)) 


Par exemple le spectre de À sera dit rattaché à 


_ Ja fonction sinus si l’on a 


S4 (€) = sin f(®). 


La multiplicité d’une valeur spectrale de B par 
- rapport à À sera ainsi définie : 19 Si une valeur b 
provient d’une seule valeur a, la multiplicité de b 
sera la même que celle de a. 

20 Si la valeur b provient de n valeurs a; dis- 
tinctes, la multiplicité de la valeur b sera la somme 
des multiplicités de chaque valeur a. 

* 80 Si la valeur b provient d’une infinité de 
valeurs a; distinctes, la multiplicité de la valeur b 
sera infinie. 
Soit alors une grandeur B ; la multiplicité de la 

valeur spectrale b par rapport à l’ensemble des gran- 
deurs sera la borne supérieure de la multiplicité 
de la valeur b par rapport à chacune des grandeurs 
physiques dont B est fonction. 


4. Spectres rattachés à la fonction [. — Les 
exemples les plus simples de spectres sont expri- 
mables par la fonction lou la fonction sinus. La 
formule de Weierstrass [2] (où C est la constante 
d’Euler) 
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1 = tee (145) e * 
P (0) n=1 n} © 


sera la fonction indicatrice 


1 
P (©) 
typique d’un spectre ayant pour valeurs spec- 
trales 0 et tous les entiers négatifs : — 1, — 2,. 

Changeons € en — € dans la formule de Weiers- 
trass. Nous obtenons 


montre que 


sera fonction indicatrice typique d’un spectre 
ayant pour valeurs spectrales 0 et tous les entiers 
positifs : 1, 2, Nous désignerons par NW, une 
grandeur ayant un tel spectre. 

D’après la propriété de réunion des spectres [1] 
nous voyons que 

1 
Ni co 

sera fonction indicatrice typique d’un spectre 
ayant pour valeurs spectrales non dégénérées tous 
les entiers positifs et négatifs, et 0 comme valeur 
spectrale dégénérée d’ordre 2, car de la formule de 
Weierstrass il résulte que l’on a 


sb =— "T0 (1 Ve 


= (1 JU+7- erb 


IT’ signifiant que l’on fait le produit pour toutes les 
valeurs de n sauf zéro. 
Mais la formule des compléments 


II 
et la formule 
T{z +1) = z3.T(:) 


qui devient si l’on change z en — z 
— 3. L'(— 2) = T1 —:) 
nous donnent : 
1 té à 
PO ENG €) 
Nous aurons alors comme fonction indicatrice 
typique d’un spectre ayant comme valeurs spec- 


trales tous les entiers positifs et noue et zéro, 
valeurs spectrales non dégénérées : 


Aie = sin ne 


— 1 


 CTO: TE) 


A NE ta 
CO 2 
=snmt=c. I (15). 
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Nous désignerons par N une grandeur ayant un tel 
spectre. 

Si l’on a B — f(A), f étant une fraction ration- 
nelle, alors une fonction indicatrice du spectre de 
B sera 


58 (Ë) = 54 (Gt) S4 (GES 2 G IE) 


Donc si À admet pour spectre les entiers posi- 
tifs, on aura pour B : 


1 
re GE). TE GE) TC GE) 


Si B admet pour spectre les entiers positifs ou 
négatifs, on aura pour À : 


SAONE 


Sa (©) = Lesin (T 


(1 (6) sin (r CG (6)).. sin (r Gr (6). 
les valeurs spectrales de À étant données par &, 


tels que Ë, = f(n). 


5. Exemples de spectres. — Considérons les 
spectres les plus simples de l’ancienne théorie des 
quanta et de la mécanique ondulatoire : 


a) OSCILLATEUR DE PLANK EN ANCIENNE 
THÉORIE DES QUANTA : 
EN hNNEnR= 10; 


Ce spectre se rattache à la fonction F ; on aura 
comme fonction indicatrice typique : 


er Ce es 
a =— + ee P (1). 
( hy 
b) ROTATEUR PLAN : 
n? h? h? 
AE 2 ne 
EN ane —\Kin? avec  K— SEL. 
On aura 
Se(g = t. (1, 
È nt 2 al 


fonction indicatrice typique sera 


mA ere €) 
“Env 
Ce spectre est donc rattaché: à la fonction sinus. 


c) ROTATEUR SPHÉRIQUE : 


2 


En = pren +1) = Æn(n +1). 


Considérons la grandeur N, ayant comme 
spectre les entiers positifs et zéro. On aura : 


E= K(N+ + N;) 


N°11. 
d’où D ESS 
De 
Alors d’après ce qui a été établi plus haut, 


comme est la fonction indicatrice du spec- * 


eo | 
tre de la grandeur V,, nous avons comme fonc- M 
tion indicatrice typique du spectre de l’énergie : 


1 


HV TE VS 
Sp €) = cos r\/E +1. 


Donc ce spectre se rattache à la fonction I". 


Sz (Ë) — 


et finalement 


d) OSCILLATEUR HARMONIQUE : 
En = Ko(n +3) cb PEN TI 
Alors ici | 
a le (+ re 5 d 
d’où la fonction indicatrice typique 


S (©) = 


e) ATOME D'HYDROGÈNE : Cette fois nous avons 
un spectre ayant un point d’accumulation ; la 
fonction indicatrice n’est pas une fonction entière 
de € mais s’exprimera encore par l'intermédire 
d’une fonction entière [1]. On a : 


CHE ki 
PAR Tale À] 
n° k°? A 


Ici aussi on peut considérer la grandeur N, et 
nous avons £ = — K N°? d'où 


K 
N+ = +4/——> 
+ VE 
et d’après ce qui précède, nous avons comme 
fonction indicatrice du spectre 


FE «(VE 


fonction analytique uniforme qui a £ — 0 comme 
point singulier essentiel et point d’accumulation 
dés zéros qui sont tous négatifs. (£ — 0 est l’ori- 
gine du spectre continu qui comprend tous les 
nombres réels positifs). 

On voit que tous les spectres d’énergie que nous 
venons de considérer sont exprimables soit par la 


* fonction sinus, soit par la fonction l. Pour avoir 
des exemples rattachés à une autre fonction 
entière, il faut considérer des cas moins connus. 


6. Exemple de spectre rattaché à la fonc- 


Iltion J,,,. — Traitons en mécanique ondulatoire 
|usuelle le cas d’un corpuscule soumis à une force 
“constante et rebondissant sur un plan orthogonal à 
cette force. Soit OZ l’axe suivant la normale au 
. “plan, O étant dans le plan, 


le potentiel est 
V = + Fz; on peut séparer les variables et 


l'on a comme équation de Schrôdinger 


ae F:) dep 
e 9m dz° NT 
Posons pour simplifier l'écriture 
* 2mF 2mE 
D ss AU 


il vient 


2 
+ (Bs + Cu = 0. 


Effectuons le changement de variable et de 


- fonction : 


13 
PMR Cire: Br-E C— Be) 


|" on obtient ainsi comme équation 


d?® no? 4 és 
. RE — 5)? -0 


| | et l’on reconnaît une équation de Bessel d’ordre 1/3. 


La solution générale est 
px) = ais (x) e BJ 173 (x) 


—… Jes conditions aux limites sont que pour z—>0o, on 


ait U(z) &V/Bz. ? —+ 0 et que pour z = 0, on 
D: : CAE 
» ait (0) — 0, d’où pour x, = 5 2 doit 


avoir œ(x,) — 0. En outre la fonction (x) doit 
être finie dans l'intervalle (0, + co) parcouru 
par z, donc o(x) doit être finie dans l’inter- 
valle correspondant. 

La fonction J,,; est régulière, par contre 
Jus a une singularité pour x — 0, done pour 


TB» 0, pour cette valeur d(x) doit 


tre finie d’où 8 = 0 et on aura finalement pour 


SA 


2092 
3B ‘ 


(x) = Jys(æ) et o (x) —=-0 pour x — 


C3/2 

3B 
J,yx ; cette fois nous connaissons immédiatement 
la fonction indicatrice d’une grandeur À, fonction 


est égal à un zéro de 


Cela implique que 
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biunivoque de l’énergie, c’est : 


Salt) = CB = ur 2 


où 7, désigne le n° zéro de J,. 
Le spectre de l’énergie s’en déduit ; en effet si a, 
2C3/2 


3B 


est le zeme zéro de J,, nous avons — 0} 


d’où en remplaçant C par sa définition : 


3 2 pe 
AT 
4 m 


Ey = Aal8 où 


E.\2/8 
om = (+) 


d’où nous tirons l’expression de la fonction indi- 
catrice du spectre de l’énergie : 


£ 
AN ee OntE) 
( A)" 


Cette fois le spectre de l’énergie est rattaché à la 
fonction J,,,. Mais la fonction S, ne peut s’expri- 
mer au moyen de J,,. 

Le même exemple peut être traité en ancienne 
théorie des quanta avec la condition de quantifi- 
cation de Sommerfeld 6pdq = nh. On trouve par 
un calcul élémentaire que 


no AT n)21° 


et 


(ce) 
ral 
= 


Sy (Ë) 


qui fournit des valeurs spectrales assez rapprochées 
de celles de la mécanique ondulatoire. Ces valeurs 
spectrales étant fonctions algébriques de n, ce 
spectre de l’ancienne théorie des quanta se rat- 
tache à la fonction l. En passant d’une théorie 
à l’autre, on change donc de type de spectre 
puisqu'on change la fonction entière à laquelle le 
spectre est rattaché. 


7. Même exemple en théorie de la double solu- 
tion. — Traitons le même exemple en théorie de 
la double solution de M. Louis de Broglie. Les con- 
ditions de quantification sont imposées à l’onde 
physique w qui doit être uniforme et nulle aux 
limites du domaine, puis on se borne à considérer 
son approximation linéaire 


US Up + D 


où 4, est une fonction singulière de l’équation 
linéaire qui tombe à zéro rapidement quand on 
s’écarte de la singularité tandis que v est une solu- 
tion régulière qui seule compte au voisinage de la 
frontière du domaine. Alors la condition de quanti- 
fication appliquée à v ne diffère qu’extrêmement 
peu de la condition appliquée aux ondes à de la 
mécanique ondulatoire usuelle car # étant une solu- 
tion régulière de l’équation d’ondes est propor- 
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tionnelle à d, on a d — Cv. Alors les spectres 
sont pratiquement identiques dans les deux théo- 
ries ; seule l'interprétation change : ici la quantifi- 
cation s'exprime par des conditions imposées à des 
éléments physiques objectifs (fonctions u) contrai- 
rement au cas de la mécanique ondulatoire usuelle. 

Dans notre exemple, le corpuscule est décrit à 
l’approximation linéaire par 


u(z) À uo(z) + v (2) = PJ) + «173 (x) 


où les constantes « et B sont fixées par l'équation 
non linéaire. La fonction J_,,,(x) présente une 
singularité pour x = 0 qui correspond à la cote 
osent. où le corpuscule se tient immobile 
dans tout mouvement quantifié. Cette cote dépend 
de la valeur £,; de l’énergie. On a ainsi : 

ue 
SOI UERS F F 


Si petit mais fini que soit 8 par rapport à «, 
u, + v présente toujours une singularité en 29% ; 
par contre aux limites du domaine, c’est-à-dire 
pour z—> co et pour 4 = 0, la fonction u, + v 
différera très peu de v si l’on à 8  «. Donc le 
spectre d'énergie, sans être rigoureusement égal à 
celui donné par la mécanique ondulatoire, en 
différera extrêmement peu car uw = 8 J_,y 
n’est pas nul pour 4 = 0 mais très faible, il y aura 
une très légère perturbation des niveaux qui 


dépend de Ê et qu’on ne peut évaluer pour le 


B 
(24 
non encore connus. En effet en mécanique ondu- 
latoire usuelle la condition de quantification est 
que pour z — 0 on a (2) = CJ,yx,) tandis 
qu'ici la condition de quantification est pour 


20, u +o = BJ_to) + a ist) = 0 


qui diffère d'autant moins de la précédente que À 


moment car —- dépend des termes non linéaires 


est plus petit. 

Dans le changement de théorie (passage de la 
mécanique ondulatoire à la théorie de la double 
solution) le spectre de l’énergie a encore changé de 
type, comme dans le passage de l’ancienne théorie 
des quanta à la mécanique ondulatoire : dans le 
problème considéré ici, en ancienne théorie des 
quanta le spectre est rattaché à la fonction F ; en 
mécanique ondulatoire, il est rattaché à la fonc- 
tion J,7,: en théorie de la double solution, il 
est rattaché à une fonction entière G(£) dont les 
zéros sont peu écartés de J,/(0), mais cette fonc- 
tion n’est pas la fonction J,,,. On a : 


zéros de [G(£)] —\:zéros de [x Jus (©) + BJ 8 (Ÿ)]. 


Si au lieu de traiter le problème à l’approximation 


1 


+ 
No11 


linéaire comme on l’a fait jusqu'ici, on le traitait 
rigoureusement avec une équation non linéaire, 
« et B ne seraient pas des constantes mais des 
fonctions lentement variables et le spectre de 
l'énergie changerait de type, sa fonction indicatricew. 
n’aurait plus aucun lien avec une fonction de Bessel. 

8. Fonction indicatrice du spectre d’un système 
radioactif. — Quand, en mécanique ondulatoire, on 
traite le cas d’un corpuscule contenu dans une” 
cuvette de potentiel, le corpuscule a une certaine 
probabilité d'évasion. Si cette probabilité est assez, 
faible, le système peut être considéré comme ayant 
des états quasi stationnaires avec une constante. 
d'amortissement : on a ainsi une fonction d'onde 
de la forme [3]: 


î À 
Sn 
UM # = ar(M).et à 


(1) 


On a en somme une fonction d'onde monochro- 
matique à valeur propre complexe, car on peut… 
poser 

Gt 
br (M, t) = ar(M).e” 
avec 


1% 
E=E+SX À 


La partie réelle de la valeur propre 6 donne 
l’énergie, et la constante d’amortissement est « 
proportionnelle à la partie imaginaire : | 


E RE: 1 = 36 


À doit être petit pour que la méhode soit appli- M 
cable. ï 
Le nombre n(t) de corpuscules dans la cuvette « 
est donné par “4. 


He no f, Los (ANNE de 


dry désignant un élément de volume autour du pi 
point M, l'intégrale étant étendue à la cuvette C. 
D’après (1) on voit que 


nt — men [| lar (M)|? drw. 


On peut normer la fonction a; qui est définie 
à un facteur constant près ; alors l’intégrale est 
égale à l’unité et », est le nombre de corpuscules 
dans la cuvette à l'instant initial t — O0 d’où 
n(t) = nyeT*. À est bien la constante d’amor- 
tissement ou encore constante d’évasion. Cette 


constante est l’inverse de la durée moyenne D 
d'emprisonnement du corpuscule : 


+ CO 
A te—At dt 
= 0 
D= + oo 
Le e—At dt 
0 


% 


_ Dest encore la vie moyenne du corpuscule dans la 
‘ _ ee 

- Cette méthode des valeurs propres complexes 
est utilisée notamment dans la théorie de la radio- 
activité « [4]. À apparaît alors comme la constante 

_de désintégration. 

… On peut étendre sans difficulté la notion de fonc- 
tion indicatrice de spectre pour le cas de valeurs 
complexes. Soit un système admettant un spectre 
-de niveaux quasi stationnaires E; avec une vie 
| “ moyenne +, d'où une constante d'évasion 


——. Les niveaux complexes seront par défi- 
x 


nition 
6; = Ej + - À 
d’où 
E,= RE; et x =? 56, 


que ces niveaux et les vies moyennes soient connus 
- par l’expérience ou par une théorie de SR 
cation. 

Dans le cas d’un spectre sans point d'accumu- 
lation, une fonction indicatrice du spectre sera 
définie par. 


D 


C\ ot) 

1S =IH{1—-—)".e" 
Ho 
ne. 

où k; désigne l’ordre de dégénérescence du niveau 
&;, et où g(Ÿ) est un polynôme de degré v(7) fixé 
par la méthode de Weierstrass [5]. S(£) est une 
fonction entière à zéros complexes. 

Dans le cas d’un spectre admettant un point 
d’accumulation &, on aura, d’après la méthode 


[1] Desroucnes (J. L.), J. Physique Rad., 1956, 17, 475. 

[2] DE Broczre (Louis), a) La physique quantique 
restera-t-elle indéterministe ? (Gauthier-Villars, Paris, 
4953). b) Cours oral à l’Institut Henri-Poincaré, Paris, 
1945-1955. 

[8] DE Broceztre (Louis), Passage des corpuscules électrisés à 
travers les barrières de potentiel. (Ann. Inst. Henri- 
Poincaré), 1933, 8, 421-446. 
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fonction entière de 
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d'Émile Picard [6]: 
ALL 

S(&) = (1 


6; — &G 
J @ 3 
—— 7e) .ehjt@), 
o 


6 — & 
La fonction S(£) qu’on obtient ainsi est une 


TR SOIL 


Ed) 
6 — 6 

Si le spectre comprend plusieurs points d’accu- 
mulation &:,,, 62... On décomposera les niveaux 
6; en séries spectrales Gym, la suite des Gym 
pour m fixé ayant 6,,, comme point d’accu- 
mulation. Pour chaque série spectrale on est 
ramené au cas précédent et on aura 


SCIE 


Il 1 
À Gen): 

La théorie des fonctions indicatrices de spectres 
se développe exactement de la même façon dans 
le cas de spectres à niveaux complexes que dans le 
cas de spectres à niveaux réels. 

En mécanique ondulatoire usuelle la méthode de 
niveaux complexes correspondant aux états quasi 
stationnaires n’est applicable que si la constante 
d'évasion est très petite. Mais les fonctions indi- 
catrices sont définies quelles que soient les 
valeurs À, pourvu que les niveaux Æ; soient 
définis. Les fonctions indicatrices de spectres à 
zéros complexes peuvent donc être utilisées pour 
tout système susceptible de se dissocier ou de se 
désintégrer. Elles pourraient notamment être uti- 
lisées en physique nucléaire. 


See 


Manuscrit reçu le 18 janvier 1956. 
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SUR LA DÉTERMINATION DES FONCTIONS D'’ONDES DU CORPUSCULE DE SPIN # 
EN INTERACTION AVEC UN CHAMP MAGNÉTIQUE OU ÉLECTRIQUE CONSTANT 


Par GÉRARD PETIAU, 


Institut Henri-Poincaré. 


Sommaire. — Résolution des équations d’ondes du corpuscule de spin # (méson vectoriel) dans 
le cas où ce corpuscule chargé est soumis à l’action d’un champ magnétique constant ou d’un champ 
électrique constant. Comparaison avec les résultats correspondants la théorie de l’électron de 


Dirac. 


Dans ce travail je détermine explicitement les 
fonctions d’ondes représentant un corpuscule de 
spin # (méson dit de type vectoriel) en interaction 
avec un champ magnétique ou électrique d’inten- 
sité constante. 

Les solutions obtenues fournissent deux 
exemples de solutions complètes des équations 
d'ondes de la particule de spin % dans le cas 
d'interaction avec un champ extérieur. 

Afin de permettre une comparaison entre les 
résultats obtenus ici pour le spin % et ceux obtenus 
antérieurement dans le cas du corpuscule de 
spin */2 je rappellerai brièvement les calculs et 
les résultats correspondants dans le cas de l’élec- 
tron de Dirac. 


I. Champ magnétique constant. 


1. Électron de Dirac dans le champ magnétique. 
HEC HP, 0 

Ce problème a été discuté par de nombreux 
auteurs (M. S: Plessett {1}: 1: D.:Huff:[21 
M. H. Johnson et B. A. Lippman [3]). Je rappel- 
lerai brièvement les calculs de Huff. 

Avec un choix convenable de la jauge, le poten- 
tiel vecteur correspondant au champ #4, = H, — 0 
H, — H peut s’écrire 


AP AR ARTE 


Les fonctions d’ondes représentant l’électron 
sont alors solutions de l’équation de Dirac 


[Po + Pet + (Py + EH) + pas + Mmocald = 0 (1) 
À Ù ps 
Po Ten 


Nous adoptons pour représentation des matrices 
æ, &4 la représentation 


&— OPys a — Oops» 
So OÙ Po = lol; 61 OÙ P1 = |93l, 
C2 OÙ po = [9 |, 53 OÙ Ps — |6.9|, 
Op Og — — ln, PpPa — —— Pr (P, gr = 1, 2, 8). 


On voit immédiatement que Py Py P: SOnt 
intégrales premières, 


Introduisant les constantes W, 8, y telles que 


pod = T0, pb = 180 pad = ny, 


De, y, 2,0 = peer 7 TA 7? 
l'équation (1) s'écrit 

7 (2) 
É Peer ARE D) ane RYAs + morts] dx A0 


Posant y = SY’ et déterminant S par les 
conditions 


SES — 1, Sans = 01, Ska —"œ, 
S*(%yas + Mocu)S = pd 


0€ 


m 
ui = ui + v? (#0 = me), 


— ee 


2u 


soit 


nous sommes ramenés à l'équation 


W 
[= + Pas + (8h + eHx)o + is] d' = 0. 


On peut donc dans (2) poser y — 0 sans perte 
de généralité. 


Si nous décomposons l’équation (2) en deux -! 
équations en écrivant 
qi) pour Ya Vas (0) pour Ÿs et Ya 


(j = 1, 2), nous obtenons le système 


+ moe] 6 + (pa: + (Br + eHx) où] o%) = 0. (3) 


le 
LE — moe] 90 + [pes + (Br + He) 03] ef 


I 
ee 


Nous en tirons 


gr —— lpass + (En + eHz) o] 6, 


w? 
[— mot? — [p2 + (Pr + eHx)? + aaHo;]| où = 0, 


Posant 


W mo cH 
7 0 EE NET ir (4) 


Non 


| _ nous écrivons ces équations 


L 5 
lu) (x) = — Hu (0251 —— i(B + Ax)o>] œ2), : 
- (RE Ra) a be k5,| pen = 0 


En effectuant le changement de variables 


DFA 


| 11 
fl 7 
| d'où À — QUES VE 2 te + kx], il vient (6) 


|: | V2k d ne 

@ ss Eh ee 2 
| Merle el. 
We CE CNT SR Le 
a lent 0 


L’équation (8) déterminant @?(X) est du type 


de Weber 
É Fe “ 
f+(v+s-r)u=0 () 


dont les solutions sont les fonctions du cylindre 
“ parabolique ou fonctions de Weber-Hermite 
BAOO. | oo. 
- Ces fonctions sont étudiées en détail dans divers 
ouvrages de mathématiques et notamment dans 
- Jes traités de F. G. Tricomi [4] et de H. Buch- 
holz [5]. 

Les fonctions D,(x) satisfont notamment aux 
relations de récurrence. 


(È FE S)Dvtx %) Pi: 


(2 —$)psta = — D, x). 


(10) 


Exprimée au moyen de fonctions hypergéomé- 
- triques confluentes la fonction D,(x) s'écrit encore 


Mn re 1 v 4 +? 
=Vr22e Le" 2° 7° :) 


Dir) 
1 
2 
no) Lee 3: a? 
5 æ PF 9 r-109 À :)] (11) 


E(3) 


Asymptotiquement, on a 


| Dte) se mt 


Introduisant les représentations explicites des 
. matrices 61, 62, 63 et les relations de récurrence 
ci-dessus, on obtient immédiatement les solutions 


x os (1 v + :_%)u 1{+ 0 (07h US): 
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(2) Fa da me C; D &—U? (X), 
2 
pt) = Va ES CoDor us, (X), 
2% 
EVORENO AT. 
ob GX : ; 
(& — 1) 5:43) 
= — = C: De 2 x 
iV2Rk À _ a (#) 
DONS 1 JR € 
q) FEES 2 
PS Vo Le | 
i V2k 
_ Ci D ous (X). 
© Te do 1 . A ) 
EN) T== k AS 
SI it sont des entiers négatifs les 


développements des fonctions hypergéométriques 
sont limités et l’on a alors soit v — 2n, soit 
v = 2n + 1. La fonction D,(x) est alors associée 
au polynome d’'Hermite par la relation 


x? 


a “ æ? 
Dh (x) = (1 6 au 2 )=eT Hot 


Cette condition conduit à une quantification de 
l'énergie W par la relation 


w?— hu, = 2kn, 


W? | 
es mèc? + 2eñHn. 


(14) 

On trouvera dans les mémoires cités [1], [2], [3], 
une étude détaillée des grandeurs associées aux 
fonctions d’ondes ci-dessus. 


2. Corpuscule de spin * en interaction avec le 
champ magnétique 77, — 4, = 0,H, — H — Cie, 

Le corpuscule de spin % peut être représenté 
selon plusieurs formalismes équivalents. 

Dans le premier ou formalisme broglien [6], il 
est représenté par les 410 fonctions ©; 
T1 Lo = À, 2, 3, 4, symétriques par rapport aux 
indices 14, te (Di, = Di), solutions du système 
des 10 équations aux dérivées partielles cons- 
truites au moyen de deux systèmes de matrices de 
Dirac 


[P03 ACAPRACAEA + (CAE (œo)istal AT 
ze (À 5 [( EAN + (ali 2) 


+ moc(ay)in (ais) Du = 0. (15) 


Ce système ne détermine directement l’évolu- 
tion que de six des dix fonctions d’ondes. 

Si l’on élimine les quatre fonctions non direc- 
tement évolutives [7]. il reste un système de 
6 fonctions que nous écrirons 


UD, k = 1,2, 


958 


solutions du système 


1 
Pilosim sue (L Æ pxlem + | dj — 


2mM0€ 
louem (SNL REZ a TE 
_ 8 )s + De (+ slam (AS) 0 (16) 
avec 
nn LU eut 0 
(SE 4 0! P3 —= d'A 
s 110000 & 1 0 1 cD 
LE) AT RA 2= = |—i 0 | 
V2]0 4 0 VO À 0 25200 
4 
S3 = 0 
— 1 


Si l’on remplace les 10 fonctions ®,;, par des 
combinaisons linéaires soient 


= — (Raghisi, Disiss 


dy 5 (Raphisis Disiss 
JC, = [Rapatqalisis Disiss Ep = — [Rapota) Uais Diséss 


(17) 
R désignant une matrice telle que 

at = RayR tt, RF=R=R, 
(P, g,r = 1,2, 3), 


le système (16) est remplacé par le système tenso” 
riel équivalent 


> ee > 
PA — PV — my& = 0, 


ie = mel 1à) 


> > > (18) 
P, &-+ (P de) — Mmoc & = 0, 


UN 
v= 1 (B&). 


> + 
Dans ce système seules les 6 fonctions & et & 


— 
évoluent directement. L’élimination de 4 et 7: 


conduit à un système de 6 équations équivalent à 
(16) soit 


L — moc & = 0, 
moC 
= + > 
P,è— jp + men à + 2 P(Pà)+ (19) 
mo mo 
ARE 
BE (HA) = 0 
mo 


Nous allons déterminer les solutions de ce sys- 
tème dans le cas où 
HS = Hi =0;)H} =H = te, 
Nous poserons encore avec un choix convenable 
de la jauge 


An = 0,0 A — Hz, VA, =0 
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IN 11ù 


Po Py et P. sont encore intégrales premières et. 


nous Re écrire 
W> Er R POEL 1 ÉEe 

à = Éd 

Po c ’ Po c , | 
> > _ Le 20 

Py À = BE, Puy Ë = mé, ). 
a > = > 

PaX—=IYA,  p, É =nYÉ 


On voit que l’on peut sans perte de généralité” 


_. y—= 0 [en remplaçant W par 
= VW? ÿ?c?%?|. On a alors 
î 
ct Six EU à 
see tee Aug Lo See Wen À type ou ail 


Si nous tenons compte de la relation 
> > Ne - a - 
PP?— prP—_9oÿ{f1 P)— ne Rot À 


l’élimination de à entre les HR équations (2% 
nous donne l’équation 
FE — P2 — mice] È & ii te A à) 


4 GDF) -0 


e2) ï 


Séparant les équations relatives aux différentes" 
composantes &,, &,, &,, &, » En G,, nous obte- » 
nons successivement 


19 Pour les composantes Œ&,{x) et &,(x), 


Ta, = — moe 6, 
72 (23) 
EE ue 4 


Ww2? : À J dre 
[—rP- méc?] Ge — HE, + rs à a +) =0 | 
; ne M2 
Ge is LUE, &)= 0: 
Ta - = LR ) + met 8] 
Nous poserons 
GE É, DA, = 5) 
Us —= Ge — 16y, & = BE — A, 
d’où nous déduisons 
(FA 2) 1 , ; 
P 6) — 5 [(Px — iPy) ui + (Px + 1Pyluol. 3 


D'autre part 
(Pe+ iPy) (Pa + iP,) = P3 + Pi+ieH, 


Le système (24) s’écrit alors + 
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ñicH 


' cH \ à D? X2 wo? — 8 ue 
LS, DARÉPEE NE} 2 
= met — (1 . à) (P oi u, En 
4 ñeH 0} ù X1? 
ne deu. 9/00 î D | no 
icH La première de ces équations est du type de 
he) pe H)] Rey pe 
eE al Us Weber-Hermite (9). Nous avons donc immédia- 
ñeH NC tement 
D ne cou 0 CE De ie Al (33) 
2k 2 


C, désignant une constante arbitraire. 

Pour résoudre les deux autres équations (22) 
: nous tiendrons compte des relations (10) liant les 
MoCU»x + —— [(P? — reH)u, + (P, — iPj)?u;]. fonctions de Weber-Hermite. 


I 


mocus + mo [CP? + %eH)u, + (P; + iP,)?uo] 


se Nous chercherons une solution de la forme 
résoudre ces équations, no ron 
e ces équations, nous poserons CD een, 4) 
. =, Te = Lo _ = + Vi Va Cr) Ca désignant des constantes à déter- 
n : (28) miner. 
A = Ke de Par substitution, nous sommes conduit aux 
EE relations 
2—= 35 
Nous obtenons ainsi pour les Eros (23), HE JA 59 
où) et (27). | m2 — 2 | 5 Gi 
| RE > jai F1 2k(1 — 6) C1 + 1e Va(V2 1)C2 0 
|: 1 —(B + Rx)? + {u? — pi)] 6; = 0 nn] 5 La 
| 4 [a ni ler (36) 
“A Pom : 
h: . RER tu Nous en déduisons 
| A 0 a: m— pu? À w 
: (à + (B+ ka) w 0 (29) nie (37) 
DE B+kt+R+T “a, k 
‘4 1 + G è 2kC: = (Lo Æ w) (vo EE W + a. Ce, (38) 
6 ù Fe 
pie Fe Eee) tu a Nous obtenons donc deux systèmes de solutions 
qui s’écrivent si nous posons 
are al w?—u? 1 w x? Lu 1 w 
| ONE 2 VE EN , mo = À. (39) 
à; — pa — (Et kx)? — DEA 2k 2 24 < BI 2 2h 
| 0 ï ; 10 uD = CD \} mx) 
à Ge meer ) ve] 2kCr a 
D. à “ uD — re Dvx (X). (40) 
4 WE = Uoüe — es (B + kx)? + #] ur (Ho + #) (u ARR NE ni) 
i 2 HE Re D. 4 
Se (2 = lEtux ke)) un] Et He MD a (A) 
dx es 2kCrx D=(X). (4) 
. . . (72 = ; Va 2 
l k 
Nous introduisons la variable 2 RATE (2e D id =] 
0 
EL, 2 é 3 . 
eee VE (B + kx). (31) Introduisant ces valeurs dans les expressions de 
, Œ, et &, tenant compte des relations (10), nous 
… Nous écrivons alors pour (29) obtenons 
| E Ho rs 
Dé X? © —4u]l. ‘2 X (Æ) fes Co D &°—uto __ 1 (X), (42) 
Las HT of ]s 4 0; (32) + L 
é RE SN RAI O PRE (43) 


>? x? HAE 
[ bo Ja ei 


DAME ONE 0c) 2 
À à à > 12 an = + ul (X), an re un Ge (44) 
+ Mr 5 ua = 0, aan Re ul (X), an + ul) (X). 
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Les conditions d’uniformité et de régularité 
des fonctions d’ondes conduisent à introduire des 
nombres entiers ñ6, 1, n, tels que 


19 Pour l’onde &,, 6, 


pè—u5 1 
Aer 45 
20 Pour l’onde &, u® 
nr EU 
— = 2. 46 
EYA Ro U NT (46) 
3° Pour l’onde AM, 
HP AUG dal w 
= 2. 
2k F9 Jo No + (47) 


On en conclut que si les 


valeurs de l’énergie 
sont les valeurs quantifiées 


1 
PRET LEE (ro ne) 


la solution acceptable des équations d'ondes est 
l'onde &,(X), 6, (x) avec A, = À, = E, = 6, = 0. 

Au contraire si les valeurs de l’énergie æ sont 
telles que : 


k 8 

* de Ù de A 
w rs HR Are 
ou 


k 3 
2 er — — |] 
aura uè + 2k (ns +2) 
les ondes acceptables sont les ondes &,, &,, &,, 


&, correspondant respectivement aux fonctions 
ŒD, ut ou AD, avec e,; =0, = 0. 


Nous compléterons cette étude par le calcul des 
fonctions d'ondes non évolutives Ÿ et 9€. 
Nous avons 


(Bee 


[Pr — iPyjus + (Pr + 1Py)uo] 


” 2m0C 

Avec les expressions précédentes de , et de 
u, les relations de récurrence (10) entre fonctions 
de Weber, nous donnent immédiatement 


V = ne [to + #] CD, (48) 
 V2RC 
TRANS EIRE ET DURE 
+ } — 
[ue w + | 
On obtient done 
Pour l’onde (I) 
, V2r Cw 
Var = — TE Di 1. (49) 
Horaires 
Ho 
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Pour l’onde (II) 


i V2k Can 
k 
— + = 
Fo 
Nous allons maintenant évaluer les 
A { 
ae — (PA à). 


Non avons 


Van D\ 


Re Te lt), 


er 
y, =, = PE s 
mu y mc aa (x) 


ee 7 (Pay — PyA). 


Pour les ondes 4€, et dt, en posant 
m — p.? 
2k 


WA = HLoË:(X), 


{ 
Vo —= GATE à partir de 
& = CoDy(Æ), 


nous obtenons 


PATES ONE TP. 
Fe 7 2 


= 14/É cu 107 +2, PER 
1, ER on 
= /E cn tx) 20 0e 
Pour l’onde 4£,, on a successivement 
3e: = = | È BRU) al (54) 


V2k 
Var C 
= 7 x Dv 12) 
LE WE — 
Ho Le 
Nous avons donc 
pour l’onde (I) 
jé SRE F Di (4) (55) 
Bo + # + Lo 
pour l’onde (IT) 
peu = VRON, pen (56) 
NE CPC 


II. Champ électrique constant. 


|. 1. Électron de Dirac dans le champ électriqu® 
constant E, — FE, FE, — E, — 0. 

L’étude des fonctions d'ondes de l’électron de 
…— Dirac en interaction avec un champ électrique 
I constant à été développée il y a longtemps par 
Î À F. Sauter [8] et S. Szczeniowski [9]. Je rappellerai 
+ brièvement la solution donnée par ces auteurs. 
L’électron de Dirac est représenté par les solu- 
tions de l'équation 


+ Pats + Mocaa] Lx, y, z,t) = 0 (1) 


Sat 
Po = — À 


ot p=ü— 


dx 


Û Le potentiel scalaire V est une fonction de la 
+ variable x seule, V — V(x). Dans le cas du champ 
… électrique constant V(x) = — Ex. 

é L’équation (1) admet les intégrales premières 


poŸ = 4 pub = BU piÿ = vid, 


et nous écrivons 


À ort-prv-vr0 


px, y, 3, t) = Yixjef 
. (x) est déterminé par le système 
M [7 + ePtc) + peus + Bas + via + mea] de) = 0 (2) 
On montre facilement qu’il existe une transfor- 
mation S telle que d(x) = S4'(x) et que 
SES EL OT Se 4 
18" (Bue + Yas + Uo%)S = tou 
| avec : f 


Uo = moc/h, = ui + B? + Y?. 
La fonction d’onde 4’ est alors déterminée 
nn par l'équation ô 
| n 
LE + ee) + pou + wa] do = 0 (9) 


à Dans ce qui suit nous écrirons d(x) pour d'(x). 
… Avec la représentation adoptée pour les matrices «, 
ce système s’écrit 
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Posant 
Up = Yi — M; 
en Ve Tee Ÿs, 


Ui = Yi + da, 
Mir Vo de Qi , 


ce système s’écrit encore 


(5) 


[re (© + ep) Ua — us = 0 


et les mêmes équations entre fonctions v, et 0. 
Il nous suffit donc de résoudre le système en w,, ua. 
On en déduit immédiatement 


qu = [re (E + ev)]u 


[rs 2 (5 +ep)" ire grad V + utn]us 20 


1) 
[r? — (F— Be) + iecE + un] Us 0: 
Pa QU — k, il vient 
ñc h 
qu = fi — 16 — ti] (8) 


d? : 
LS + Ge — ka? — y? — it] Uy — 0. 
Si nous effectuons le changement de variable 
w — kx — ÀX, cette équation s'écrit 


Ê 4 12 : 
+ Htu+iblu 0 (0) 


Si lon détermine À par la condition 


ae F/ 
E = EE À — . 
F % À 5 k, oT 


nous obtenons 
2 ir 
Der 7e 2l#w— ka), 


D? X2? ip? 

ae 
Cette équation est encore du type de Weber- 

Hermite et nous écrivons la fonction w; 


;| NE (10) 


u(X) CD iut (Z). (11) 
à T2 
On en déduit l’expression de w1 : 
ik [D X m 
HE 55 GS X 12 
UE Dre 53 CD {ul } > °(42) 


962 


ou encore 
se HÉMENTE : 
Hess 40 D_iw_; (à). (13) 
On écrirait de même les fonctions #, et P2. 
2 = CaD Le (4) 
= er, one ICE (14) 
1 V2} 2: 1 ) 


Les propriétés des fonctions du type de Weber 
à indice complexe telles que , 


. 9 que 
Us — Ci Dr w[ 2e E ww — ka) | 


ne sont pas évidentes. F. Sauter les a étudiées 
d’une façon très complète. 
Si l’on détermine À par la relation 


X 1 TR 15 
ee x UE (15) 
on obtient au lieu de (10) l’équation 
MAT Ua 
+75 ua(X) = 0 (16) 
Cette équation est du type 
: 
f + (Fa) 0 (19) 
étudiée notamment par C. G. Darwin [10]. 
Si l’on détermine À par la condition 
12 
sue À= Vh: (18) 
on obtient l’équation 
d? He 
> + X? — (5 + à) [us 0: (19) 
Cette équation est du type : 
+ (x? + a)y = 0, (20) 


dont les solutions ont été étudiées très en détail par 
C. P. Wells et R. D. Spence [11]. 

Ces études mathématiques permettent de déter- 
miner complètement les caractères et les valeurs 
numériques des fonctions u, et us, v, et v2. 


2. Corpuseule de spin À en interaction avec le 
champ électrique E(x) dérivant du potentiel sea- 
laire V(æ). 


Les équations (19 I) s’écrivent ici 


[Po + eV(x)] à = ae ni &) + mic? à) 
(21) 


era en = [p%+ mc] à GG à) 
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> > Li 
Si nous désignons le couple (&, ë) par ÿ nous | 


avons encore les intégrales premières 


pod = T4, prb ag pa = ny4, 


d’où 


; (Wt-nBu— ya 


ÿ = Y(rle 


di 
Nous pouvons encore sans perte de généralité | 
poser y — 0, mais non pas à la fois 8 = 0 et. 


y = 0 et ceci introduira une différence importante | 
entre les cas du méson vectoriel et de l’électron. 


Nous séparerons dans le système (21) les ondes 
correspondant aux composantes z d’une part, x et W 


y de l’autre. Nous obtenons ainsi : 
Pour les composantes &{(x) et &,(x) 


[= + ePtaies = moéste) 


W 
[= “ ep|s, à BE [p? L mice 
Nous en déduisons pour &,{(x) 
W 
fe ere née = 
ou encore 
2 2 h2 
BE LfE + er" pa = 


Pour les composantes &,, &,, &, 6, nous 


obtenons les équations 


mo F F eV (a)] Ào = [ps + MC lEz + ABPrÈy 


W 

me [5 + eV(a]Rs = nBpeés + 1%? + mets, 
W 1 
moc [= + eV()] En = [n°6 + mc Ar — 16paly # 
(25) > 
moC Es + eV()] Ey = [pi + mic]Ay — 28pel ge ; 

Si nous posons 
À = n°8? + mic, 

ui = Az + 16y, D, = Ay + 16, (26) L 


Us 16 Re 


les équations (24) et (25) nous donnent par combi- 


naisons 

Au; = [r6p2 + imoc (F nA cv)]», 
2 
Aug = [rêre — imoC (+ cp)]es. (27) 


(pè, + méc?lv, — [ imoc Fe ie cp) à 18pal =: (8) 


LE + mets = [moe (5 + 6) + aBpel ur. 


N°11. 


Dial ne: ht reléechhe ice éd. à 


(29) à 


4) à 


9 5 à . Ê . 
—. L'élimination des fonctions w,, u, s'effectue sans 
difficulté et on obtient pour les ondes »,, v, les 


= — ( E ep) + mêc? + re] HE) =) 
| (29) 


— ee + ep) + mêc? + ne] Do (Lt) 20: 


no moe : ; or 
. Posant u, = ——,ces équations s’écrivent encore 


é [È 2e VE  gradP + al + ep] —(ui+62]e =0 


: ne sgradF + À [+ cv] tu + 82/00. 


ro 


(30) 


On voit sur ces équations que l'intervention de 

- la constante B différente de zéro entraîne la diffé- 

rence de nature des fonctions p, et p2. 

__ Jusqu'ici nous avons supposé simplement un 
. potentiel V — V(x). Nous allons maintenant 

- spécialiser celui-ci en considérant le cas du champ 

… électrique Æ constant. Nous poserons alors 


à : AP = EX; 
Écrivant 
Dee cE W pee 
ñ ñc 


EC Les fonctions v,(x), va(t) sont solutions des 
éq ations 
F5 + uw 8 + 89 — Foie = 0 
Le [+ (æ — kx)? — (uë + PR?) + À Va(a) - 0. 
“ Pi et u, sont alors déterminés par 
fu + Blu = à]8 À + vols — koi], 
fuÿ + Blue = if — malw — tons, . 
1 Nous poserons encore 
| w— kr = XX — (33) 
nous déterminerons À par la condition 
dunes 
CTE d’où AVE see 
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Nous obtenons alors les équations 


ose (SE) aa — 0 
Lx F :(É FF) JE] 112) = 0, (35) 
Re (S) + ET rcn = 0 


Ces équations sont encore du type de Weber- 
Hermite et nous en écrivons les solutions 


à &,(X) = CoD ut +6? 1 (X), 
ARS SE 
Di(X) = CD pus 64 1 (4) se) 
LE tal 
v,(X) = CGR que 1 (ZX), 
Se FE 


Co Cas C2 désignant des constantes arbitraires. 
À partir des relations entre fonctions de Weber 


2D,(x) = — Dyyi(a) + vDy=:(x), 


: OÙ 
aDy(x) = vDy= (x) + Dy+1(x), (2) 
posant 
de ua + B\ 1 
ARE PRE 2k ] ao 
RE en il Re A 38 
V1 = i( 3k 9° . (33) 
ea oc pie, B:\ "1 
LRÉEE ( 2k ei 2’ 
nous obtenons 
(39) 


à À À 
SX) = U À 4, (X) = — ge Cove Du) + D,,+1(X)], 


' AU 
Wu + Blu = — EG [ SE — ex] D, 0x, 
1e tu (40) 
Gi + plu, = — CS + de x]D, (0. 
Après calculs nous écrivons ces expressions 
AC 
NT EE 
[to — iB)Dyti(X ) + wllo + i8)Dv,:1(2)1, 
ER 1ÀCi (41) 
HAT ec 
[io + ÉB)Dui+1(X) + Vilto — 8) D, —1(Æ)], 


Rassemblant les expressions de wj, Uy, 03, Vs 
nous obtenons les fonctions d’ondes 


(A) = à li + te] (&2) 


LES À . « 
= qe pl allo — iB}Du (À) 


+ (uo + 18)D,, +1(X)] 
TA Co[lto e 16)Dy,—1(X) He (Bo QE TE 1B)D,,+1(X)], 


Ay(X) = 5 La + 03] = 5 [CaDua(X) + CDw(AN, 
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AUX) = CoDvo(X), 
E(X) = — (1/2) [01 —v» 
SylX) = (2) — ul = 
= [A2{u5 + 8?)] [Cilv 
+ (to + iB)Dy,+1(4)] 
+ Colvto + ÉB)Du1(X) + (to— 16)Dv,+1(Æ)T], 
EÂX) = (M0) ColvoDw—1 + Dve+1le 


Au moyen de ces fonctions nous calculons sans 
difficulté les fonctions d’ondes non directement 
évolutives : 


] re (1/2) [CiDya(A) GER CyDy(X)], 


1(to — iB)Dn—1(X) 


% = (1[ms)(? 8), (43) 
d’où (44) 
VA) | Cf D) — 55 Du] 
1 | 
— [2 Due (2) 2 Due si( |, 
JCg(X) = (B/u0) CoDy (X), 
2,12) = (N/u0) ColDu 1 () — Du (A) 
dC,(X) = (1 [moc) (Px dy = Py Xe) (45) 


F pe Dé Dy +12) ]|: 


Uo — 18 Wo + 16 
Dans le cas particulier où 8 = 0 ces résultats se 


simplifient notablement. La résolution directe des 
équations (21) s'effectue immédiatement 


> > 
Dans ce cas, posant ÿ — (&, ë), on à 


W ASS Dans les d idérés, ch étique. 
AT, DR 0, ÿ — den PE don (£6 ans les deux cas considérés, champ magnétique 
D GATE AE ES is A A pe constant, champ électrique constant, le calcul des… 
) [F ; Le AE dix fonctions d’ondes peut donc s’effectuer 
Fe RE eV{(a)] La e 1pe 4e 0 Et) complètement. Comme dans le cas de la théorie de, 
(47)  Dirac ces fonctions d’ondes s'expriment par des. 
FE + e V(a)] Exlx) = moc A(x) combinaisons de fonctions de Weber-Hermite, 
je mais la théorie du méson introduit trois fonctions 
W co 7 principales quand la théorie de Dirac en introduisait 
Fe Æ e V{x)] CAMERA ENTER deux, ce qui correspond aux nombres respectifs 
2 (&8) des états de spin. 
PAR RS 2 
Fe, Pie) | Eu) MC [p£ + mé c'1A (x) Manuscrit reçu le 28 février 1956. 
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[Wie + eV(a]lGs(e) = (1/m6)[pà + mêc?] À ,(a). 


Nous en déduisons immédiatement pour 


u(x) Sr Glx), AXy(x), X, (x) (50) 
[[WJe + eV(x)]? — pà — mi c?] u(x) = 0. 
Pour V(x) — — Ex, effectuant encore le change- | 
ment de variable 
(w — kx) = XX, A2 = 1k]2 


nous obtenons l’équation de Weber-Hermite. 
Dore ip ]94 (A 0) 
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Gx(X) se C D,,(X), 
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EXPOSÉ ET MISE AU POINT BIBLIOGRAPHIQUE 


ÉMISSION ÉLECTRONIQUE DES CORPS SOLIDES 
APRÈS UN TRAVAIL MÉCANIQUE OU UNE IRRADIATION 
(EXOÉLECTRONS & PHOTOÉLECTRONS) 


Par MM. E. Lours HUGUENIN et JEAN-GEORGEs VALAT,. 


Sommaire. — De nombreux métaux et isolants peuvent émettre spontanément des électrons, 
après avoir été soumis à un travail mécanique ou à une irradiation. L’appareil de mesure généra- 
lement utilisé est un compteur Geiger-Müller du type « cloche ». Les résultats montrent que cette 
émission, très faible, décroît dans le temps suivant une loi approximativement hyperbolique, 
fonction des traitements subis par l’échantillon. Une élévation de la température a une très grande 


influence sur l’émission : 


que l’exoémission provénait d’une réaction chimique (oxydation). 


elle fait apparaître de nombreux maxima. Certains auteurs ont pensé 


Une théorie récente due à 


Nassenstein établit un parallélisme entre l’exoémission et la phosphorescence. Cette théorie élec- 
tronique, qui permet d'expliquer la plupart des phénomènes observés tant pour les métaux que 
pour les isolants, semble vérifiée par les travaux récents de Bohun. 


Introduction. — Les métaux et certains isolants 
peuvent émettre spontanément des électrons à la 


… température ambiante, après avoir été soumis à un 


- limage, écrouissage, étirage, etc. 
20 Une irradiation : 
quelques matériaux, avec la lumière ultra- violette ou 


« traitement mécanique ou à une irradiation. Ils peuvent 
» aussi donner lieu à une émission photoélectrique même 


quand la longueur d’onde de l’éclairement est supérieure 


… à la longueur d’onde du seuil. Cette émission est en 


général] très faible et décroît dans le temps suivant une 


. loi approximativement hyperbolique. 


Les traitements qui favorisent une telle émission 


sont les suivants : 


19 Untraitement mécanique : polissage à l’ émeri (°), 
[17], [21 à 29]. 


avec des rayons X, et pour 


même visible, 3° Un bombardement corpusculaire : 


… électrons [45] ou ions [17]. 4° Une attaque chimique : 
… acide [8]. 5° Au cours d’une transformation allo- 


tropique, on peut observer également une émission 
électronique [8], [221, [32]. 
L'émission spontanée avait été observée par de 


— nombreux auteurs au cours de manipulations avec des 
. compteurs 


Geiger-Müller, démontables, du type 
« cloche » ; dès la mise sous tension de leur appareil, 
ils avaient constaté le passage de particules ionisantes 
dont l'intensité décroissait dans le temps. En pro- 
tégeant leur installation contre les rayons cosmiques 


et les photons, Wiedenbeck et Crane (1949) [491] ont 
montré que lPémission provenait uniquement de la 


popode récemment polie de leur compteur. 


) Pur l’émeri, 


{1 voir les réserves faites par Lohff et 
Raether [34]. 


Kramer (1949-1955) [21 à 29] a étudié systéma- 
tiquement l’émission des exo et photoélectrons. Il a 
pensé que l’émission était provoquée par des phéno- 
mènes exothermiques (changements de phase) à la 
surface du métal d’où le nom : exoélectrons. Haxel, 
Houtermans et Segger (1951)[17]ont pensé au contraire 
que l’émission provenait des réactions chimiques super- 
ficielles. Mais ces deux théories présentent de nom- 
breuses lacunes. Récemment Heinrich Nassenstein 
(1955) [38] à exposé une théorie analogue à la théorie 
des centres F et de la phosphorescence ; il a pu 
ainsi interpréter la plupart des phénomènes observés, 

Par ailleurs l'effet Kramer est certainement res- 
ponsable de deux phénomènes particuliers qui jusqu’à 
présent n’ont pas reçu d'explications satisfaisantes : 
l’effet Russel et l’effet Shaw-Russel (1897) [41] et 
Shaw (1948) [47] ont montré que les métaux récemment 
polis ou coupés provoquaient des réactions chimiques 
(formation de peroxyde d'hydrogène, par exemple). 


CHAPITRE I 
APPAREILLAGE UTILISÉ 


1° Comptage de l’émission. — L'émission des exo 
et des photoélectrons étant très faible (tout au plus 
quelques milliers d’électrons /minute /cm?, au bout 
d’une minute) il est nécessaire de disposer d’un appa- 
retllage extrêmement sensible. L’organe de mesure à la 
base de tous les montages est un compteur Geiger- 
Müller du type « cloche » où la cathode est constituée 
par l’élément à observer. 
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On peut chauffer cette cathode pour étudier l’in- 
fluence de la température. Des montages très perfec- 
tionnés ont été imaginés par Kramer (fig. 1) et par 


ee 
VEETTERRENN 
N 


F1c. 1. — (D’après Kramer [25] et [28]). 
Compteur à pointe pour l’étude des surfaces. 

1, Pointe ; 2. Diaphragme interchangeable ; 3. Raccord 
de l’amplificateur et arrivée de la haute tension ; 4. et 5. 
Lentilles de plexiglass ; 6. Glissière pour les filtres optiques ; 
7. Source de lumière ; 8. Vis de positionnement (trois) ; 
9. Ouverture pour l’essai d’étirage ; 10. Vis d’étirage ; 
11. Microscope de lecture ; 12. Déplacement du microscope ; 
13. Comparateur. 


Haxel, Houtermans et Seeger (fig. 2) [17]. Le montage 
de Kramer facilite en particulier la mesure de l’émission 
des photoélectrons (il est décrit en détail dans [28]). 


220V 80W 


Chauffage 


Circulation 
d'eau 


ST Vers 
l'ampli la pompe 


Fic. 2. — (D’après Haxel, Houtermans et Seeger [27]). 

Dispositif d’essai : Le tube de comptage possède un 
cylindre amovible à la base duquel on peut attacher la 
feuille de métal passée à l’émeri. On peut faire varier la 
température de l’échantillon par un dispositif de chauffage 
pouvant s’emboîter dans le cylindre. 


Remarque. — Lohff et Raether [33], [34] ont utilisé 
un montage différent pour étudier les échantillons 
métalliques sous vide :les exoélectrons sont recueillis 
par un photomultiplicateur à la sortie duquel se trouve 
un amplificateur suivi d’un système de comptage. 


29 Enregistrement photographique. — TL'’enregis- 
trement photographique est basé sur l’effet Russel, 
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Lorsqu'un métal récemment poli est en contact avec 
de la vapeur d’eau et de l’oxygène (cas de l’air ambiant) 
les électrons quittant la surface (effet Kramer) pro 
voquent la formation de peroxyde d'hydrogène qui 
réagit sur la plaque photographique posée sur le métal" 
Les plaques utilisées par les chercheurs sont des 
plaques « Ilford Q 3 » [8], [13], [14] ou des plaques, 
« Schumann » [2], [12]. Les photographies obtenuess 
permettent de déterminer la position des zones très. 
émissives et d'évaluer la densité des centres émetteurs 
Remarque. — Kramer [26] a imaginé un montage 
qui réduit considérablement le temps de pose (habituel-W, 
lement, on pose pendant plusieurs jours). Il déplace son 
échantillon poli devant le petit diaphragme de son“ 
compteur « cloche ». Chaque fois que le compteur est 
traversé par des électrons, il allume une petite lampe à 
décharge, dont l’image se forme sur une plaque photo-" 
graphique ; cette plaque suit fidèlement les dépla: 
cements de l'échantillon. Ce montage donne des photo- 
graphies en noir et blanc. La qualité de cette « ampli" 
fication optique » dépend de la finesse du diaphragme. 


CHAPITRE II 
RÉSULTATS EXPÉRIMENTAUX 


L'émission des exo et des photoélectrons a été étudiée! 
en détail par de nombreux auteurs (voir références) qui 
ont obtenu des résultats concordants mais en ont. 
donné des interprétations différentes. | 

Nous mentionnerons tout d’abord les expériences! 
de Tanaka [48], Haxel, Houtermans et Seeger [17], et 
Lintner et Schmidt [32] qui ont montré que les surfaces 
polies ou irradiées donnent lieu à une émission d’élec- 
trons et non à une émission de photons. Haxel 
Houtermans et Seeger ont placé une grille devant 
l’échantillon traité, celui-ci étant porté à un potentiel 
positif par rapport à la grille. Ils ont constaté que” 
l’émission était entièrement arrêtée pour 10 V. 

Cette expérience permet de conclure : 1) que l’on am 
bien à faire à des particules négatives ; 2) que ces 
particules quittent le métal avec une vitesse inférieure” 
à 10 V. 

Une mesure précise de la vitesse a été impossible” 
avec cet appareillage car le champ créé par le filament 
pénétrait dans l’espace échantillon-grille. Lintner et 
Schmidt ont utilisé un montage différent : ils ont. 
accéléré les particules sous 30 kV et leur ont fait 
traverser un écran de 1 mg /cm?, derrière lequel était 
placé le compteur. Seuls des électrons de 30 KV — et 
non des ions —peuvent traverser un tel écran. Cette 
manipulation confirme que les particules émises par 
les métaux récemment polis sont bien des électrons. 
Il serait cependant intéressant de préciser si l'émission 
spontanée des métaux ne contient pas aussi des ions 
négatifs (par exemple des ions O5). 


10 Résultats des expériences avec le compteur Geiger= 
Müller. — 1) VARIATION EN FONCTION DU TEMPS. — 
a) Cas des métaux. — L'émission des exoélectrons 


: 


SEC TLR : : A 
décroit suivant une loi hyperbolique 1 = 7; tout au … 


moins à partir d’un certain temps. Le produit Z X test 
donc constant. Dans le cas de Ja figure 7a nous avons 
It 46,105 coups em?, 
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La loi hyperbolique n’est pas toujours vérifiée : en 


. particulier, dans le cas des métaux alcalins et alcalino- 
… terreux, 1l faut admettre une loi plus générale 


1 = A. (fig. 3) avec a supérieur à 1 pour Na, Ca, 


Nombre d'électrons 
Seconde xcm? 


Fic. 3. — (D’après Lohff et Raether [34]). 
Émission de quelques métaux de différents groupes de la 
classification périodique, après passage à la brosse d’acier. 
En haut, au milieu, lire Ca au lieu de Co. 


et a inférieur à 1 pour Be [33], [34]. Quelquefois il faut 
admettre une loi plus compliquée. 

En ce qui concerne les photoélectrons, l’émission 
déeroît suivant une loi plus complexe dépendant de la 
longueur d’onde (fig. 4). Cette émission apparaît pour 


Nombre de coups parseconde 


04 06 1 
Temps écoulé après Je polissage 


2:46 10 20 40 60100 200min 


Fuel — (D’après Kramer [25]). 
Décroissance des exoélectrons (fig. 1) et des photo- 
électrons (fig. 2) pour le platine poli à l’émeri. 


des longueurs d’ondes supérieures à celle du seuil : 
Grunberg [14] trouve que le seuil du zinc se déplace 
après un polissage. 


b) Cas des non-métaux. — Les principales études ont 
été faites par Kramer [23], [25]. La courbe de décrois- 
sance de l'émission est encore une droite (en échelle log). 
Elle répond à la formule générale Z = B.ft (fig. 5 et 6). 
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… Remarque. — Dans le cas des corps non métalliques, 
il est préférable de traiter les échantillons avec des 
rayons X. 


, Nombre de coups par seconde 


1 ee 


4 6810 20 40 
Nombre de minutes apres le broyage 


60 80 100 


FiG. 5. — (D’après Kramer [23]). 
Courbe de décroissance pour le quartz : après 10 minutes 
le quartz a été porté à 55 9C en dehors du compteur à 
pointe ; puis étudié à nouveau lors de son refroidissement. 
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Fic. 6. — (D’après Kramer [23]). 
Courbe de décroissance pour le quartz : le quartz a été 
maintenu pendant 10 minutes à 26 °C, puis porté à 52 0C 
dans le compteur à pointe. 


2) INFLUENCE DE LA TEMPÉRATURE. — Lorsqu'on 
élève progressivement la température de l’échantillon, 
on constate que l’émission passe par plusieurs maxima 
avant de s’annuler. Un nouveau traitement et un 
nouveau chauffage renouvellent l'émission et donnent 
une courbe sensiblement identique. Cette possibilité 
de régénérer l’échantillon disparaît cependant après 
quelques opérations. 


968 


Les corps métalliques et non métalliques se com- 
portent sensiblement de la même manière (fig. 7b et 8). 
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Fic. 7. — (D’après Haxel, Houtermans et Seeger [17]). 
Ja) Décroissance dans le temps de l’émission électronique 
d’un échantillon d’or. Pour t = 0, celui-ci a été passé à 
l’émeri. 
La ‘courbe continue indique le cours de l’émission élec- 
tronique. Les cercles donnent le produit de l’émission par 
le temps. 


7b) Dépendance envers la température de l’émission 
électronique, lors d’une élévation continue de la tempé- 
rature. 

La courbe en tirets indique l’élévation de la température 
(échelle se trouve à l’extrême droite de la figure). La 
courbe en traits pleins indique le courant émis en fonction 
du temps et par conséquent en fonction de la température. 
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; Fic. 8. — (D’après Kramer [25]). 

Emission exoélectronique du gypse en fonction de la 
température. 

1. et 2. Gypse revenu à 640 °C ; 3. et 4. Gypse revenu 
à 940 0C ; 5. et 6. Gypse revenu à 1 240 0C ; 1, 3, 5. Sans 
irradiation ; 2, 4, 6. Avec irradiation aux rayons X. 
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L’allure des courbes obtenues est différente suivant 
que l’élévation de température a été faite progres- 
sivement ou par palier, et dans ce dernier cas, suivant. 
l'importance des « sauts » de température (fig. 9). 


Nombre d'impulsions 
surface x temps 


0 4 10 024 30 LUF 50 60 70 
Temps t, minutes 


Fic. 9. — (D’après Haxel, Houtermans et Seeger [17]). 

Émission électronique d’un échantillon d’or passé à 
l’émeri, pour un saut brusque de la température. 

Douze minutes après le polissage, on a brutalement élevé 
la température, puis on l’a maintenue constante. Cinquante 
minutes après, on l’a élevé progressivement. 

La surface utilisée mesurait 10 cm?. 


L'’aire hachurée est sensiblement constante dans 
les 2 cas de la figure 9. Ce résultat est valable pour 
tous les corps. 


3) INFLUENCE DES TRAITEMENTS. — [La valeur 
absolue de l’émission et les hauteurs relatives des 
maxima en fonction de la température dépendent d’une 
manière très marquée des traitements subis par 
l’échantillon. Cependant l’allure de la décroissance et 
la position des maxima demeurent sensiblement iden- 
tiques. 
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Nombre de coups par seconde 


10 15 29 
Ftirege % 
F1c. 10. — (D’après Kramer [25]). 

Émission photoélectronique de l’aluminium étiré. 

1. Alsans préparation, sans filtre ; 2. Al sans préparation, 
avec filtre GG 13 ; 3. Al attaqué à la soude caustique et à 
l’acide fluorhydrique, sans filtre. 4. Al attaqué à la soude 
caustique et à l’acide fluorhydrique, avec filtre GG 13. 


Lu a) Influence d’une déformation. — Kramer [25] a 
étudié l'influence de l’étirage sur l'émission. La 
figure 10 montre que l’émission des photoélectrons 
“croît très vite avec l’étirage. 

On constate aussi que plus l’étirage est grand, plus 
l'émission photoélectrique décroît rapidement (fig. 11). 


AG 

ag) : 

0,2 | j Ï 

2 46810 20 #060100 3 5 10 20 50 100150 
Temps écoulé après J'étirage, min 


Nombre de coups par sec. 


Fic. 11. — (D’après Kramer [25]). 
. Décroissance des photoélectrons, après étirage de l’alumi- 
nium. 
» 1. Étirage à 1 %, sans filtre ; 2. Étirage à 2,5 %, sans 
filtre ; 3. Étirage à 3 %, avec filtre GG 13. 


b) {Influence du polissage. — Nous retiendrons les 
observations de Haxel, Houtermans et Seeger [17] qui 
ont constaté que la valeur du produit Zt (voir plus 

- haut) dépend d’une part de la taille des grains d’émeri 
… et d’autre part de la durée du polissage. Plus l’émeri 
“est fin, plus l’émission est importante. Haxel, 
— Houtermans et Seeger polissaient habituellement avec 
de l’émeri de 0,1 mm, pendant 30 secondes (saturation 
atteinte dans le cas du plomb). Ils ont dressé une liste 
des métaux suivant l'importance de l’émission après 
… un polissage à l’émeri fin : Pb, In, Al, Au, Cu, Mo, 
… Sn, Fe, alliage de Wood, graphite. 

c) Influence de l’irradiation. — L’irradiation est un 
traitement beaucoup plus souple que le polissage. On 
peut, comme l’a fait Kramer, [25] irradier par endroits 
… l'échantillon : on constate que les zones très irradiées 
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À Fic. 12. — (D’après Kramer [25]). 
Exploration par les exoélectrons, d’une lamelle de gypse, 

irradiée 10 minutes aux rayons X, distance 10 em. 
1. Aussitôt après le rayonnement ; 2. Cinq jours après. 
Ouverture de la fenêtre : 2 X 5 mm?,. 


sont les plus émissives. Les courbes de la figure 12 
… (exoélectrons) et 13 (photoélectrons) ont été obtenues 
- en déplaçant l'échantillon devant le diaphragme 

2 X 5 mm? du compteur. On constate aussi que 
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l'émission est beaucoup plus importante quand la 
longueur d’onde de l’éclairement est plus courte. Les 
courbes obtenues par Seeger [44], [45] illustrent ce 
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20 
Distance , mm 
Fic. 13. — (D’après Kramer [25]). 

Exploration par photoélectrons d’une lamelle de gypse, 
irradiée 30 secondes aux rayons X, distance 96 cm. Ouver- 
ture de la fenêtre : 2 X 5 mm?°,. 


phénomène (fig. 14). On remarquera que les maxima 
correspondent toujours aux mêmes valeurs de latempé- 
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F1. 14. — (D’après Seeger [44]). 
Post-courant d’un échantillon de molybdène après irra- 
diation par U. V. (courbe supérieure) ; après irradiation 
par lumière visible (courbe inférieure) (t — 0 fin de l’irra- 
diation). 


d) Influence d’un revenu. — Un traitement ther- 
mique, un revenu, modifient l’émission mais cependant 
sans déplacer la position des maxima (fig. 8). Sur la 
figure 15, nous voyons que l’émission du gypse, traité 
pendant 3 heures à 1 240 0C décroît très lentement ; 
par contre l'émission au temps { — 0 est beaucoup plus 
faible que dans le cas du gypse non traité. 


e) Cas d'un traitement par bombardement corpus- 
culaire. — Haxel, Houtermans et Seeger [17] ont 
étudié l’exoémission après une décharge gazeuse 
(0,4 mm Hg, 1 à 2 sec, 0,5.10—4 A). Les courbes 
obtenues ont sensiblement la même allure que dans 
le cas d’un traitement mécanique ou d’une irradiation. 
Il est intéressant de remarquer que les intensités sont 
à présent 10 à 100 fois plus grandes que dans le cas 
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du polissage et qu’une décharge dans l’hydrogène 
donne des émissions très importantes. 

Les résultats que nous venons d’exposer montrent 
à quel point l’émission des exo et des photoélectrons 
dépend des conditions dans lesquelles elle est produite. 
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Fic. 15. — (D’après Kramer [25]). 
Gypse irradié aux rayons X. Décroissance des exoélectrons. 
1. Sans préparation ; 2. Gypse revenu pendant 3 heures 
à 4 240 0C. 


4) CAS PARTICULIERS D’EXOÉMISSION. — On peut 
observer également une émission d’exoélectrons après 
une attaque chimique ou pendant un changement de 
phase (solidification). 


a) Après avoir attaqué un échantillon d'Al poly- 
cristallin avec le mélange HF + HCI dilué à 5 %, 
Bruna, Lintner, Müller et Schmidt [8] ont pu obtenir 
une émission d’exoélectrons (fig. 16) qui suit les lois 


Ù nn 
a © 
ao © 


à 
o 
[e] 


20 80 
Temps , minutes 


Nb. de particules par minute 


Fic. 16. — (D’après Bruna, Lintner, Müller et Schmidt [8]). 
Émission électronique d’un échantillon d’Al attaqué par 
HF + HCI] dilués. 


expérimentales que nous avons énoncées (décroissance 
suivant une loi hyperbolique, maxima en fonction 
de la température). Cette méthode particulière de 
production  d’exoélectrons semble apporter une 
confirmation à la théorie de Haxel, Houtermans et 
Seeger. 


b) Au cours d’une transformation allotropique [22], 
[32] — passage du chrome h. c. au chrome €. e. — ou 
au cours d’une solidification [8], [32], on observe 
également une émission d’exoélectrons. Bruna, Lintner, 
Müller et Schmidt ont réalisé un montage spécial pour 
l'étude du mercure (fig. 17) ; ils ont constaté que, au 
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moment de la solidification, le mercure et l’alliage… 
de Wood donnent lieu à une forte émission d’exo- 


F1G.°17. — (D’après Bruna, Lintner, Müller et Schmidt [8/). | 
Appareillage d’essai pour la mesure de l’émission exoélec 
trique lors de la solidification du mercure. | 


électrons (fig. 18). Mais dans le cas inverse — passage” 
de la phase solide à la phase liquide — aucune émission 
particulière n’est observée. 
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Fi1G. 18. — (D’après Bruna, Lintner, Müller et Schmidt [8]): 


Électroémission de l’alliage de Wood, lors de sa solidi=M 
fication. ; 


20 Étude de l’effet Russel. — Les exoélectrons ont” 
des énergies beaucoup trop faibles pour impressionner 
directement une plaque photographique (vitesse 
d'émission inférieure à 10 volts). En réalité le noircis- 
sement des plaques provient de l’action du peroxyde 
d'hydrogène H,0, produit après la sortie des exo- 
électrons. On a (d’après [15]). 


e + H,0 — H,0. — H. + (: OH)— 
H. +0, — HO; 
HO;. + H,0 > H,0/ 4+10H 
e + .OH — (: OH). 


La formation d’une image photographique est par 
conséquent liée à la présence de vapeur d’eau et 
d'oxygène : Grunberg et Wright [13] et Becker et 
Wolff [2] ont montré que la formation d’une image 
photographique (par contact) est impossible dans le 
vide. 

Les études de Grunberg et Wright ont porté sur des 
films métalliques vaporisés sous vide. Ils ont mesuré 


… No! 


la densité de l’image en fonction de la pression d'air 
(avec 1,8 %, de vapeur d’eau). Grunberg [15] a étudié 
… également la formation du peroxyde d'hydrogène dans 
l’eau. 


CHAPITRE III 


ASPECT THÉORIQUE 


L’émission exoélectronique est interprétée sur des 
bases différentes par Kramer; par Haxel, Houtermans 


— et Seeger ; et par Nassenstein. La théorie de Nassenstein 


couvre un plus grand nombre de phénomènes et permet 
en particulier d'expliquer correctement l’exoémission 
des corps non métalliques et l’émission des métaux 


| après une irradiation. 
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19 Théorie de Kramer [21 à 29]. — Pendant un 
travail mécanique, la structure superficielle du métal 


mn est perturbée et Kramer pense qu'il y a apparition 
… d’un état particulier du métal qu'il a appelé « non- 
— métallique ». Le passage de l’état « non-métallique » à 
» l’état métallique normal, la recristallisation super- 


ficielle, libèrent de l’énergie qui sert à extraire les 
électrons. Dans le cas de l’émission après un traitement 
par rayons X ou par bombardement électronique, on 
admet que l'énergie est fournie par le passage des 
particules dans la matière (voir par exemple [46|). 


20 Théorie de Haxel, Houtermans et Seeger [17]. — 
Haxel, Houtermans et Seeger ont fait remarquer que 
le changement de phase de la théorie précédente 
implique la libération d’une énergie anormalement 
élevée. Ils ont pensé au contraire que l’énergie pro- 
venait d’une réaction chimique à la surface du métal, 
en l’occurrence une oxydation. Le polissage ou la 
décharge gazeuse élimine la couche d’oxyde par 
endroits. Quand un atome d'oxygène se réunit à 
nouveau avec un atome métallique, une certaine 
uantité d'énergie est libérée et sert à extraire un 
électron.-En chaque point où la pellicule d'oxyde a été 


enlevée par un traitement mécanique, apparaît ainsi 
* un centre émetteur. 


Une énergie d’activation est nécessaire pour amor- 


- cer la réaction d’oxydation. Cette énergie est apportée 
par l'agitation thermique des atomes considérés ou 
- par un chauffage (apparition des maxima de l’émission 
—… qui correspondent aux différentes valeurs de l’énergie 


d'activation). Les calculs conduisent à des énergies 
d'activation en accord avec l’expérience et montrent 
aussi que. la densité des centres émetteurs a pour 
valeur 105 /cm? — donc environ un atome superficiel 


- sur 40° est un centre émetteur. 


Les expériences de Seeger et de Lohff et Raether 
semblent confirmer cette théorie. Seeger [45] a étudié 
lPémission exoélectrique du tungstène après un bom- 
bardement électronique (50 A /0,1 cm? ; 1 000 V); 


il a observé que l’on pouvait renouveler l’émission 
- d’un échantillon épuisé en le réoxydant. Lohff et 


Raether [33], [34] ont trouvé que plus les métaux 


+ sont oxydables, plus ils émettent : dans un vide 


de 2.105 Torr., Al, Mg, Pb, Na, Ca (fig. 7) émettent. 
Mais il n’en est pas de même pour Fe, Ni, Au, Ag, Pt ; 
pour ces derniers métaux, il faut augmenter la pression 
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partielle d'oxygène pour obtenir une émission. L’exo- 
émission semble donc provenir d’une réaction 
d’oxydation. 

La théorie de Haxel, Houtermans et Seeger est 
cependant incomplète car elle ne permet pas d’inter- 
préter l’émission des exoélectrons après une irra- 
diation X ou U. V. (à moins d'admettre que le rayon- 
nement « arrache » des atomes d'oxygène). 


39 Théorie de Nassenstein [38]. Nassenstein 
interprète l’émission des exoélectrons d’une manière 
totalement différente. Il fait appel à la théorie des 
centres F [36] et de la phosphorescence [401, et compare 
en particulier la décroissance de l’exoémission à la 
décroissance de la phosphorescence. Ce parallélisme 
entre l’exoémission et la phosphorescence avait déjà 
été souligné par Kramer [27]. 

L'hypothèse fondamentale de cette théorie est la 
suivante : les impuretés et les défauts du réseau créent 
des niveaux énergétiques relativement élevés. Ces 
niveaux (ou pièges), sur lesquels l’électron n’est pas 
libre, sont remplis par une irradiation ou par un 
travail mécanique [30], [39]. En certains points (ou 
centres) de la surface le travail de sortie ® est abaïssé 
et prend une valeur plus faible € ; l’électron a donc 
une probabilité plus grande de s'évader. Deux cas 
peuvent se présenter alors : 


a) A la température ambiante, l’agitation thermique 
communique à l’électron suffisamment d'énergie pour 
qu’il puisse franchir « spontanément » la hauteur du 
potentiel e. En ce cas apparaït l’émission exoélectrique. 


b) La hauteur € reste malgré tout importante et il 
faut chauffer légèrement l’échantillon pour obtenir 
une émMISSION. 

En ce qui concerne l’émission photoélectrique, nous 
voyons que le seuil est à présent décalé vers les grandes 
longueurs d’onde : 


L'émission est nulle quand tous les centres sont 
épuisés. Une nouvelle irradiation ou un nouveau 
chauffage pourra les remplir à nouveau. Cette possi- 
bilité de renouveler l’émission disparaît après quelques 
opérations quand les défauts du réseau ont disparu 
ou quand des couches absorbantes se sont formées à 
la surface de l’échantillon. 


DÉCROISSANCE EN FONCTION DU TEMPS. 
Nassenstein étudie en détail l’émission exoélectronique 
des corps métalliques et non métalliques en tenant 
compte de la distribution énergétique des centres 
émetteurs.* Soit : J le nombre d'électrons émis par 
seconde ; € le niveau énergétique d’un centre en 
eV ; t le temps en secondes ; T la température abso- 
lue ; Æ la constante de Boltzmann. 

L’électron a une certaine probabilité de franchir la 
hauteur de potentiel £ : cette probabilité est donnée 


par la formule p = s exp (x) avéc s #10+8 sec. 


Voir [36]. 
Il est Rors nécessaire de faire une hypothèse sur 
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le spectre de distribution énergétique et Nassenstein 
‘distingue quatre cas : 

4) Spectre de valeurs discrètes : Soit n le nombre 
d'électrons dans les centres de niveau e, et seulement e, 
et soit », le nombre au temps t = 0. On obtient alors : 


TJ=—"n _ = nos € ÀT exp [— st e—€/kT] 


ou encore 
J = Jo exp [— stet/#T]. 


Si les électrons proviennent directement des centres’ 
n = 1 et n, est égal au nombre de centres se trouvant 
suffisamment près de la surface pour permettre à leurs 
électrons d'abandonner directement le cristal. Mais 
si les électrons passent tout d’abord dans la bande de 
conduction, il faut faire intervenir des processus de 
« retrapping » et ce cas se rattache à la théorie des 
cathodes à oxyde [18]. On suppose ici n = 1 et l’on a 
alors la figure 19. 


Ets 


Fic. 19. — (D’après Nassenstein [38]). 
Décroissance de l’émission pour une seule valeur discrète 


du spectre énergétique des centres : 
Et=10,00 6 TI =S000K 


Au lieu d’une seule valeur discrète €, on peut en 
avoir plusieurs. On obtient alors la figure 20 que l’on 
comparera à la courbe expérimentale de la figure 21. 


2) Spectre uniforme : Soit N, le nombre de centres 
au temps & — 0, de niveau énergétique compris entre € 
et + de. On obtient alors : 


= IN Ne s e—E/kT exp [— st.e—#/kT] de. 


N, étant constant dans tout le domaine de variation 
de €, on obtient : 


J = Na.(KT ft).(1 — est). 


Pour t > 1077 sec, cette expression devient prati- 
quement 
J _ Ne RT 
t 


et l’on obtient la loi souvent observée en 
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Nous 
3) Spectre exponentiel : On a : 

Ne — Aer 
FE [  Ais.exp [— ste #T].exp [— e/kT — aejde. 


L'intégration donne 
J = js AT). BtCRT +1) 
c’est-à-dire on obtient la loi en 4 avec 
a = «kT +1. 


F1G. 20. — (D’après Nassenstein [38]). 
Décroissance de l’émission à 300 °K, pour deux valeurs 
discrètes du spectre : &, — 0,55 eV €, — 0,65 eV. 
Courbe a :n, — n; — 104 (ordonnée de gauche). 
Courbe b : n; — 10%, n; — 105 (ordonnée de gauche). 
Courbe c : n; — 105, n,; — 104 (ordonnée de droite). 
Ordonnée : Nombre d’électrons émis par seconde. 


10° sec 


t — 


Fic. 21. — (D’après Nassenstein [38]). 
Électroémission de SO,K, après concassage du cristal. 
Ordonnée J : Nombre d’électrons émis pour 5 secondes. 
Abscisse : Temps t après le broyage. 
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4) Spectre en créneaux. On a alors : L’allure des courbes dépend très fortement de € 
KT et es comme l'indique la figure 22. Remarquons 

EE Er LE —€a/KT| DE y—E1/T É 1 7 ; 2 
“OS t jexp DR également l'influence de la température. On a supposé 
puisque V, = constante pour & < € < € au début l'existence d’un seul créneau, mais on peut 
Ne = Opoure <e, et pour e > ee. : en considérer plusieurs, comme dans les figures 23, 24, 


25. Nous voyons combien l’allure des courbes dépend 
de la répartition de l’énergie, ce qui permet la recherche 
d’une « structure fine ». 


Fic. 22. — (D’après Nassenstein[38]). 5 40 
Décroissance de l’émission pour les répartitions éner- 1 10 10 
gétiques suivantes : ere 
Ne — constante pour €, <e<e,; Ne = 0 pour € < €; Fic. 24. — (D’après Nassenstein [38]). 
et pour € > e. Décroissance de l’émission à 300 °K pour le cas:: 
Courbe a : €, — 0,55eV, €, — 0,75 eV,. T — 300 °K. Me = 1,6.105 (eV)—1 pour 0,58 < € < 0,60 eV. 
Courbe b: e — €; > 0,4 eV, T — 300 °K. Ne = 1,6.106 (eV)—1 pour 0,63 < € < 0,65 eV. 
Courbelc : €, — 0,55 eV, . €, — 0,75 eV, T = 350 °K. Ne = 1,6.107 (eV)—1 pour 0,68 < € < 0,70 eV. 
Courbe d': e — €, > 0,4eV, T = 350 0K. Ordonnée : Nombre d’électrons émis par seconde. 


10 sec 


10° 


10° 


10° 


+ — 10°! 
Fc. 23. — (D’après Nassenstein [38]). 1 10° 10 & 10 sec 10° 
Décroissance de l’émission pour les distributions éner- 
gétiques suivantes : Fic. 25. — (D’après Nassenstein [38]). 
= —=1,6.106 (eV)—1 pour 0,58 < € < 0,60 eV Décroissance de l’émission à 300 °K pourle cas : 
0,36 < e < 0,65 eV Ne = 1,6.107 (eV)—2 pour 0,58< € < 0,60 eV. 
0,68 < € < 0,70 eV Ne = 1,6.106 (eV)—1 pour 0,63< € < 0,65 eV. 
T = 3000K Ne — 1,6.105 (eV)—1 pour 0,68 < e < 0,70 eV. 


Ordonnée : Nombre d'électrons émis par seconde. Ordonnée : Nombre d’électrons émis par seconde. 
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ÉMISSION EN FONCTION DE LA TEMPÉRATURE. — 
Nassenstein suit ici Randall et Wilkins [40]. Soit 6 la 
vitesse de l’élévation de température : d7 = Pdi. 

On a alors, en supposant une seule valeur discrète 
du niveau énergétique : 


T 
J = nn(T).exp EJ Ses ar] .s.e78/kT, 


Lors d’une élévation de température, l'émission croît 
exponentiellement dans un certain domaine de tempé- 
rature, fonction de €, atteint un maximum, décroit et 
s’annule quand tous les centres sont vides. La tempé- 
rature 7» du maximum dépend étroitement de la 
vitesse 8. Pour B voisin de 1° par seconde on a : 
er 25 kTn c’est-à-dire 50 °K correspondent à 
environ 0,1 eV. 

Si la distribution des pièges est plus compliquée, il 
apparaît plusieurs maxima. L’analyse de la courbe 
donne donc directement la répartition énergétique des 
centres actifs pour l’émission. Toutefois il ne faut pas 
s’attendre à une concordance des énergies d’activation 
optiques €, avec les énergies d’activation thermique er, 
du fait du principe de Franck-Condon [36]. 


VÉRIFICATION DE LA THÉORIE DE NASSENSTEIN. — 
Nassenstein reconnait que les considérations mathé- 
matiques indiquées plus haut ne doivent pas être 
considérées comme une preuve de la théorie. Mais il 
fait remarquer que les théories basées sur un chan- 
gement de phase ou sur une réaction chimique super- 
ficielle peuvent très difficilement interpréter certains 
faits expérimentaux, en particulier pour les corps non 
métalliques (?). 

On a déjà vu par exemple le décalage du seuil 
photoélectrique. Kramer [25] a également observé le 
phénomène suivant : après l’excitation optique d’un 
corps non métallique, les maxima de l’émission par 
chauffage de l’échantillon ont lieu à des températures 
d'autant plus basses que l’on emploie des longueurs 
d’ondes plus courtes. Ceci est compréhensible car plus 
courtes sont les longueurs d’ondes, plus grands sont les 
quanta d’énergie, donc plus élevés les niveaux occupés, 
et par conséquent plus basses les températures suscep- 
tibles de vider ces niveaux. On comprend également 
cet autre fait expérimental : les substances ayant une 
émission exoélectronique à la température ambiante 
possèdent toujours une émission photoélectrique, mais 
l'inverse ne se produit pas. 


TRAVAUX DE BoHuN. — Dans quelques cas on peut 
même obtenir la preuve expérimentale que l’exo- 
émission est liée d’une façon fondamentale à la présence 
de centres F. Il faut ici se tourner vers les travaux 
récents de Bohun [3], [4], [6] en particulier sur le 
chlorure de sodium NaCl [5]. Les résultats de Bohun 
peuvent se résumer ainsi : 


1) Un cristal de NaCl non coloré ne présente ni 
photoémission par illumination avec de la lumière 
visible, ni thermoémission. 


2) Après coloration du cristal par irradiation aux 
rayons X, on obtient par chauffage de nombreux 


(2) Pour les métaux, la formation d’une couche d’oxyde 
ne semble pas incompatible avec la théorie de Nassenstein. 
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Not 


maxima ; par illumination on obtient aussi une photo- 


émission. 


3) Dès que le cristal change de couleur du fait du : 
chauffage ou de l’illumination, l’émission disparaît et « 


n'apparaît plus même par réchauffement. Cet effet peut 
se répéter à volonté. 

Bohun a poursuivi ces recherches ; son dernier 
travail s’est effectué sur le cuivre [7]. Il faut remarquer 
qu’il tient Cu:O et CuO comme cause des pièges 
responsables de l’émission. 

Nassenstein a retrouvé les résultats de Bohun sur le 


chlorure de sodium NaCI. Des travaux encore plus ! 


récents de Seeger [45] parlent en faveur de la théorie 
des niveaux. D’autres résultats parus au même moment 


et dus à Lohff et Raether [33] sont interprétés par ces . 


auteurs dans le cadre de la théorie chimique. Mais ils 
peuvent, d’après Nassenstein, s’interpréter également 
par la théorie des niveaux énergétiques. 


PRÉCISIONS SUR LE MÉCANISME D’'ÉMISSION. 
Bohun suppose que l’émission est due à des chocs de 
seconde espèce entre des excitons et des centres F 
excités thermiquement, chocs ayant lieu à la surface 
du cristal. Cette hypothèse semble étayée par des 
travaux théoriques non publiés dus à Ansbacher [43] ; 
d’après ce dernier, lors d’un choc inélastique entre un 
électron excité provenant d’un- défaut et un électron 
de conduction, l’énergie d’excitation peut se commu- 
niquer à l’électron de conduction sous forme d’énergie 
cinétique, tandis que l’électron provenant du défaut 
retourne à l’état fondamental. 

D’autres recherches sont nécessaires, en particulier 
il, serait intéressant d'entreprendre des recherches 
parallèles avec d’autres méthodes, par exemple des 
mesures de phosphorescence et de conductibilité. Un 
premier travail dans cette direction a été publié par 
Lepper [31 ]sur le phosphore, où il étudie simultanément 
l'émission d’électrons et l’émission de lumière, en 
particulier lors d’un chauffage. Lepper admet le 
mécanisme suivant : les électrons sont tout d’abord 
élevés des pièges jusqu’à la bande de conduction. Une 
partie des électrons, se trouvant dans le voisinage de 
la surface, diffuse dans une bande superficielle, d’où 
les électrons sortent par un nouvel apport d’énergie 
thermique. Si la sortie des électrons se produit direc- 
tement en sautant la bande de conduction, le travail 
de sortie joue alors nécessairement un grand rôle : ne 


peuvent alors émettre que les cristaux dont le travail 


de sortie est suffisamment faible. Il est toutefois 
intéressant de noter que Bohun a trouvé de nombreux 
maxima de l’émission par chauffage même pour le 


chlorure d’argent AgCl, dont le travail de sortie est. 


égal à 3,5 eV. 


CONCLUSION 


L'émission retardée lumineuse et l’émission retardée 
électronique présentent de nombreuses analogies : non 
seulement ces émissions obéissent aux mêmes lois 
(décroissance en fonction du temps, influence de la 
température), mais encore sont provoquées par les 
mêmes traitements. 

Ce parallélisme entre l’exoémission et la phos- 
phorescence [27] parle en faveur de la théorie « élec- 
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_ tronique » de Nassenstein, beaucoup plus générale que d’ailleurs s’étendre, d’après Nassenstein, à la catalyse 
M la théorie «chimique » de Haxel, Houtermans et Seeger. hétérogène et à la métallurgie des poudres [9], [16]. 
D'LA théorie de Nassenstein en effet s’applique indif- Remarquons enfin que cette méthode d'étude des 
“iéremment aux métaux ou aux corps non métalliques, corps solides, dont les travaux de Kramer sont à 
explique facilement certains phénomènes et permet une l’origine, permet d’étudier le grand nombre de maté- 
recherche de la structure fine de la distribution éner-  riaux qui ne sont pas phosphorescents. 
gétique de la surface des corps solides. 


Le domaine d’application de cette théorie pourrait Manuscrit reçu le 12 juillet 1956. 
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LETTRE A LA RÉDACTION 


SPECTROGRAPHE A VIDE 
ET RÉSEAU CONCAVE 
EN INCIDENCE NORMALE 
POUR L’ULTRAVIOLET LOINTAIN 


Par Mme Simone Rogin et M. Stéphane Roi, 
Institut des Hautes Études de Dakar. 


Ce spectrographe prévu pour un réseau Bausch et 
Lomb de 1 mètre de rayon de courbure et portant 
1 200 traits au mm (:), a été dessinéen collaboration avec 
M. B. Vodar et le Laboratoire des Hautres Pressions 
de Bellevue et réalisé par les Etablissements Beaudoin 
(Paris). Il bénéficie des perfectionnements techniques 
apportés à divers appareils de ce genre construits et 
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La fente réglable de l'extérieur (K) et la vanne àn 
lumière (V) (étanche au vide dans les 2 sens) sont d’un 
type précédemment décrit [3] mais ils ont été réunis 
en un bloc (B) fixé à l’extérieur du spectrographe. Par} 
suite de cette disposition, la fente se trouve à l’exté| 
rieur du cercle de Rowland mais on sait que, si lai 
plaque est convenablement disposée (légèrement vers 
lintérieur du cercle de Rowland), les qualités du! 
spectre ne sont pratiquement pas diminuées en inéi=M{ 
dence normale [4]. Ce bloc, très facilement démontable,k | 
est d’ailleurs réversible, ce qui permet de placer la 
fente en arrière de la vanne à lumière (pour être plus" 
proche du cercle de Rowland) ou en avant pour faciliter 
certains montages, notamment l’utilisation des sources | 
sans fenêtres [2]. di 

Ce spectrographe est actuellement équipé d’un 
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mis au point au Laboratoire de Physique-Enseignement 
à la Sorbonne [1], [2], [3]. Il est entièrement réalisé en 
acier pour éviter toute porosité et utilise exclusivement 
des joints toriques pour les assemblages, fermetures et 
entrées de mouvement ; le problème de la recherche 
des fuites s’est trouvé, par suite, presque totalement 
éliminé et l’appareil a pu être mis en service en quelques 
ours. 

La tête du spectrographe porte un orifice de pom- 
page (P) le châssis (C), la fente et les entrées de mouve- 
ment (E) ; le porte-réseau (R) est fixé à l’intérieur d’un 
berceau fait d’un tube d’acier (T). fixé à la tête du 
spectrographe et sur lequel vient se placer une cloche 
horizontale (D) mobile sur rails ; cette disposition 
permet un accès facile à l’intérieur de l’appareil ; des 
tubulures (G et H) sont prévues pour l'installation des 
jauges de mesure du vide. Le réglage est assuré par 
rotation et translation du réseau (E,, E,) et par rota- 
tion du châssis (E,) qui sont commandés de l’extérieur 
ainsi que le déplacement d’un cache (E,) et un système 
de translation (E;) du châssis pour la prise des spectres ; 
le châssis pour plaques 9 X 12 permet la prise de 
24 spectres de 2 mm de largeur. 


(*) Ce réseau concentre la lumière aux environs de 1 700 
dans le 1er ordre. 
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Le 
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groupe de pompage Edwards comprenant une pompe 
primaire à palettes 1S 50 et une pompe secondaire à 
huile F 203 ; ce groupe permet d'obtenir un vide 
inférieur à 104 mm de Hg en 15 minutes environ. 
Pour permettre son utilisation à des mesures d'absOrp-« 
tion, il est équipé en outre d’un dispositif de tourelles 
classiques (cf. [1], [3]) renfermant 2 lentilles de fluorine 
qui focalisent la lumière de la source sur la fente après 
passage à travers l’absorbant ou 5 grilles d'étalonnage 
dont le déplacement est commandé de l’extérieur. Avec 
les lentilles de fluorine dont nous disposons (?), l’appa- 
reil en service depuis quelques mois nous à permis 
d'enregistrer des spectres s'étendant jusqu’à 1 400 À À 
environ. "+ 
(2) Cette fluorine provient de la Société de Recherches et. f 
d'Exploitation minière (Pau). \ 
Manuscrit reçu le 7 août 1956, 
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